Measuring the variability of local characteristics in complex networks: Empirical and analytical analysis

We examine the dynamics for the average degree of a node's neighbors in complex networks. It is a Markov stochastic process, and at each moment of time, this quantity takes on its values in accordance with some probability distribution. We are interested in some characteristics of this distribution: its expectation and its variance, as well as its coefficient of variation. First, we look at several real communities to understand how these values change over time in social networks. The empirical analysis of the behavior of these quantities for real networks shows that the coefficient of variation remains at high level as the network grows. This means that the standard deviation and the mean degree of the neighbors are comparable. Then, we examine the evolution of these three quantities over time for networks obtained as simulations of one of the well-known varieties of the Barabási-Albert model, the growth model with nonlinear preferential attachment (NPA) and a fixed number of attached links at each iteration. We analytically show that the coefficient of variation for the average degree of a node's neighbors tends to zero in such networks (albeit very slowly). Thus, we establish that the behavior of the average degree of neighbors in Barabási-Albert networks differs from its behavior in real networks. In this regard, we propose a model based on the NPA mechanism with the rule of random number of edges added at each iteration in which the dynamics of the average degree of neighbors is comparable to its dynamics in real networks.

tatDB: a database of Ago1-mediated targets of transfer RNA fragments

tRNA-derived fragments (tRFs) are a class of emerging post-transcriptional regulators of gene expression likely binding to the transcripts of target genes. However, only a few tRFs targets have been experimentally validated, making it hard to extrapolate the functions or binding mechanisms of tRFs. The paucity of resources supporting the identification of the targets of tRFs creates a bottleneck in the fast-developing field. We have previously analyzed chimeric reads in crosslinked Argonaute1-RNA complexes to help infer the guide-target pairs and binding mechanisms of multiple tRFs based on experimental data in human HEK293 cells. To efficiently disseminate these results to the research community, we designed a web-based database tatDB (targets of tRFs DataBase) populated with close to 250 000 experimentally determined guide-target pairs with ∼23 000 tRF isoforms. tatDB has a user-friendly interface with flexible query options/filters allowing one to obtain comprehensive information on given tRFs (or targets). Modes of interactions are supported by secondary structures of potential guide-target hybrids and binding motifs, essential for understanding the targeting mechanisms of tRFs. Further, we illustrate the value of the database on an example of hypothesis-building for a tRFs potentially involved in the lifecycle of the SARS-CoV-2 virus. tatDB is freely accessible at https://grigoriev-lab.camden.rutgers.edu/tatdb.

SARSNTdb database: Factors affecting SARS-CoV-2 sequence conservation

