Ecosystem impacts of three sequential hurricanes (Dennis, Floyd, and Irene) on the United States' largest lagoonal estuary, Pamlico Sound, NC

Three sequential hurricanes, Dennis, Floyd, and Irene, affected coastal North Carolina in September and October 1999. These hurricanes inundated the region with up to 1 m of rainfall, causing 50- to 500-year flooding in the watershed of the Pamlico Sound, the largest lagoonal estuary in the United States and a key West Atlantic fisheries nursery. We investigated the ecosystem-level impacts on and responses of the Sound to the floodwater discharge. Floodwaters displaced three-fourths of the volume of the Sound, depressed salinity by a similar amount, and delivered at least half of the typical annual nitrogen load to this nitrogen-sensitive ecosystem. Organic carbon concentrations in floodwaters entering Pamlico Sound via a major tributary (the Neuse River Estuary) were at least 2-fold higher than concentrations under prefloodwater conditions. A cascading set of physical, chemical, and ecological impacts followed, including strong vertical stratification, bottom water hypoxia, a sustained increase in algal biomass, displacement of many marine organisms, and a rise in fish disease. Because of the Sound's long residence time ( approximately 1 year), we hypothesize that the effects of the short-term nutrient enrichment could prove to be multiannual. A predicted increase in the frequency of hurricane activity over the next few decades may cause longer-term biogeochemical and trophic changes in this and other estuarine and coastal habitats.

Colony-stimulating factor 2 (CSF-2) improves development and posttransfer survival of bovine embryos produced in vitro

In this study, we tested the role of colony-stimulating factor 2 (CSF2) as one of the regulatory molecules that mediate maternal effects on embryonic development during the preimplantation period. Our objective was to verify effects of CSF2 on blastocyst yield, determine posttransfer survival, and evaluate properties of the blastocyst formed after CSF2 treatment. In vitro, CSF2 increased the percentage of oocytes that became morulae and blastocysts. Blastocysts that were treated with CSF2 tended to have a greater number of inner cell mass cells and had a higher ratio of inner cell mass to trophectoderm cells. There was no effect of CSF2 on the incidence of apoptosis. Treatment with CSF2 from d 5 to 7 after insemination increased embryonic survival as indicated by improved pregnancy rate at d 30-35 of gestation. Moreover, treatment with CSF2 from either d 1-7 or 5-7 after insemination reduced pregnancy loss after d 30-35. Results indicate that treatment with CSF2 can affect embryonic development and enhance embryo competence for posttransfer survival. The fact that treatment with CSF2 during such a narrow window of development altered embryonic function much later in pregnancy suggests that CSF2 may exert epigenetic effects on the developing embryo that result in persistent changes in function during the embryonic and fetal periods of development.

Differences in thermoregulatory ability between slick-haired and wild-type lactating Holstein cows in response to acute heat stress

Animals inheriting the slick hair gene have a short, sleek, and sometimes glossy coat. The objective of the present study was to determine whether slick-haired Holstein cows regulate body temperature more effectively than wild-type Holstein cows when exposed to an acute increase in heat stress. Lactating slick cows (n = 10) and wild-type cows (n = 10) were placed for 10 h in an indoor environment with a solid roof, fans, and evaporative cooling or in an outdoor environment with shade cloth and no fans or evaporative cooling. Cows were exposed to both environments in a single reversal design. Vaginal temperature, respiration rate, surface temperature, and sweating rate were measured at 1200, 1500, 1800, and 2100 h (replicate 1) or 1200 and 1500 h (replicate 2), and blood samples were collected for plasma cortisol concentration. Cows in the outdoor environment had higher vaginal and surface temperatures, respiration rates, and sweating rates than cows in the indoor environment. In both environments, slick-haired cows had lower vaginal temperatures (indoor: 39.0 vs. 39.4 degrees C; outdoor 39.6 vs. 40.2 degrees C; SEM = 0.07) and respiration rate (indoor: 67 vs. 79 breaths/ min; outdoor 97 vs. 107 breaths/min; SEM = 5.5) than wild-type cows and greater sweating rates in unclipped areas of skin (indoor: 57 vs. 43 g x h(-1)/m(2); outdoor 82 vs. 61 g x h(-1)/m(2); SEM = 8). Clipping the hair at the site of sweating measurement eliminated the difference between slick-haired and wild-type cows. Results indicate that slick-haired Holstein cows can regulate body temperature more effectively than wild-type cows during heat stress. One reason slick-haired animals are better able to regulate body temperature is increased sweating rate.

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.

