Reprever: resolving low-copy duplicated sequences using template driven assembly

Genomic sequence duplication is an important mechanism for genome evolution, often resulting in large sequence variations with implications for disease progression. Although paired-end sequencing technologies are commonly used for structural variation discovery, the discovery of novel duplicated sequences remains an unmet challenge. We analyze duplicons starting from identified high-copy number variants. Given paired-end mapped reads, and a candidate high-copy region, our tool, Reprever, identifies (a) the insertion breakpoints where the extra duplicons inserted into the donor genome and (b) the actual sequence of the duplicon. Reprever resolves ambiguous mapping signatures from existing homologs, repetitive elements and sequencing errors to identify breakpoint. At each breakpoint, Reprever reconstructs the inserted sequence using profile hidden Markov model (PHMM)-based guided assembly. In a test on 1000 artificial genomes with simulated duplication, Reprever could identify novel duplicates up to 97% of genomes within 3 bp positional and 1% sequence errors. Validation on 680 fosmid sequences identified and reconstructed eight duplicated sequences with high accuracy. We applied Reprever to reanalyzing a re-sequenced data set from the African individual NA18507 to identify >800 novel duplicates, including insertions in genes and insertions with additional variation. polymerase chain reaction followed by capillary sequencing validated both the insertion locations of the strongest predictions and their predicted sequence.     

Seasonal acclimatization to the hot summer over 60 days in the Republic of Korea suppresses sweating sensitivity during passive heating