SARSNTdb offers a curated, nucleotide-centric database for users of varying levels of SARS-CoV-2 knowledge. Its user-friendly interface enables querying coding regions and coordinate intervals to find out the various functional and selective constraints that act upon the corresponding nucleotides and amino acids. Users can easily obtain information about viral genes and proteins, functional domains, repeats, secondary structure formation, intragenomic interactions, and mutation prevalence. Currently, many databases are focused on the phylogeny and amino acid substitutions, mainly in the spike protein. We took a novel, more nucleotide-focused approach as RNA does more than just code for proteins and many insights can be gleaned from its study. For example, RNA-targeted drug therapies for SARS-CoV-2 are currently being developed and it is essential to understand the features only visible at that level. This database enables the user to identify regions that are more prone to forming secondary structures that drugs can target. SARSNTdb also provides illustrative mutation data from a subset of ~25,000 patient samples with a reliable read coverage across the whole genome (from different locations and time points in the pandemic. Finally, the database allows for comparing SARS-CoV-2 and SARS-CoV domains and sequences. SARSNTdb can serve the research community by being a curated repository for information that gives a jump start to analyze a mutation’s effect far beyond just determining synonymous/non-synonymous substitutions in protein sequences.

COVID-19 and neurodegeneration: The mitochondrial connection

There is still a significant lack of knowledge regarding many aspects of the etiopathology and consequences of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in humans. For example, the variety of molecular mechanisms mediating this infection, and the long-term consequences of the disease remain poorly understood. It first seemed like the SARS-CoV-2 infection primarily caused a serious respiratory syndrome. However, over the last years, an increasing number of studies also pointed towards the damaging effects of this infection has on the central nervous system (CNS). In fact, evidence suggests a possible disruption of the blood-brain barrier and deleterious effects on the CNS, especially in patients who already suffer from other pathologies, such as neurodegenerative disorders. The molecular mechanisms behind these effects on the CNS could involve the dysregulation of mitochondrial physiology, a well-known early marker of neurodegeneration and a hallmark of aging. Moreover, mitochondria are involved in the activation of the inflammatory response, which has also been broadly described in the CNS in COVID-19. Here, we critically review the current bibliography regarding the presence of neurodegenerative symptoms in COVID-19 patients, with a special emphasis on the mitochondrial mechanisms of these disorders.

Large-Scale Computational Discovery of Binding Motifs in tRNA Fragments

Accumulating evidence has suggested that tRNA-derived fragments (tRFs) could be loaded to Argonaute proteins and function as regulatory small RNAs. However, their mode of action remains largely unknown, and investigations of their binding mechanisms have been limited, revealing little more than microRNA-like seed regions in a handful of tRFs and a few targets. Here, we identified such regions of potential interaction on a larger scale, using in vivo formed hybrids of guides and targets in crosslinked chimeric reads in two orientations. We considered "forward pairs" (with guides located on the 5' ends and targets on the 3' ends of hybrids) and "reverse pairs" (opposite orientation) and compared them as independent sets of biological constructs. We observed intriguing differences between the two chimera orientations, including the paucity of tRNA halves and abundance of polyT-containing targets in forward pairs. We found a total of 197 quality-ranked motifs supported by ∼120,000 tRF-mRNA chimeras, with 103 interacting motifs common in forward and reverse pairs. By analyzing T→C conversions in human and mouse PAR-CLIP datasets, we detected Argonaute crosslinking sites in tRFs, conserved across species. We proposed a novel model connecting the formation of asymmetric pairs in two sets to the potential binding mechanisms of tRFs, involving the identified interaction motifs and crosslinking sites to Argonaute proteins. Our results suggest the way forward for further experimental elucidation of tRF-binding mechanisms.

Computational meta-analysis of ribosomal RNA fragments: potential targets and interaction mechanisms

The most abundant cellular RNA species, ribosomal RNA (rRNA), appears to be a source of massive amounts of non-randomly generated fragments. We found rRNA fragments (rRFs) in immunoprecipitated Argonaute (Ago-IP) complexes in human and mouse cells and in small RNA sequencing datasets. In human Ago1-IP, guanine-rich rRFs were preferentially cut in single-stranded regions of mature rRNAs between pyrimidines and adenosine, and non-randomly paired with cellular transcripts in crosslinked chimeras. Numerous identical rRFs were found in the cytoplasm and nucleus in mouse Ago2-IP. We report specific interaction motifs enriched in rRF-target pairs. Locations of such motifs on rRFs were compatible with the Ago structural features and patterns of the Ago-RNA crosslinking in both species. Strikingly, many of these motifs may bind to double-stranded regions on target RNAs, suggesting a potential pathway for regulating translation by unwinding mRNAs. Occurring on either end of rRFs and matching intronic, untranslated or coding regions in targets, such interaction sites extend the concept of microRNA seed regions. Targeting both borders of certain short introns, rRFs may be involved in their biogenesis or function, facilitated by Ago. Frequently dismissed as noise, rRFs are poised to greatly enrich the known functional spectrum of small RNA regulation.

Age-Related Argonaute Loading of Ribosomal RNA Fragments

BACKGROUND

Accumulating evidence points to the functional roles of rRNA derived Fragments (rRFs), often considered degradation byproducts. Small RNAs, including miRNAs and tRNA-derived Fragments (tRFs), have been implicated in the aging process and we considered rRFs in this context.

OBJECTIVE

We performed a computational analysis of Argonaute-loaded rRFs in Drosophila melanogaster to study rRF changes with age. We determined rRF abundance in Ago1 and Ago2 at 3 and 30 days to identify Ago1-guided and Ago2-guided fragments. We searched for putative seed sequences in rRFs based on frequent matches of sliding k-mer windows to the conserved regions of 12 Drosophila genomes. We predicted putative targets (containing matches to seeds identified in four rRFs) and studied their functional enrichments using Gene Ontology.

RESULTS

We identified precise cleavage sites of distinct rRF isoforms from both nuclear and mitochondrial rRNAs. The most prominent rRF isoforms were enriched in Ago2 at 3 days and that loading strongly decreased with age. For less abundant rRFs, loading of Ago2-guided rRFs generally increased in Ago2, whereas Ago1-guided rRFs revealed diverse patterns. The distribution of seed matches in targets suggested that rRFs may bind to various conserved regions of many genes, possibly via miRNA-like seed-based mechanisms.

CONCLUSION

Our observations suggest that rRFs may be functional molecules, with age-dependent Argonaute loading, comparable to that of miRNAs and tRFs. The putative rRF targets showed significant enrichment in developmental processes and biological regulation, similar to tRFs and consistent with a possible involvement of these newly identified small RNAs in the Drosophila aging.

Circulating Plasma microRNA to Differentiate Cushing's Disease From Ectopic ACTH Syndrome

Corticotropinomas and adrenocorticotropic hormone (ACTH)-secreting neuroendocrine tumors exhibit differential levels of some microRNAs (miRs) compared to normal tissue. Because miRs can be released from tissues into circulation, they offer promise as novel disease biomarkers. Objective: To evaluate whether miRs are differentially detected in plasma samples of patients with ACTH-dependent Cushing's syndrome (CS). Design: Case-control study. Methods: Morning fasting plasma samples were collected from 41 consecutive patients with confirmed ACTH-dependent CS and 11 healthy subjects and stored at -80°C. Twenty-one miRs previously reported to be differentially expressed in ACTH-secreting tumors vs. healthy tissue samples were quantified in plasma by qPCR. Results: Among enrolled subjects, 28 were confirmed to have Cushing's disease (CD), 13 had ectopic ACTH secretion (EAS) and 11 were healthy controls. We found statistically significant differences in the circulating levels of miR-16-5p [45.04 (95% CI 28.77-61.31) in CD vs. 5.26 (2.65-7.87) in EAS, P < 0.001; q = 0.001], miR-145-5p [0.097 (0.027-0.167) in CD vs. undetectable levels in EAS, P = 0.008; q = 0.087] and differences in miR-7g-5p [1.842 (1.283-2.400) in CD vs. 0.847 (0.187-1.507) in EAS, P = 0.02; q = 0.14]. The area under the receiver-operator (ROC) curve was 0.879 (95% CI 0.770-0.987), p < 0.001, when using miR-16-5p to distinguish between CD and EAS. Circulating levels of miR-16-5p in the healthy control group differed from that of both the CD and EAS groups. Conclusions: Plasma miR levels differ in patients with CD and EAS. In particular, miR-16-5p, miR-145-5p and miR-7g-5p are promising biomarkers for further research to differentiate ACTH-dependent CS.

Inferring targeting modes of Argonaute-loaded tRNA fragments

Transfer RNA fragments (tRFs) are an emerging class of small RNA molecules derived from mature or precursor tRNAs. They are found across a wide range of organisms and tissues, in small RNA fraction or loaded to Argonaute in numbers comparable to microRNAs. Their functions and mechanisms of action are largely unknown, and results obtained on individual tRFs are often hard to generalize. Here we predicted binding mechanisms and specific target interaction sites of 26 human Argonaute-loaded tRFs of different types using large-scale meta-analyses of available experimental data. Strikingly, our findings matched all interaction sites detected in a recent experimental screen, confirming the validity of our computational approach. Such sites are primarily located on the 5' end of tRFs and often involve additional binding along the tRF length, similar to microRNAs. Indicative of multiple layers of regulation, diverse regulatory non-coding RNAs comprised a third of the tRF targets, with the rest being protein-coding transcripts. In the latter, coding sequence and 3' UTRs were the likely primary target regions, although we observed interactions of tRFs with 5' UTRs. Another novel phenomenon we report, a large number of putative interactions between tRFs and introns, is compatible with the roles of Argonaute in the nucleus. Further, observed tRF-intron binding modes suggest a mechanism of interaction of tRFs with Argonaute-dependent introns, and we predict here >20 candidate introns of this type. Taken together, these results present tRFs as regulatory molecules with a rich functional spectrum.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Anticarcinogenic activity of blue fluorescent hexagonal boron nitride quantum dots: as an effective enhancer for DNA cleavage activity of anticancer drug doxorubicin

Blue fluorescent hexagonal boron nitride quantum dots (h-BNQDs) of ∼10 nm size as an effective enhancer for DNA cleavage activity of anticancer drug doxorubicin (DOX) were synthesized using simple one-step hydrothermal disintegration of exfoliated hexagonal boron nitride at very low temperature ∼ 120 °C. Boron nitride quantum dots (BNQDs) at a concentration of 25 μg/ml enhanced DNA cleavage activity of DOX up to 70% as checked by converting supercoiled fragment into nicked circular PBR322 DNA. The interaction of BNQDs with DOX is proportional to the concentration of BNQDs, with binding constant K b ∼0.07338 μg/ml. In addition, ab initio theoretical results indicate that DOX is absorbed on BNQDs at the N-terminated edge with binding energy -1.075 eV and prevented the normal replication mechanisms in DNA. BNQDs have been shown to kill the breast cancer cell MCF-7 extensively as compared with the normal human keratinocyte cell HaCaT. The cytotoxicity of BNQDs may be correlated with reduced reactive oxygen species level and increased apoptosis in MCF-7 cells, which may be liable to enhance the anticancerous activity of DOX. The results provide a base to develop BNQD-DOX as a more effective anticancer drug.

ARIADNA: machine learning method for ancient DNA variant discovery

Ancient DNA (aDNA) studies often rely on standard methods of mutation calling, optimized for high-quality contemporary DNA but not for excessive contamination, time- or environment-related damage of aDNA. In the absence of validated datasets and despite showing extreme sensitivity to aDNA quality, these methods have been used in many published studies, sometimes with additions of arbitrary filters or modifications, designed to overcome aDNA degradation and contamination problems. The general lack of best practices for aDNA mutation calling may lead to inaccurate results. To address these problems, we present ARIADNA (ARtificial Intelligence for Ancient DNA), a novel approach based on machine learning techniques, using specific aDNA characteristics as features to yield improved mutation calls. In our comparisons of variant callers across several ancient genomes, ARIADNA consistently detected higher-quality genome variants with fast runtimes, while reducing the false positive rate compared with other approaches.

Lightning-fast genome variant detection with GROM

Current human whole genome sequencing projects produce massive amounts of data, often creating significant computational challenges. Different approaches have been developed for each type of genome variant and method of its detection, necessitating users to run multiple algorithms to find variants. We present Genome Rearrangement OmniMapper (GROM), a novel comprehensive variant detection algorithm accepting aligned read files as input and finding SNVs, indels, structural variants (SVs), and copy number variants (CNVs). We show that GROM outperforms state-of-the-art methods on 7 validated benchmarks using 2 whole genome sequencing (WGS) data sets. Additionally, GROM boasts lightning-fast run times, analyzing a 50× WGS human data set (NA12878) on commonly available computer hardware in 11 minutes, more than an order of magnitude (up to 72 times) faster than tools detecting a similar range of variants. Addressing the needs of big data analysis, GROM combines in 1 algorithm SNV, indel, SV, and CNV detection, providing superior speed, sensitivity, and precision. GROM is also able to detect CNVs, SNVs, and indels in non-paired-read WGS libraries, as well as SNVs and indels in whole exome or RNA sequencing data sets.

Radiation situation dynamics at the Andreeva Bay site for temporary storage of spent nuclear fuel and radioactive waste over the period 2002-2016

The Coastal Technical Base (CTB) №569 at Andreeva Bay was established in the early 1960s and intended for the refueling of nuclear submarine reactors and temporary storage of spent nuclear fuel (SNF) and radioactive waste (RW). In 2001, the base was transferred to the Russian Ministry for Atomic Energy and the site remediation began. The paper describes in detail the radiation situation change at the technical site in Andreeva Bay from 2002-2016, the period of preparation for the most critical phase of remedial work: removal of spent fuel assemblies. The analysis of aggregated indicators and data mining were used. The article suggests the best number and location of checkpoints needed to ensure sufficient accuracy of the radiation situation description. The fractal properties of the radiation field are studied using the Hurst index. The relationship between checkpoints was assessed using the method of searching for checkpoint communities. The decrease in the integral of the ambient dose equivalent rate (ADER) at the technical site was evaluated by the method of time series decomposition. Three components of time series were identified: trend, seasonal and residual. The trend of the ADER integral over the technical site is a monotonic decreasing function, where the initial and final values differ tenfold. Taking into account that ¹³⁷Cs dominates the radiation situation on-site, it is clear that the ADER due to the radionuclide decay will have decreased by 1.4 times. It is estimated that only a small proportion of ¹³⁷Cs has migrated off-site. Therefore, approximately a sevenfold decrease in dose rate is mainly due to remediation activities of personnel. During the year, the seasonal component varies the ADER integral by a factor of two, due to snowfall. The residual component reflects the uncertainty of the ADER integral calculation and phases of active SNF and RW management. The methods developed are used to support the optimization of remediation work as well as regulatory supervision of occupational radiation protection.

Effects of active acromegaly on bone mRNA and microRNA expression patterns

OBJECTIVE

To evaluate the response of bone to chronic long-term growth hormone (GH) and insulin-like growth factor-1 (IGF1) excess by measuring the expression of selected mRNA and microRNA (miR) in bone tissue samples of patients with active acromegaly.

DESIGN

Case-control study.

METHODS

Bone tissue samples were obtained during transsphenoidal adenomectomy from the sphenoid bone (sella turcica) from 14 patients with clinically and biochemically confirmed acromegaly and 10 patients with clinically non-functioning pituitary adenoma (NFPA) matched by sex and age. Expression of genes involved in the regulation of bone remodeling was studied using quantitative polymerase chain reaction (qPCR).

RESULTS

Of the genes involved in osteoblast and osteoclast activity, only alkaline phosphatase (ALP) mRNA was 50% downregulated in patients with acromegaly. GH excess caused increased expression of the Wnt signaling antagonists (DKK1) and agonists (WNT10B) and changes in the levels of miR involved in mesenchymal stem cell commitment to chondrocytes (miR-199a-5p) or adipocytes (miR-27-5p, miR-125b-5p, miR-34a-5p, miR-188-3p) P < 0.05; q < 0.1. Relevant compensatory mechanisms were found through the changes in miR involved in osteoblastogenesis (miR-210-5p, miR-135a-5p, miR-211, miR-23a-3p, miR-204-5p), but the expression of TWIST1 was 50% downregulated and RUNX2 was unchanged.

CONCLUSIONS

Acromegaly had minimal effects on tested mRNAs specific to osteoblast or osteoclast function except for downregulated ALP expression. The expressions of miR known to be involved in mesenchymal stem cell commitment and downregulated TWIST1 expression suggest acromegaly has a negative effect on osteoblastogenesis.

Evolutionary adaptation revealed by comparative genome analysis of woolly mammoths and elephants

Comparative genomics studies typically limit their focus to single nucleotide variants (SNVs) and that was the case for previous comparisons of woolly mammoth genomes. We extended the analysis to systematically identify not only SNVs but also larger structural variants (SVs) and indels and found multiple mammoth-specific deletions and duplications affecting exons or even complete genes. The most prominent SV found was an amplification of RNase L (with different copy numbers in different mammoth genomes, up to 9-fold), involved in antiviral defense and inflammasome function. This amplification was accompanied by mutations affecting several domains of the protein including the active site and produced different sets of RNase L paralogs in four mammoth genomes likely contributing to adaptations to environmental threats. In addition to immunity and defense, we found many other unique genetic changes in woolly mammoths that suggest adaptations to life in harsh Arctic conditions, including variants involving lipid metabolism, circadian rhythms, and skeletal and body features. Together, these variants paint a complex picture of evolution of the mammoth species and may be relevant in the studies of their population history and extinction.

Dynamics of tRNA fragments and their targets in aging mammalian brain

Background: The progress of next-generation sequencing technologies has unveiled various non-coding RNAs that have previously been considered products of random degradation and attracted only minimal interest. Among small RNA families, microRNA (miRNAs) have traditionally been considered key post-transcriptional regulators. However, recent studies have reported evidence for widespread presence of fragments of tRNA molecules (tRFs) across a range of organisms and tissues, and of tRF involvement in Argonaute complexes.  Methods:To elucidate potential tRF functionality, we compared available RNA sequencing datasets derived from the brains of young, mid-aged and old rats. Using sliding 7-mer windows along a tRF, we searched for putative seed sequences with high numbers of conserved complementary sites within 3' UTRs of 23 vertebrate genomes. We analyzed Gene Ontology term enrichment of predicted tRF targets and compared their transcript levels with targets of miRNAs in the context of age.  Results and Discussion: We detected tRFs originating from 3’- and 5’-ends of tRNAs in rat brains at significant levels. These fragments showed dynamic changes: 3’ tRFs monotonously increased with age, while 5’ tRFs displayed less consistent patterns. Furthermore, 3’ tRFs showed a narrow size range compared to 5’ tRFs, suggesting a difference in their biogenesis mechanisms. Similar to our earlier results in Drosophila and compatible with other experimental findings, we found “seed” sequence locations on both ends of different tRFs. Putative targets of these fragments were found to be enriched in neuronal and developmental functions. Comparison of tRFs and miRNAs increasing in abundance with age revealed small, but distinct changes in brain target transcript levels for these two types of small RNA, with the higher proportion of tRF targets decreasing with age. We also illustrated the utility of tRF analysis for annotating tRNA genes in sequenced genomes.

Age-driven modulation of tRNA-derived fragments in Drosophila and their potential targets

Background

Development of sequencing technologies and supporting computation enable discovery of small RNA molecules that previously escaped detection or were ignored due to low count numbers. While the focus in the analysis of small RNA libraries has been primarily on microRNAs (miRNAs), recent studies have reported findings of fragments of transfer RNAs (tRFs) across a range of organisms.

Results

Here we describe Drosophila melanogaster tRFs, which appear to have a number of structural and functional features similar to those of miRNAs but are less abundant. As is the case with miRNAs, (i) tRFs seem to have distinct isoforms preferentially originating from 5’ or 3’ end of a precursor molecule (in this case, tRNA), (ii) ends of tRFs appear to contain short “seed” sequences matching conserved regions across 12 Drosophila genomes, preferentially in 3’ UTRs but also in introns and exons; (iii) tRFs display specific isoform loading into Ago1 and Ago2 and thus likely function in RISC complexes; (iii) levels of loading in Ago1 and Ago2 differ considerably; and (iv) both tRF expression and loading appear to be age-dependent, indicating potential regulatory changes from young to adult organisms.

Conclusions

We found that Drosophila tRF reads mapped to both nuclear and mitochondrial tRNA genes for all 20 amino acids, while previous studies have usually reported fragments from only a few tRNAs. These tRFs show a number of similarities with miRNAs, including seed sequences. Based on complementarity with conserved Drosophila regions we identified such seed sequences and their possible targets with matches in the 3’UTR regions. Strikingly, the potential target genes of the most abundant tRFs show significant Gene Ontology enrichment in development and neuronal function. The latter suggests that involvement of tRFs in the RNA interfering pathway may play a role in brain activity or brain changes with age.

Reviewers

This article was reviewed by Eugene Koonin, Neil Smalheiser and Alexander Kel.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0081-6) contains supplementary material, which is available to authorized users.

