Dusty Planetary Systems

Extensive photometric stellar surveys show that many main sequence stars show emission at infrared and longer wavelengths that is in excess of the stellar photosphere; this emission is thought to arise from circumstellar dust. The presence of dust disks is confirmed by spatially resolved imaging at infrared to millimeter wavelengths (tracing the dust thermal emission), and at optical to near infrared wavelengths (tracing the dust scattered light). Because the expected lifetime of these dust particles is much shorter than the age of the stars (>10 Myr), it is inferred that this solid material not primordial, i.e. the remaining from the placental cloud of gas and dust where the star was born, but instead is replenished by dust-producing planetesimals. These planetesimals are analogous to the asteroids, comets and Kuiper Belt objects (KBOs) in our Solar system that produce the interplanetary dust that gives rise to the zodiacal light (tracing the inner component of the Solar system debris disk). The presence of these "debris disks" around stars with a wide range of masses, luminosities, and metallicities, with and without binary companions, is evidence that planetesimal formation is a robust process that can take place under a wide range of conditions. This chapter is divided in two parts. Part I discusses how the study of the Solar system debris disk and the study of debris disks around other stars can help us learn about the formation, evolution and diversity of planetary systems by shedding light on the frequency and timing of planetesimal formation, the location and physical properties of the planetesimals, the presence of long-period planets, and the dynamical and collisional evolution of the system. Part II reviews the physical processes that affect dust particles in the gas-free environment of a debris disk and their effect on the dust particle size and spatial distribution.Comment: 68 pages, 25 figures. To be published in "Solar and Planetary Systems" (P. Kalas and L. French, Eds.), Volume 3 of the series "Planets, Stars and Stellar Systems" (T.D. Oswalt, Editor-in-chief), Springer 201

svclassify: a method to establish benchmark structural variant calls

The human genome contains variants ranging in size from small single nucleotide polymorphisms (SNPs) to large structural variants (SVs). High-quality benchmark small variant calls for the pilot National Institute of Standards and Technology (NIST) Reference Material (NA12878) have been developed by the Genome in a Bottle Consortium, but no similar high-quality benchmark SV calls exist for this genome. Since SV callers output highly discordant results, we developed methods to combine multiple forms of evidence from multiple sequencing technologies to classify candidate SVs into likely true or false positives. Our method (svclassify) calculates annotations from one or more aligned bam files from many high-throughput sequencing technologies, and then builds a one-class model using these annotations to classify candidate SVs as likely true or false positives. We first used pedigree analysis to develop a set of high-confidence breakpoint-resolved large deletions. We then used svclassify to cluster and classify these deletions as well as a set of high-confidence deletions from the 1000 Genomes Project and a set of breakpoint-resolved complex insertions from Spiral Genetics. We find that likely SVs cluster separately from likely non-SVs based on our annotations, and that the SVs cluster into different types of deletions. We then developed a supervised one-class classification method that uses a training set of random non-SV regions to determine whether candidate SVs have abnormal annotations different from most of the genome. To test this classification method, we use our pedigree-based breakpoint-resolved SVs, SVs validated by the 1000 Genomes Project, and assembly-based breakpoint-resolved insertions, along with semi-automated visualization using svviz. We find that candidate SVs with high scores from multiple technologies have high concordance with PCR validation and an orthogonal consensus method MetaSV (99.7 % concordant), and candidate SVs with low scores are questionable. We distribute a set of 2676 high-confidence deletions and 68 high-confidence insertions with high svclassify scores from these call sets for benchmarking SV callers. We expect these methods to be particularly useful for establishing high-confidence SV calls for benchmark samples that have been characterized by multiple technologies.https://doi.org/10.1186/s12864-016-2366-

Accurate detection of complex structural variations using single-molecule sequencing