Male-specific Fruitless isoforms have different regulatory roles conferred by distinct zinc finger DNA binding domains

Background

Drosophila melanogaster adult males perform an elaborate courtship ritual to entice females to mate. fruitless (fru), a gene that is one of the key regulators of male courtship behavior, encodes multiple male-specific isoforms (Fru^M). These isoforms vary in their carboxy-terminal zinc finger domains, which are predicted to facilitate DNA binding.

Results

By over-expressing individual Fru^M isoforms in fru-expressing neurons in either males or females and assaying the global transcriptional response by RNA-sequencing, we show that three Fru^M isoforms have different regulatory activities that depend on the sex of the fly. We identified several sets of genes regulated downstream of Fru^M isoforms, including many annotated with neuronal functions. By determining the binding sites of individual Fru^M isoforms using SELEX we demonstrate that the distinct zinc finger domain of each Fru^M isoforms confers different DNA binding specificities. A genome-wide search for these binding site sequences finds that the gene sets identified as induced by over-expression of Fru^M isoforms in males are enriched for genes that contain the binding sites. An analysis of the chromosomal distribution of genes downstream of Fru^M shows that those that are induced and repressed in males are highly enriched and depleted on the X chromosome, respectively.

Conclusions

This study elucidates the different regulatory and DNA binding activities of three Fru^M isoforms on a genome-wide scale and identifies genes regulated by these isoforms. These results add to our understanding of sex chromosome biology and further support the hypothesis that in some cell-types genes with male-biased expression are enriched on the X chromosome.

Developmental changes in expression of genes involved in regulation of apoptosis in the bovine preimplantation embryo

The early bovine preimplantation embryo is resistant to proapoptotic signals until around the 8- to 16-cell stage. We hypothesized that 2-cell embryos have higher amounts of antiapoptotic proteins and lower amounts of proapoptotic proteins when compared to embryos ≥16 cells. Steady-state concentrations of mRNA for the antiapoptotic genes BCL2 and HSPA1A were higher for MII oocytes, 2-cell embryos, and 2-cell embryos treated with alpha-amanitin as compared to ≥16-cell embryos. Steady-state concentrations of mRNA for the proapoptotic gene BAD increased in embryos ≥16 cells. There was no significant effect of stage of development on steady-state mRNA concentrations of BCL2L1, DFFA, or BAX. Using immunohistochemistry, it was found that BCL2 was present in greater relative concentrations for 2-cell embryos than for embryos ≥16 cells. These results were confirmed by Western blotting. Relative amounts of immunoreactive BAX detected by immunofluorescence were lower for 2-cell embryos than for embryos ≥16 cells. Using Western blotting, a high molecular weight (46 kDa) form of BAX was highest in ≥16-cell embryos, intermediate in 2-cell embryos, and lowest in MII oocytes. There were no effects of stage of development on relative amounts of immunoreactive BCL2L1, HSPA1A, or BAD, as determined by immunofluorescence. Treatment of embryos with alpha-amanitin from Day 0 to Day 5 or Day 4 to Day 5 after insemination reduced activation of group II caspases and terminal deoxynucleotidyl transferase dUTP nick end labeling after treatment with the proapoptotic signal C(2) ceramide at Day 5 after fertilization. Thus, transcription of BAX or other proteins is required for acquisition of the capacity for apoptosis. Results support the idea that changes in amounts of BCL2 family members are important for the inhibition of apoptosis in the 2-cell embryo and in the establishment of the capacity for apoptosis later in development.

Dynamic sex chromosome expression in Drosophila male germ cells

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

A flexible Bayesian method for detecting allelic imbalance in RNA-seq data

Background