GROM-RD: resolving genomic biases to improve read depth detection of copy number variants

Amplifications or deletions of genome segments, known as copy number variants (CNVs), have been associated with many diseases. Read depth analysis of next-generation sequencing (NGS) is an essential method of detecting CNVs. However, genome read coverage is frequently distorted by various biases of NGS platforms, which reduce predictive capabilities of existing approaches. Additionally, the use of read depth tools has been somewhat hindered by imprecise breakpoint identification. We developed GROM-RD, an algorithm that analyzes multiple biases in read coverage to detect CNVs in NGS data. We found non-uniform variance across distinct GC regions after using existing GC bias correction methods and developed a novel approach to normalize such variance. Although complex and repetitive genome segments complicate CNV detection, GROM-RD adjusts for repeat bias and uses a two-pipeline masking approach to detect CNVs in complex and repetitive segments while improving sensitivity in less complicated regions. To overcome a typical weakness of RD methods, GROM-RD employs a CNV search using size-varying overlapping windows to improve breakpoint resolution. We compared our method to two widely used programs based on read depth methods, CNVnator and RDXplorer, and observed improved CNV detection and breakpoint accuracy for GROM-RD. GROM-RD is available at http://grigoriev.rutgers.edu/software/.

Visualization of nucleotide substitutions in the (micro)transcriptome

BACKGROUND

RNA-related applications of the next-generation sequencing (NGS) technologies require context-specific interpretations: e.g., sequence mismatches may indicate sites of RNA editing, or uneven read coverage often points to mature form of microRNA. Existing visualization tools traditionally show RNA molecules in two dimensions, with their base pairing and the resulting secondary structure. However, it is not straightforward to combine a linear NGS data display with the 2-D RNA depictions.

RESULTS

We present a novel approach for interactive representation of nucleotide substitutions and modifications in the transcribed genome. With the focus on RNA secondary structure in the context of NGS data, it provides intuitive visualization of genomic environment, sequence reads, nucleotide polymorphisms and editing events integrated with the structural and functional elements of both coding and non-coding RNA molecules. Using our approach we present and discuss examples and general trends of polymorphisms and editing in the context of the secondary structure of microRNAs. As expected, most of the substitutions comprised A to G and C to T events, consistent with typical RNA editing patterns. However, we did not observe prevalence of editing in double-stranded regions of the microRNA stem-loop. We describe novel prominent editing event candidates, observed across several small RNA libraries of Drosophila melanogaster.

CONCLUSIONS

In contrast to the existing general tools for NGS data visualization, the power of our approach is not only in the display of read alignments and their counts, but the integration of RNA secondary structure, sequencing depth, and rates/patterns of editing or other modifications. It provides a comprehensive picture, important for large-scale studies and detailed analyses, helping to gain insight into the intricate relationships between different events in RNA biogenesis.

Comparative metagenome analysis of an Alaskan glacier

The temperature in the Arctic region has been increasing in the recent past accompanied by melting of its glaciers. We took a snapshot of the current microbial inhabitation of an Alaskan glacier (which can be considered as one of the simplest possible ecosystems) by using metagenomic sequencing of 16S rRNA recovered from ice/snow samples. Somewhat contrary to our expectations and earlier estimates, a rich and diverse microbial population of more than 2,500 species was revealed including several species of Archaea that has been identified for the first time in the glaciers of the Northern hemisphere. The most prominent bacterial groups found were Proteobacteria, Bacteroidetes, and Firmicutes. Firmicutes were not reported in large numbers in a previously studied Alpine glacier but were dominant in an Antarctic subglacial lake. Representatives of Cyanobacteria, Actinobacteria and Planctomycetes were among the most numerous, likely reflecting the dependence of the ecosystem on the energy obtained through photosynthesis and close links with the microbial community of the soil. Principal component analysis (PCA) of nucleotide word frequency revealed distinct sequence clusters for different taxonomic groups in the Alaskan glacier community and separate clusters for the glacial communities from other regions of the world. Comparative analysis of the community composition and bacterial diversity present in the Byron glacier in Alaska with other environments showed larger overlap with an Arctic soil than with a high Arctic lake, indicating patterns of community exchange and suggesting that these bacteria may play an important role in soil development during glacial retreat.

Impact of age-associated increase in 2'-O-methylation of miRNAs on aging and neurodegeneration in Drosophila

miRNAs exhibit heterogeneity in length and sequence in different biological contexts. Abe et al. found that Drosophila miRNAs undergo isoform pattern changes with age, and an increase of some miRNAs reflects increased 2′-O-methylation of select isoforms. Loading of miRNAs into Ago2, but not Ago1, increased with age. Hen1 and Ago2 mutations caused accelerated neurodegeneration and shorter life span, suggesting that the age-associated increase of 2′-O-methylation of miRNAs affects age-associated processes.