The main objective of this study was to determine the central mechanisms involved in suppression of thermal sweating after seasonal acclimatization (SA) during passive heating (immersing the legs in 43 °C hot water for 30 min). Testing was performed in July (before-SA) and August (after-SA) [25.2±2.2 °C, 73.9±10.3% relative humidity (RH), Cheonan (Chungnam,126° 52′N, 33.38′E), in the Republic of Korea. All experiments were carried out in an automated climatic chamber (25.0±0.5 °C and RH 60.0±3.00%). Twelve healthy men (height, 174.6±5.40 cm; weight, 65.4±5.71 kg; age, 22.7±2.90 yr) participated. The local sweat onset time was delayed in the after-SA compared to that in the before-SA (p<0.001). The local sweat rate and whole body sweat loss volume decreased in the after-SA compared to those in the before-SA (p<0.001). In addition, evaporative loss volume decreased significantly in the after-SA compared to that in the before-SA [chest, upper-back, thigh and forearm (p<0.001)]. Changes in tympanic temperature and mean body temperature were significantly lower (p<0.05) and the basal metabolic rate decreased significantly in the after-SA compared to those in the before-SA (p<0.001). These results suggest that maintenance of a lower body temperature and basal metabolic rate can occur and blunt the central sudomotor mechanisms following seasonal acclimatization, which suppresses sweating sensitivity.     

Algivirga pacifica gen. nov., sp. nov., a Novel Agar-Degrading Marine Bacterium of the Family Flammeovirgaceae Isolated from Micronesia

An aerobic, Gram-negative, coccoid to short rod-shaped and non-flagellated marine bacterial strain S354^T was isolated from seawater of Micronesia. The strain was capable to degrade agar-forming slight depression into agar plate. Growth occurred at a temperature range of 12–44 °C, a pH range of 5–9, and a salinity range of 1–7 % (w/v) NaCl. Phylogenetic analyses based on 16S rRNA gene sequences suggested that S354^T belongs to the family Flammeovirgaceae. The novel strain was most closely related to Limibacter armeniacum YM 11-185^T with similarity of 92.5 %. The DNA G+C content was 43.8 mol%. The major fatty acids (>10 %) were iso-C_15:0 and C_16:1 ω5c. The predominant isoprenoid quinone was determined to be MK-7. Polar lipid profile of S354^T consisted of phosphatidylethanolamine, unknown polar lipid, and unknown glycolipids. Based on the phenotypic, phylogenetic, biochemical, and physiological tests conducted in this study, S354^T is proposed to represent a type strain of a novel genus and species. The 16S rRNA gene sequence of S354^T is registered in GenBank under the accession number JQ639084. The type of strain Algivirga pacifica gen. nov., sp. nov. is S354^T (=KCCM 90107^T=JCM 18326^T).     

MuPIT interactive: webserver for mapping variant positions to annotated, interactive 3D structures

Mutation position imaging toolbox (MuPIT) interactive is a browser-based application for single-nucleotide variants (SNVs), which automatically maps the genomic coordinates of SNVs onto the coordinates of available three-dimensional (3D) protein structures. The application is designed for interactive browser-based visualization of the putative functional relevance of SNVs by biologists who are not necessarily experts either in bioinformatics or protein structure. Users may submit batches of several thousand SNVs and review all protein structures that cover the SNVs, including available functional annotations such as binding sites, mutagenesis experiments, and common polymorphisms. Multiple SNVs may be mapped onto each structure, enabling 3D visualization of SNV clusters and their relationship to functionally annotated positions. We illustrate the utility of MuPIT interactive in rationalizing the impact of selected polymorphisms in the PharmGKB database, somatic mutations identified in the Cancer Genome Atlas study of invasive breast carcinomas, and rare variants identified in the exome sequencing project. MuPIT interactive is freely available for non-profit use at http://mupit.icm.jhu.edu.     

Effects of Ethanol Exposure During Early Pregnancy in Hyperactive, Inattentive and Impulsive Behaviors and MeCP2 Expression in Rodent Offspring

Prenatal exposure to alcohol has consistently been associated with adverse effects on neurodevelopment, which is collectively called fetal alcohol spectrum disorder (FASD). Increasing evidence suggest that prenatal exposure to alcohol increases the risk of developing attention deficit/hyperactivity disorder-like behavior in human. In this study, we investigated the behavioral effects of prenatal exposure to EtOH in offspring mice and rats focusing on hyperactivity and impulsivity. We also examined changes in dopamine transporter and MeCP2 expression, which may underlie as a key neurobiological and epigenetic determinant in FASD and hyperactive, inattentive and impulsive behaviors. Mouse or rat offspring born from dam exposed to alcohol during pregnancy (EtOH group) showed hyper locomotive activity, attention deficit and impulsivity. EtOH group also showed increased dopamine transporter and norepinephrine transporter level compared to control group in the prefrontal cortex and striatum. Prenatal exposure to EtOH also significantly decreased the expression of MeCP2 in both prefrontal cortex and striatum. These results suggest that prenatal exposure to EtOH induces hyperactive, inattentive and impulsive behaviors in rodent offspring that might be related to global epigenetic changes as well as aberration in catecholamine neurotransmitter transporter system.     

Effect of GABA Receptor Agonists or Antagonists Injected Spinally on the Blood Glucose Level in Mice

The possible roles of gamma-amino butyric acid (GABA) receptors located in the spinal cord for the regulation of the blood glucose level were studied in ICR mice. We found in the present study that intrathecal (i.t.) injection with baclofen (a GABA_B receptor agonist; 1–10 μg/5 μl) or bicuculline (a GABA_A receptor antagonist; 1–10 μg/5 μl) caused an elevation of the blood glucose level in a dose-dependent manner. The hyperglycemic effect induced by baclofen was more pronounced than that induced by bicuculline. However, muscimol (a GABA_A receptor agonist; 1–5 μg/5 μl) or phaclofen (a GABA_B receptor antagonist; 5–10 μg/5 μl) administered i.t. did not affect the blood glucose level. Baclofen–induced elevation of the blood glucose was dose-dependently attenuated by phaclofen. Furthermore, i.t. pretreatment with pertussis toxin (PTX; 0.05 or 0.1 μg/5 μl) for 6 days dose-dependently reduced the hyperglycemic effect induced by baclofen. Our results suggest that GABA_B receptors located in the spinal cord play important roles for the elevation of the blood glucose level. Spinally located PTX-sensitive G-proteins appear to be involved in hyperglycemic effect induced by baclofen. Furthermore, inactivation of GABA_A receptors located in the spinal cord appears to be responsible for tonic up-regulation of the blood glucose level.     

Time-course Changes in Immunoreactivities of Glucokinase and Glucokinase Regulatory Protein in the Gerbil Hippocampus Following Transient Cerebral Ischemia

Glucose is a main energy source for normal brain functions. Glucokinase (GK) plays an important role in glucose metabolism as a glucose sensor, and GK activity is modulated by glucokinase regulatory protein (GKRP). In this study, we examined the changes of GK and GKRP immunoreactivities in the gerbil hippocampus after 5 min of transient global cerebral ischemia. In the sham-operated-group, GK and GKRP immunoreactivities were easily detected in the pyramidal neurons of the stratum pyramidale of the hippocampus. GK and GKRP immunoreactivities in the pyramidal neurons were distinctively decreased in the hippocampal CA1 region (CA), not CA2/3, 3 days after ischemia–reperfusion (I–R). Five days after I–R, GK and GKRP immunoreactivities were hardly detected in the CA1, not CA2/3, pyramidal neurons; however, at this point in time, GK and GKRP immunoreactivities were newly expressed in astrocytes, not microglia, in the ischemic CA1. In brief, GK and GKRP immunoreactivities are changed in pyramidal neurons and newly expressed in astrocytes in the ischemic CA1 after transient cerebral ischemia. These indicate that changes of GK and GKRP expression may be related to the ischemia-induced neuronal damage/death.     

The novel MER transposon-derived miRNAs in human genome

MicroRNAs (miRNAs) are small RNA molecules (~ 20–30 nucleotides) that generally act in gene silencing and translational repression through the RNA interference pathway. They generally originate from intergenic genomic regions, but some are found in genomic regions that have been characterized such as introns, exons, and transposable elements (TE). To identify the miRNAs that are derived from palindromic MERs, we analyzed MER paralogs in human genome. The structures of the palindromic MERs were similar to the hairpin structure of miRNA in humans. Three miRNAs derived from MER96 located on chromosome 3, and MER91C paralogs located on chromosome 8 and chromosome 17 were identified in HeLa, HCT116, and HEK293 cell lines. The interactions between these MER-derived miRNAs and AGO1, AGO2, and AGO3 proteins were validated by immunoprecipitation assays. The data suggest that miRNAs derived from transposable elements could widely affect various target genes in the human genome.     

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae)

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.