One method of identifying cis regulatory differences is to analyze allele-specific expression (ASE) and identify cases of allelic imbalance (AI). RNA-seq is the most common way to measure ASE and a binomial test is often applied to determine statistical significance of AI. This implicitly assumes that there is no bias in estimation of AI. However, bias has been found to result from multiple factors including: genome ambiguity, reference quality, the mapping algorithm, and biases in the sequencing process. Two alternative approaches have been developed to handle bias: adjusting for bias using a statistical model and filtering regions of the genome suspected of harboring bias. Existing statistical models which account for bias rely on information from DNA controls, which can be cost prohibitive for large intraspecific studies. In contrast, data filtering is inexpensive and straightforward, but necessarily involves sacrificing a portion of the data.

Results

Here we propose a flexible Bayesian model for analysis of AI, which accounts for bias and can be implemented without DNA controls. In lieu of DNA controls, this Poisson-Gamma (PG) model uses an estimate of bias from simulations. The proposed model always has a lower type I error rate compared to the binomial test. Consistent with prior studies, bias dramatically affects the type I error rate. All of the tested models are sensitive to misspecification of bias. The closer the estimate of bias is to the true underlying bias, the lower the type I error rate. Correct estimates of bias result in a level alpha test.

Conclusions

To improve the assessment of AI, some forms of systematic error (e.g., map bias) can be identified using simulation. The resulting estimates of bias can be used to correct for bias in the PG model, without data filtering. Other sources of bias (e.g., unidentified variant calls) can be easily captured by DNA controls, but are missed by common filtering approaches. Consequently, as variant identification improves, the need for DNA controls will be reduced. Filtering does not significantly improve performance and is not recommended, as information is sacrificed without a measurable gain. The PG model developed here performs well when bias is known, or slightly misspecified. The model is flexible and can accommodate differences in experimental design and bias estimation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-920) contains supplementary material, which is available to authorized users.

SECIMTools: a suite of metabolomics data analysis tools

Background

Metabolomics has the promise to transform the area of personalized medicine with the rapid development of high throughput technology for untargeted analysis of metabolites. Open access, easy to use, analytic tools that are broadly accessible to the biological community need to be developed. While technology used in metabolomics varies, most metabolomics studies have a set of features identified. Galaxy is an open access platform that enables scientists at all levels to interact with big data. Galaxy promotes reproducibility by saving histories and enabling the sharing workflows among scientists.

Results

SECIMTools (SouthEast Center for Integrated Metabolomics) is a set of Python applications that are available both as standalone tools and wrapped for use in Galaxy. The suite includes a comprehensive set of quality control metrics (retention time window evaluation and various peak evaluation tools), visualization techniques (hierarchical cluster heatmap, principal component analysis, modular modularity clustering), basic statistical analysis methods (partial least squares - discriminant analysis, analysis of variance, t-test, Kruskal-Wallis non-parametric test), advanced classification methods (random forest, support vector machines), and advanced variable selection tools (least absolute shrinkage and selection operator LASSO and Elastic Net).

Conclusions

SECIMTools leverages the Galaxy platform and enables integrated workflows for metabolomics data analysis made from building blocks designed for easy use and interpretability. Standard data formats and a set of utilities allow arbitrary linkages between tools to encourage novel workflow designs. The Galaxy framework enables future data integration for metabolomics studies with other omics data.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2134-1) contains supplementary material, which is available to authorized users.

Sperm Head-Tail Linkage Requires Restriction of Pericentriolar Material to the Proximal Centriole End

The centriole, or basal body, is the center of attachment between the sperm head and tail. While the distal end of the centriole templates the cilia, the proximal end associates with the nucleus. Using Drosophila, we identify a centriole-centric mechanism that ensures proper proximal end docking to the nucleus. This mechanism relies on the restriction of pericentrin-like protein (PLP) and the pericentriolar material (PCM) to the proximal end of the centriole. PLP is restricted proximally by limiting its mRNA and protein to the earliest stages of centriole elongation. Ectopic positioning of PLP to more distal portions of the centriole is sufficient to redistribute PCM and microtubules along the entire centriole length. This results in erroneous, lateral centriole docking to the nucleus, leading to spermatid decapitation as a result of a failure to form a stable head-tail linkage.

The Wright stuff: reimagining path analysis reveals novel components of the sex determination hierarchy in Drosophila melanogaster

BACKGROUND

The Drosophila sex determination hierarchy is a classic example of a transcriptional regulatory hierarchy, with sex-specific isoforms regulating morphology and behavior. We use a structural equation modeling approach, leveraging natural genetic variation from two studies on Drosophila female head tissues--DSPR collection (596 F1-hybrids from crosses between DSPR sub-populations) and CEGS population (75 F1-hybrids from crosses between DGRP/Winters lines to a reference strain w1118)--to expand understanding of the sex hierarchy gene regulatory network (GRN). This approach is completely generalizable to any natural population, including humans.

RESULTS

We expanded the sex hierarchy GRN adding novel links among genes, including a link from fruitless (fru) to Sex-lethal (Sxl) identified in both populations. This link is further supported by the presence of fru binding sites in the Sxl locus. 754 candidate genes were added to the pathway, including the splicing factors male-specific lethal 2 and Rm62 as downstream targets of Sxl which are well-supported links in males. Independent studies of doublesex and transformer mutants support many additions, including evidence for a link between the sex hierarchy and metabolism, via Insulin-like receptor.

CONCLUSIONS

The genes added in the CEGS population were enriched for genes with sex-biased splicing and components of the spliceosome. A common goal of molecular biologists is to expand understanding about regulatory interactions among genes. Using natural alleles we can not only identify novel relationships, but using supervised approaches can order genes into a regulatory hierarchy. Combining these results with independent large effect mutation studies, allows clear candidates for detailed molecular follow-up to emerge.

Reliable Detection of Herpes Simplex Virus Sequence Variation by High-Throughput Resequencing

High-throughput sequencing (HTS) has resulted in data for a number of herpes simplex virus (HSV) laboratory strains and clinical isolates. The knowledge of these sequences has been critical for investigating viral pathogenicity. However, the assembly of complete herpesviral genomes, including HSV, is complicated due to the existence of large repeat regions and arrays of smaller reiterated sequences that are commonly found in these genomes. In addition, the inherent genetic variation in populations of isolates for viruses and other microorganisms presents an additional challenge to many existing HTS sequence assembly pipelines. Here, we evaluate two approaches for the identification of genetic variants in HSV1 strains using Illumina short read sequencing data. The first, a reference-based approach, identifies variants from reads aligned to a reference sequence and the second, a de novo assembly approach, identifies variants from reads aligned to de novo assembled consensus sequences. Of critical importance for both approaches is the reduction in the number of low complexity regions through the construction of a non-redundant reference genome. We compared variants identified in the two methods. Our results indicate that approximately 85% of variants are identified regardless of the approach. The reference-based approach to variant discovery captures an additional 15% representing variants divergent from the HSV1 reference possibly due to viral passage. Reference-based approaches are significantly less labor-intensive and identify variants across the genome where de novo assembly-based approaches are limited to regions where contigs have been successfully assembled. In addition, regions of poor quality assembly can lead to false variant identification in de novo consensus sequences. For viruses with a well-assembled reference genome, a reference-based approach is recommended.

NetREX-CF integrates incomplete transcription factor data with gene expression to reconstruct gene regulatory networks

The inference of Gene Regulatory Networks (GRNs) is one of the key challenges in systems biology. Leading algorithms utilize, in addition to gene expression, prior knowledge such as Transcription Factor (TF) DNA binding motifs or results of TF binding experiments. However, such prior knowledge is typically incomplete, therefore, integrating it with gene expression to infer GRNs remains difficult. To address this challenge, we introduce NetREX-CF-Regulatory Network Reconstruction using EXpression and Collaborative Filtering-a GRN reconstruction approach that brings together Collaborative Filtering to address the incompleteness of the prior knowledge and a biologically justified model of gene expression (sparse Network Component Analysis based model). We validated the NetREX-CF using Yeast data and then used it to construct the GRN for Drosophila Schneider 2 (S2) cells. To corroborate the GRN, we performed a large-scale RNA-Seq analysis followed by a high-throughput RNAi treatment against all 465 expressed TFs in the cell line. Our knockdown result has not only extensively validated the GRN we built, but also provides a benchmark that our community can use for evaluating GRNs. Finally, we demonstrate that NetREX-CF can infer GRNs using single-cell RNA-Seq, and outperforms other methods, by using previously published human data.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.