MicroRNAs (miRNAs) are 20- to ∼24-nucleotide (nt) small RNAs that impact a variety of biological processes, from development to age-associated events. To study the role of miRNAs in aging, studies have profiled the levels of miRNAs with time. However, evidence suggests that miRNAs show heterogeneity in length and sequence in different biological contexts. Here, by examining the expression pattern of miRNAs by Northern blot analysis, we found that Drosophila miRNAs show distinct isoform pattern changes with age. Surprisingly, an increase of some miRNAs reflects increased 2′-O-methylation of select isoforms. Small RNA deep sequencing revealed a global increase of miRNAs loaded into Ago2, but not into Ago1, with age. Our data suggest increased loading of miRNAs into Ago2, but not Ago1, with age, indicating a mechanism for differential loading of miRNAs with age between Ago1 and Ago2. Mutations in Hen1 and Ago2, which lack 2′-O-methylation of miRNAs, result in accelerated neurodegeneration and shorter life span, suggesting a potential impact of the age-associated increase of 2′-O-methylation of small RNAs on age-associated processes. Our study highlights that miRNA 2′-O-methylation at the 3′ end is modulated by differential partitioning of miRNAs between Ago1 and Ago2 with age and that this process, along with other functions of Ago2, might impact age-associated events in Drosophila.

A computational study of potential molecular switches that exploit Baird's rule on excited-state aromaticity and antiaromaticity

A series of tentative single-molecule conductance switches which could be triggered by light were examined by computational means using density functional theory (DFT) with non-equilibrium Green's functions (NEGF). The switches exploit the reversal in electron counting rules for aromaticity and antiaromaticity upon excitation from the electronic ground state (S₀) to the lowest ππ* excited singlet and triplet states (S₁ or T₁), as described by Hückel's and Baird's rules, respectively. Four different switches and one antifuse were designed which rely on various photoreactions that either lead from the OFF to the ON states (switches 1, 2 and 4, and antifuse 5) or from the ON to the OFF state (switch 3). The highest and lowest ideal calculated switching ratios are 1175 and 5, respectively, observed for switches 1 and 4. Increased thermal stability of the 1-ON isomer is achieved by benzannulation (switch 1B-OFF/ON). The effects of constrained electrode–electrode distances on activation energies for thermal hydrogen back-transfer from 1-ON to 1-OFF and the relative energies of 1-ON and 1-OFF at constrained geometries were also studied. The switching ratio is strongly distance-dependent as revealed for 1B-ON/OFF where it equals 711 and 148 when the ON and OFF isomers are calculated in electrode gaps with distances confined to either that of the OFF isomer or to that of the ON isomer, respectively.

Identification of vibrational signatures from short chains of interlinked molecule-nanoparticle junctions obtained by inelastic electron tunnelling spectroscopy

Short chains containing a series of metal-molecule-nanoparticle nanojunctions are a nano-material system with the potential to give electrical signatures close to those from single molecule experiments while enabling us to build portable devices on a chip. Inelastic electron tunnelling spectroscopy (IETS) measurements provide one of the most characteristic electrical signals of single and few molecules. In interlinked molecule-nanoparticle (NP) chains containing typically 5-7 molecules in a chain, the spectrum is expected to be a superposition of the vibrational signatures of individual molecules. We have established a stable and reproducible molecule-AuNP multi-junction by placing a few 1,8-octanedithiol (ODT) molecules onto a versatile and portable nanoparticle-nanoelectrode platform and measured for the first time vibrational molecular signatures at complex and coupled few-molecule-NP junctions. From quantum transport calculations, we model the IETS spectra and identify vibrational modes as well as the number of molecules contributing to the electron transport in the measured spectra.

Assessment of a nanoparticle bridge platform for molecular electronics measurements

A combination of electron beam lithography, photolithography and focused ion beam milling was used to create a nanogap platform, which was bridged by gold nanoparticles in order to make electrical measurements and assess the platform under ambient conditions. Non-functionalized electrodes were tested to determine the intrinsic response of the platform and it was found that creating devices in ambient conditions requires careful cleaning and awareness of the contributions contaminants may make to measurements. The platform was then used to make measurements on octanethiol (OT) and biphenyldithiol (BPDT) molecules by functionalizing the nanoelectrodes with the molecules prior to bridging the nanogap with nanoparticles. Measurements on OT show that it is possible to make measurements on relatively small numbers of molecules, but that a large variation in response can be expected when one of the metal-molecule junctions is physisorbed, which was partially explained by attachment of OT molecules to different sites on the surface of the Au electrode using a density functional theory calculation. On the other hand, when dealing with BPDT, high yields for device creation are very difficult to achieve under ambient conditions. Significant hysteresis in the I-V curves of BPDT was also observed, which was attributed primarily to voltage induced changes at the interface between the molecule and the metal.

Electron transport in stretched monoatomic gold wires

The conductance of monoatomic gold wires containing 3-7 gold atoms has been obtained from ab initio calculations. The transmission is found to vary significantly depending on the wire stretching and the number of incorporated atoms. Such oscillations are determined by the electronic structure of the one-dimensional (1D) part of the wire between the contacts. Our results indicate that the conductivity of 1D wires can be suppressed without breaking the contact.

Solar Wind-Induced Atmospheric Erosion at Mars: First Results from ASPERA-3 on Mars Express

The Analyzer of Space Plasma and Energetic Atoms (ASPERA) on board the Mars Express spacecraft found that solar wind plasma and accelerated ionospheric ions may be observed all the way down to the Mars Express pericenter of 270 kilometers above the dayside planetary surface. This is very deep in the ionosphere, implying direct exposure of the martian topside atmosphere to solar wind plasma forcing. The low-altitude penetration of solar wind plasma and the energization of ionospheric plasma may be due to solar wind irregularities or perturbations, to magnetic anomalies at Mars, or both.

Hormonal changes during long-term isolation