Structural variations are the greatest source of genetic variation, but they remain poorly understood because of technological limitations. Single-molecule long-read sequencing has the potential to dramatically advance the field, although high error rates are a challenge with existing methods. Addressing this need, we introduce open-source methods for long-read alignment (NGMLR; https://github.com/philres/ngmlr ) and structural variant identification (Sniffles; https://github.com/fritzsedlazeck/Sniffles ) that provide unprecedented sensitivity and precision for variant detection, even in repeat-rich regions and for complex nested events that can have substantial effects on human health. In several long-read datasets, including healthy and cancerous human genomes, we discovered thousands of novel variants and categorized systematic errors in short-read approaches. NGMLR and Sniffles can automatically filter false events and operate on low-coverage data, thereby reducing the high costs that have hindered the application of long reads in clinical and research settings

Hippocampal CA3 Transcriptome Signature Correlates with Initial Precipitating Injury in Refractory Mesial Temporal Lobe Epilepsy

Background: Prolonged febrile seizures constitute an initial precipitating injury (IPI) commonly associated with refractory mesial temporal lobe epilepsy (RMTLE). in order to investigate IPI influence on the transcriptional phenotype underlying RMTLE we comparatively analyzed the transcriptomic signatures of CA3 explants surgically obtained from RMTLE patients with (FS) or without (NFS) febrile seizure history. Texture analyses on MRI images of dentate gyrus were conducted in a subset of surgically removed sclerotic hippocampi for identifying IPI-associated histo-radiological alterations.Methodology/Principal Findings: DNA microarray analysis revealed that CA3 global gene expression differed significantly between FS and NFS subgroups. An integrative functional genomics methodology was used for characterizing the relations between GO biological processes themes and constructing transcriptional interaction networks defining the FS and NFS transcriptomic signatures and its major gene-gene links (hubs). Co-expression network analysis showed that: i) CA3 transcriptomic profiles differ according to the IPI; ii) FS distinctive hubs are mostly linked to glutamatergic signalization while NFS hubs predominantly involve GABAergic pathways and neurotransmission modulation. Both networks have relevant hubs related to nervous system development, what is consistent with cell genesis activity in the hippocampus of RMTLE patients. Moreover, two candidate genes for therapeutic targeting came out from this analysis: SSTR1, a relevant common hub in febrile and afebrile transcriptomes, and CHRM3, due to its putative role in epilepsy susceptibility development. MRI texture analysis allowed an overall accuracy of 90% for pixels correctly classified as belonging to FS or NFS groups. Histological examination revealed that granule cell loss was significantly higher in FS hippocampi.Conclusions/Significance: CA3 transcriptional signatures and dentate gyrus morphology fairly correlate with IPI in RMTLE, indicating that FS-RMTLE represents a distinct phenotype. These findings may shed light on the molecular mechanisms underlying refractory epilepsy phenotypes and contribute to the discovery of novel specific drug targets for therapeutic interventions

Nanopore sequencing and assembly of a human genome with ultra-long reads

We report the sequencing and assembly of a reference genome for the human GM12878 Utah/Ceph cell line using the MinION (Oxford Nanopore Technologies) nanopore sequencer. 91.2 Gb of sequence data, representing ~30× theoretical coverage, were produced. Reference-based alignment enabled detection of large structural variants and epigenetic modifications. De novo assembly of nanopore reads alone yielded a contiguous assembly (NG50 ~3 Mb). Next, we developed a protocol to generate ultra-long reads (N50 > 100kb, up to 882 kb). Incorporating an additional 5×-coverage of these data more than doubled the assembly contiguity (NG50 ~6.4 Mb). The final assembled genome was 2,867 million bases in size, covering 85.8% of the reference. Assembly accuracy, after incorporating complementary short-read sequencing data, exceeded 99.8%. Ultra-long reads enabled assembly and phasing of the 4 Mb major histocompatibility complex (MHC) locus in its entirety, measurement of telomere repeat length and closure of gaps in the reference human genome assembly GRCh38

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios.Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature.Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays).Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties