Improving Power Spectral Estimation using Multitapering: Precise asteroseismic modeling of stars, exoplanets, and beyond

Asteroseismic time-series data have imprints of stellar oscillation modes, whose detection and characterization through time-series analysis allows us to probe stellar interiors physics. Such analyses usually occur in the Fourier domain by computing the Lomb-Scargle (LS) periodogram, an estimator of the \textit{power spectrum} underlying unevenly-sampled time-series data. However, the LS periodogram suffers from the statistical problems of (1) inconsistency (or noise) and (2) bias due to high spectral leakage. In addition, it is designed to detect strictly periodic signals but is unsuitable for non-sinusoidal periodic or quasi-periodic signals. Here, we develop a multitaper spectral estimation method that tackles the inconsistency and bias problems of the LS periodogram. We combine this multitaper method with the Non-Uniform Fast Fourier Transform (\texttt{mtNUFFT}) to more precisely estimate the frequencies of asteroseismic signals that are non-sinusoidal periodic (e.g., exoplanet transits) or quasi-periodic (e.g., pressure modes). We illustrate this using a simulated and the Kepler-91 red giant light curve. Particularly, we detect the Kepler-91b exoplanet and precisely estimate its period, $6.246 \pm 0.002$ days, in the frequency domain using the multitaper F-test alone. We also integrate \texttt{mtNUFFT} into the \texttt{PBjam} package to obtain a Kepler-91 age estimate of $3.96 \pm 0.48$ Gyr. This $36$\% improvement in age precision relative to the $4.27 \pm 0.75$ Gyr APOKASC-2 (uncorrected) estimate illustrates that \texttt{mtNUFFT} has promising implications for Galactic archaeology, in addition to stellar interiors and exoplanet studies. Our frequency analysis method generally applies to time-domain astronomy and is implemented in the public Python package \texttt{tapify}, available at \url{https://github.com/aaryapatil/tapify}.Comment: 32 pages (3 pages in the Appendix), 14 figures, 2 tables, Submitted to A

GLUE: a flexible software system for virus sequence data

Background: Virus genome sequences, generated in ever-higher volumes, can provide new scientific insights and inform our responses to epidemics and outbreaks. To facilitate interpretation, such data must be organised and processed within scalable computing resources that encapsulate virology expertise. GLUE (Genes Linked by Underlying Evolution) is a data-centric bioinformatics environment for building such resources. The GLUE core data schema organises sequence data along evolutionary lines, capturing not only nucleotide data but associated items such as alignments, genotype definitions, genome annotations and motifs. Its flexible design emphasises applicability to different viruses and to diverse needs within research, clinical or public health contexts. Results: HCV-GLUE is a case study GLUE resource for hepatitis C virus (HCV). It includes an interactive public web application providing sequence analysis in the form of a maximum-likelihood-based genotyping method, antiviral resistance detection and graphical sequence visualisation. HCV sequence data from GenBank is categorised and stored in a large-scale sequence alignment which is accessible via web-based queries. Whereas this web resource provides a range of basic functionality, the underlying GLUE project can also be downloaded and extended by bioinformaticians addressing more advanced questions. Conclusion: GLUE can be used to rapidly develop virus sequence data resources with public health, research and clinical applications. This streamlined approach, with its focus on reuse, will help realise the full value of virus sequence data

Hepatitis C and the absence of genomic data in low-income countries: a barrier on the road to elimination?

Following the development of highly effective direct acting antiviral (DAA) compounds for the treatment of the hepatitis C virus (HCV), WHO has set out plans for disease eradication by 2030. Many barriers must be surmounted before this can be achieved, including buy-in from governments and policy makers, reduced drug costs, and improved infrastructure for the pathway from diagnosis to treatment. A comprehensive set of guidelines was produced by WHO in 2014, updated in 2016, and they are due to be revised later this year

Rice Galaxy: An open resource for plant science

Background: Rice molecular genetics, breeding, genetic diversity, and allied research (such as rice-pathogen interaction) have adopted sequencing technologies and high-density genotyping platforms for genome variation analysis and gene discovery. Germplasm collections representing rice diversity, improved varieties, and elite breeding materials are accessible through rice gene banks for use in research and breeding, with many having genome sequences and high-density genotype data available. Combining phenotypic and genotypic information on these accessions enables genome-wide association analysis, which is driving quantitative trait loci discovery and molecular marker development. Comparative sequence analyses across quantitative trait loci regions facilitate the discovery of novel alleles. Analyses involving DNA sequences and large genotyping matrices for thousands of samples, however, pose a challenge to non−computer savvy rice researchers. Findings: The Rice Galaxy resource has shared datasets that include high-density genotypes from the 3,000 Rice Genomes project and sequences with corresponding annotations from 9 published rice genomes. The Rice Galaxy web server and deployment installer includes tools for designing single-nucleotide polymorphism assays, analyzing genome-wide association studies, population diversity, rice−bacterial pathogen diagnostics, and a suite of published genomic prediction methods. A prototype Rice Galaxy compliant to Open Access, Open Data, and Findable, Accessible, Interoperable, and Reproducible principles is also presented. Conclusions: Rice Galaxy is a freely available resource that empowers the plant research community to perform state-of-the-art analyses and utilize publicly available big datasets for both fundamental and applied science

Interpreting viral deep sequencing data with GLUE

Using deep sequencing technologies such as Illumina’s platform, it is possible to obtain reads from the viral RNA population revealing the viral genome diversity within a single host. A range of software tools and pipelines can transform raw deep sequencing reads into Sequence Alignment Mapping (SAM) files. We propose that interpretation tools should process these SAM files, directly translating individual reads to amino acids in order to extract statistics of interest such as the proportion of different amino acid residues at specific sites. This preserves per-read linkage between nucleotide variants at different positions within a codon location. The samReporter is a subsystem of the GLUE software toolkit which follows this direct read translation approach in its processing of SAM files. We test samReporter on a deep sequencing dataset obtained from a cohort of 241 UK HCV patients for whom prior treatment with direct-acting antivirals has failed; deep sequencing and resistance testing have been suggested to be of clinical use in this context. We compared the polymorphism interpretation results of the samReporter against an approach that does not preserve per-read linkage. We found that the samReporter was able to properly interpret the sequence data at resistance-associated locations in nine patients where the alternative approach was equivocal. In three cases, the samReporter confirmed that resistance or an atypical substitution was present at NS5A position 30. In three further cases, it confirmed that the sofosbuvir-resistant NS5B substitution S282T was absent. This suggests the direct read translation approach implemented is of value for interpreting viral deep sequencing data

mTOR Controls Ovarian Follicle Growth by Regulating Granulosa Cell Proliferation

We have shown that inhibition of mTOR in granulosa cells and ovarian follicles results in compromised granulosa proliferation and reduced follicle growth. Further analysis here using spontaneously immortalized rat granulosa cells has revealed that mTOR pathway activity is enhanced during M-phase of the cell cycle. mTOR specific phosphorylation of p70S6 kinase and 4E-BP, and expression of Raptor are all enhanced during M-phase. The predominant effect of mTOR inhibition by the specific inhibitor Rapamycin (RAP) was a dose-responsive arrest in the G1 cell cycle stage. The fraction of granulosa cells that continued to divide in the presence of RAP exhibited a dose-dependent increase in aberrant mitotic figures known as anaphase bridges. Strikingly, estradiol consistently decreased the incidence of aberrant mitotic figures. In mice treated with RAP, the mitotic index was reduced compared to controls, and a similar increase in aberrant mitotic events was noted. RAP injected during a superovulation regime resulted in a dose-dependent reduction in the numbers of eggs ovulated. Implications for the real-time regulation of follicle growth and dominance, including the consequences of increased numbers of aneuploid granulosa cells, are discussed

Chromogenic detection of yam mosaic virus by closed-tube reverse transcription loop-mediated isothermal amplification (CT-RT-LAMP)

A closed-tube reverse transcription loop-mediated isothermal amplification (CT-RT-LAMP) assay was developed for the detection of yam mosaic virus (YMV, genus Potyvirus) infecting yam (Dioscorea spp.). The assay uses a set of six oligonucleotide primers targeting the YMV coat protein region, and the amplification products in YMV-positive samples are visualized by chromogenic detection with SYBR Green I dye. The CT-RT-LAMP assay detected YMV in leaf and tuber tissues of infected plants. The assay is 100 times more sensitive in detecting YMV than standard RT-PCR, while maintaining the same specificity

Sublithospheric diamond ages and the supercontinent cycle.

Subduction related to the ancient supercontinent cycle is poorly constrained by mantle samples. Sublithospheric diamond crystallization records the release of melts from subducting oceanic lithosphere at 300-700 km depths1,2 and is especially suited to tracking the timing and effects of deep mantle processes on supercontinents. Here we show that four isotope systems (Rb-Sr, Sm-Nd, U-Pb and Re-Os) applied to Fe-sulfide and CaSiO3 inclusions within 13 sublithospheric diamonds from Juína (Brazil) and Kankan (Guinea) give broadly overlapping crystallization ages from around 450 to 650 million years ago. The intracratonic location of the diamond deposits on Gondwana and the ages, initial isotopic ratios, and trace element content of the inclusions indicate formation from a peri-Gondwanan subduction system. Preservation of these Neoproterozoic-Palaeozoic sublithospheric diamonds beneath Gondwana until its Cretaceous breakup, coupled with majorite geobarometry3,4, suggests that they accreted to and were retained in the lithospheric keel for more than 300 Myr during supercontinent migration. We propose that this process of lithosphere growth-with diamonds attached to the supercontinent keel by the diapiric uprise of depleted buoyant material and pieces of slab crust-could have enhanced supercontinent stability

The Rapid ASKAP Continuum Survey I: Design and First Results

The Rapid ASKAP Continuum Survey (RACS) is the first large-area survey to be conducted with the full 36-antenna Australian Square Kilometre Array Pathfinder (ASKAP) telescope. RACS will provide a shallow model of the ASKAP sky that will aid the calibration of future deep ASKAP surveys. RACS will cover the whole sky visible from the ASKAP site in Western Australia, and will cover the full ASKAP band of $700-1800$ MHz. The RACS images are generally deeper than the existing NRAO VLA Sky Survey (NVSS) and Sydney University Molonglo Sky Survey (SUMSS) radio surveys and have better spatial resolution. All RACS survey products will be public, including radio images (with $\sim 15$ arcsecond resolution) and catalogues of about three million source components with spectral index and polarisation information. In this paper, we present a description of the RACS survey and the first data release of 903 images covering the sky south of declination $+41^\circ$ made over a 288 MHz band centred at 887.5 MHz.Comment: 24 pages, 17 figures, 4 tables. For associated data see https://data.csiro.au/collections/domain/casdaObservation/results/PRAS110%20-%20The%20Rapid%20ASKAP%20Continuu