Confinement and inactivity induce considerable psychological and physiological modifications through social and sensory deprivation. The aim of the SFINCSS-99 experiment was to determine the cardiovascular and hormonal pattern of blood volume regulation during long-term isolation and confinement. Simulation experiments were performed in pressurized chambers similar in size to the volumes of modern space vehicles. Group I consisted of four Russian male volunteers, who spent 240 days in a 100-m(3 )chamber. Group II included four males (one German and three Russians) who spent 110 days in isolation (200-m(3) module). The blood samples, taken before, during and after the isolation period, were used to determine haematocrit (Ht), growth hormone (GH), active renin, aldosterone, and osmolality levels. From the urine samples, electrolytes, osmolality, nitrites, nitrates, cortisol, antidiuretic hormone (ADH), aldosterone, normetanephrine and metanephrine levels were determined. The increase in plasma volume (PV) that is associated with a tendency for a decrease in plasma active renin is likely to be due to decreased sympathetic activity, and concords with the changes in urinary catecholamine levels during confinement. Urinary catecholamine levels were significantly higher during the recovery period than during confinement. This suggests that the sympathoadrenal system was activated, and concords with the increase in heart rate. Vascular resistance is determined by not only the vasoconstrictor but also vasodilator systems. The ratio of nitrite/nitrate in urine, as an indicator of nitric oxide release, did not reveal any significant changes. Analysis of data suggests that the duration of the isolation was a main factor involved in the regulation of hormones.

Deuterium-free water (1H2O) in complex life-support systems of long-term space missions

Heavy water containing deuterium displays toxic property. It is stated that any quantity of a heavy isotope of hydrogen--deuterium--is undesirable to animals and plants. It was earlier shown by us that physical-chemical life support systems on board the "MIR" station fractionate (change) isotopes of hydrogen, oxygen and carbon. Therefore, the problem of regenerative systems in habitable space objects should include removal, from water, of a heavy stable isotope of hydrogen--deuterium. In this article we consider one method of obtaining deuterium-free water--decomposition of distillate water in an electrolyser to hydrogen and oxygen with subsequent synthesis in a catalytic or high-temperature reactor. The influence of deuterium-free water on the growth and development of Arabidopsis thaliana and Japanese quail is investigated. It is shown that with the help of the electrolysis method it is possible to fabricate water containing 80% less deuterium in comparison with SMOW. Experimentally, it is proved on a culture of Arabidopsis thaliana and Japanese quail that water with reduced contents of deuterium (80%) displays positive biological activity.

Influence of hypergravity on hypothalamic vasopressin and oxytocin neurons in rats

This study was aimed to evaluate the reaction of the vasopressin (VP) and oxytocin (OT) neurons of the supraoptic nucleus (SON) in rats to single or repeated hypergravity (HG). Special attention was paid to the tyrosine hydroxylase (TH) expression in VP neurons as a marker of the neuron activation. Rats were revolved in a centrifuge with overloading 2G for 5 days or 34 days as well as for 34 days plus 5 days with an interval of 39 days between two rotations. Control rats were kept in a centrifuge room. Radioimmununoassay, quantitative and semi-quantitative immunocytochemistry and in situ hybridization were used to evaluate: a) VP concentration in the pituitary posterior lobe (PL) and in plasma; b) the number of VP-, OT- and TH-immunoreactive neurons in the SON; c) the optic density of VP-, OT- and TH-immunoreactive materials in cell bodies (SON) and distal axons (PL), d) the optic density of VP and OT mRNAs signals (S35) in the whole SON on microfilms. According to our data, VP neurons were strongly activated during HG (5 days or 34 days) that was manifested in the functional hypertrophy of the neurons, greatly increased concentrations of VP mRNA in the SON and VP in plasma, the onset of the TH expression. The neurons showed initially (5 days) the functional insufficiency (VP release > VP synthesis) followed by their adaptation (subsequent 29 days) to the increased need in VP (VP release < VP synthesis). No reaction of VP neurons was observed to repeated HG. In contrast to VP neurons, OT neurons did not react to short-term HG or showed functional depression after the long-term treatment.

A relationship between gene expression and protein interactions on the proteome scale: analysis of the bacteriophage T7 and the yeast Saccharomyces cerevisiae

The relationship between the similarity of expression patterns for a pair of genes and interaction of the proteins they encode is demonstrated both for the simple genome of the bacteriophage T7 and the considerably more complex genome of the yeast Saccharomyces cerevisiae. Statistical analysis of large-scale gene expression and protein interaction data shows that protein pairs encoded by co-expressed genes interact with each other more frequently than with random proteins. Furthermore, the mean similarity of expression profiles is significantly higher for respective interacting protein pairs than for random ones. Such coupled analysis of gene expression and protein interaction data may allow evaluation of the results of large-scale gene expression and protein interaction screens as demonstrated for several publicly available datasets. The role of this link between expression and interaction in the evolution from monomeric to oligomeric protein structures is also discussed.

Effects of exposure to space flight on endocrine regulations in experimental animals

This minireview summarizes the results of the observations on changes in endocrine functions of rats exposed to space flights for various periods. The results found after space flights are compared with those obtained from rats in acute or repeated restrain stress. A slight increase of plasma catecholamine levels was observed in rats after space flight of longer duration (>14 days), but no changes in catecholamine content in the activity of catecholamine synthesizing enzymes were noted in adrenal medulla and in hypothalamus. The norepinephrine content was, however, decreased in several nuclei selected from hypothalamus of flight rats. Plasma corticosterone levels were increased after space flight and morphological examination of pituitary showed elevated activity of corticotrophs. However, the plasma levels of ACTH were not increased in rats 6 hours after space flight. These changes in plasma hormone levels affected the activity of enzymes involved in metabolism of amino acids in liver and lipolysis in adipose tissue. The plasma levels of testosterone and triiodothyronine were diminished after space flight suggesting the suppression of the thyroid and gonadal activity. Increase of plasma insulin and glucose levels were found in rats after space flight, but the glucagon values were not changed. Comparing these results from flight rats with the animals exposed to acute or repeated stress indicate that long stay in microgravity do not represent very intensive stressogenic stimulus for adrenocortical and sympatho-adrenomedullar systems, and hormone alterations observed after space flight may be due to acute gravitational stress resulting from a return to Earth gravity. Therefore further studies including the inflight animal experiments on a board of International Space Station are necessary for elucidation of the effects of microgravity on endocrine functions.

Mutations induced by bacteriophage T7 RNA polymerase and their effects on the composition of the T7 genome

We show here that transcription by the bacteriophage T7 RNA polymerase increases the deamination of cytosine bases in the non-transcribed strand to uracil, causing C to T mutations in that strand. Under optimal conditions, the mutation frequency increases about fivefold over background, and is similar to that seen with the Escherichia coli RNA polymerase. Further, we found that a mutant T7 RNA polymerase with a slower rate of elongation caused more cytosine deaminations than its wild-type parent. These results suggest that promoting cytosine deamination in the non-transcribed strand is a general property of transcription in E. coli and is dependent on the length of time the transcription bubble stays open during elongation. To see if transcription-induced mutations have influenced the evolution of bacteriophage T7, we analyzed its genome for a bias in base composition. Our analysis showed a significant excess of thymine over cytosine bases in the highly transcribed regions of the genome. Moreover, the average value of this bias correlated well with the levels of transcription of different genomic regions. Our results indicate that transcription-induced mutations have altered the composition of bacteriophage T7 genome and suggest that this may be a significant force in genome evolution.

Strand-specific compositional asymmetries in double-stranded DNA viruses

Analysis of 22 complete sequences of double-stranded DNA viruses reveals striking compositional asymmetries between leading and lagging, and between transcribed and non-transcribed strands. In all bi-directionally replicated genomes analyzed, the observed leading strand GC skew (measuring relative excess of guanines versus cytosines) is different from that in the lagging strand. In most of these genomes GC skew switches polarity close to replication origins. GC skew changes linearly across adenovirus linear genomes, which replicate from one end. In papillomavirus, GC skew is positive in one half of the genome where transcription and replication proceed in the same direction, and is close to zero in the other half with divergent transcription and replication. Possible contributions of these two processes (and associated repair mechanisms) as well as other potential sources of strand bias in the observed asymmetries are discussed. Use of cumulative skew plots for genome comparisons is demonstrated on the example of herpes simplex virus.

Reusable graphical interface to genome information resources

This paper describes a prototype genome display and query system for the World Wide Web, which could play the role of a graphical interactive gateway to online genome information services. It provides a uniform interface to display mapping and sequencing data for the human, mouse and yeast genomes and could be easily extended to accommodate more information as it becomes available. This system uses a Java applet, DerBrowser, for delivering interactive content to an end user. The architecture and functionality of this applet are described, with respect to views of both users and data providers.

Influence of bisphosphonates on calcium metabolism and bone tissue during simulation of the physiological effects of microgravity

At present a significant experience has been gathered in experimental and clinical use of synthetic derivatives of bisphosphonic acid. Purpose of our study was to investigate the effects of bisphosphonates on calcium metabolism and bone tissue in humans and animals under the conditions simulating the physiological effects of microgravity. In a series of experiments with animals the effectiveness of various preparations and treatment plans was compared Action of 1-hydroxyethyliden-1, 1-bisphosphonic acid (EHBP, xydiphon) on calcium metabolism and bone tissue in human was the subject of experiments with long-term head-down tilt (HDT).

A distributed environment for physical map construction

MOTIVATION

With the main focus of the Human Genome Project shifting to sequencing, bioinformatics support for constructing large-scale genomic maps of other organisms is still required. We attempt to provide for this with our work, aimed at the delivery of robust and user-friendly contig-building software on the WWW.

RESULTS

We present a prototype distributed analytical environment for molecular biologists working in the area of genomic mapping. It consists of the WWW server for constructing contigs from users' data with a hypertext output connected to Java-based map visualization software.

AVAILABILITY

Freely available on http://www.mpimg-berlin-dahlem.mpg. de/ approximately andy/server/

CONTACT

andy@rag3.rz-berlin.mpg.de

Analyzing genomes with cumulative skew diagrams

A novel method of cumulative diagrams shows that the nucleotide composition of a microbial chromosome changes at two points separated by about a half of its length. These points coincide with sites of replication origin and terminus for all bacteria where such sites are known. The leading strand is found to contain more guanine than cytosine residues. This fact is used to predict origin and terminus locations in other bacterial and archaeal genomes. Local changes, visible as diagram distortions, may represent recent genome rearrangements, as demonstrated for two strains of Escherichia coli . Analysis of the diagrams of viral and mitochondrial genomes suggests a link between the base composition bias and the time spent by DNA in a single stranded state during replication.

IXDB, an X chromosome integrated database

The integrated X chromosome database (IXDB) is a repository for physical mapping data of the human X chromosome. Its current content is the result of a strict integration of data stemming from many different sources. The main features of IXDB include a flexible and extendible schema, a comfortable and fully cross-referenced WWW interface (http://ixdb.mpimg-berlin-dahlem.mpg.de ) and a graphical map viewer implemented in JAVA. The database stores objects used in physical mapping as well as the maps resulting from this work, but a strong emphasis is placed on recording experiments that connect objects together. This should greatly contribute to fulfilling one of the major goals of the database: to support the construction of an integrated physical, genetic, transcript and sequence map of the human X chromosome.

Effects of three days of dry immersion on heart rate and blood pressure variabilities during head-up tilting in humans

In this study we examined the effects of 3 days of dry immersion on heart rate variability (HRV) and blood pressure variability (BPV) in response to head-up tilting in 4 healthy young subjects aged 21 to 36 years. Resting value of the high-frequency (HF) power of BPV and HRV decreased, while the low-frequency (LF)/HF ratio of HRV and blood pressure increased after the dry immersion. The HF power of HRV decreased, while the LF/HF ratio of HRV and heart rate increased during head-up tilting before and after dry immersion. All powers of BPV and BP were unchanged during head-up tilting. The values of the decrease in the HF power of HRV and the increase in the LF/HF ratio during head-up tilting after dry immersion were larger than those before. These results suggest that 3 days of dry immersion altered the autonomic balance toward the sympathetic dominant at rest, and that cardiac function may play a crucial role in BP maintenance during head-up tilting compared with vasoconstrictor function.