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.

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.

Targeted Mutagenesis of the Multicopy Myrosinase Gene Family in Allotetraploid Brassica juncea Reduces Pungency in Fresh Leaves across Environments.

Recent breeding efforts in Brassica have focused on the development of new oilseed feedstock crop for biofuels (e.g., ethanol, biodiesel, bio-jet fuel), bio-industrial uses (e.g., bio-plastics, lubricants), specialty fatty acids (e.g., erucic acid), and producing low glucosinolates levels for oilseed and feed meal production for animal consumption. We identified a novel opportunity to enhance the availability of nutritious, fresh leafy greens for human consumption. Here, we demonstrated the efficacy of disarming the 'mustard bomb' reaction in reducing pungency upon the mastication of fresh tissue-a major source of unpleasant flavor and/or odor in leafy Brassica. Using gene-specific mutagenesis via CRISPR-Cas12a, we created knockouts of all functional copies of the type-I myrosinase multigene family in tetraploid Brassica juncea. Our greenhouse and field trials demonstrate, via sensory and biochemical analyses, a stable reduction in pungency in edited plants across multiple environments. Collectively, these efforts provide a compelling path toward boosting the human consumption of nutrient-dense, fresh, leafy green vegetables.

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.

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.

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.

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.

Sex Differences in Drosophila Somatic Gene Expression: Variation and Regulation by doublesex.

Sex differences in gene expression have been widely studied in Drosophila melanogaster Sex differences vary across strains, but many molecular studies focus on only a single strain, or on genes that show sexually dimorphic expression in many strains. How extensive variability is and whether this variability occurs among genes regulated by sex determination hierarchy terminal transcription factors is unknown. To address these questions, we examine differences in sexually dimorphic gene expression between two strains in Drosophila adult head tissues. We also examine gene expression in doublesex (dsx) mutant strains to determine which sex-differentially expressed genes are regulated by DSX, and the mode by which DSX regulates expression. We find substantial variation in sex-differential expression. The sets of genes with sexually dimorphic expression in each strain show little overlap. The prevalence of different DSX regulatory modes also varies between the two strains. Neither the patterns of DSX DNA occupancy, nor mode of DSX regulation explain why some genes show consistent sex-differential expression across strains. We find that the genes identified as regulated by DSX in this study are enriched with known sites of DSX DNA occupancy. Finally, we find that sex-differentially expressed genes and genes regulated by DSX are highly enriched on the fourth chromosome. These results provide insights into a more complete pool of potential DSX targets, as well as revealing the molecular flexibility of DSX regulation.

Buffering of Genetic Regulatory Networks in Drosophila melanogaster.

Regulatory variation in gene expression can be described by cis- and trans-genetic components. Here we used RNA-seq data from a population panel of Drosophila melanogaster test crosses to compare allelic imbalance (AI) in female head tissue between mated and virgin flies, an environmental change known to affect transcription. Indeed, 3048 exons (1610 genes) are differentially expressed in this study. A Bayesian model for AI, with an intersection test, controls type I error. There are ∼200 genes with AI exclusively in mated or virgin flies, indicating an environmental component of expression regulation. On average 34% of genes within a cross and 54% of all genes show evidence for genetic regulation of transcription. Nearly all differentially regulated genes are affected in cis, with an average of 63% of expression variation explained by the cis-effects. Trans-effects explain 8% of the variance in AI on average and the interaction between cis and trans explains an average of 11% of the total variance in AI. In both environments cis- and trans-effects are compensatory in their overall effect, with a negative association between cis- and trans-effects in 85% of the exons examined. We hypothesize that the gene expression level perturbed by cis-regulatory mutations is compensated through trans-regulatory mechanisms, e.g., trans and cis by trans-factors buffering cis-mutations. In addition, when AI is detected in both environments, cis-mated, cis-virgin, and trans-mated-trans-virgin estimates are highly concordant with 99% of all exons positively correlated with a median correlation of 0.83 for cis and 0.95 for trans We conclude that the gene regulatory networks (GRNs) are robust and that trans-buffering explains robustness.

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.

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.

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.

Cheating death at the dawn of life: developmental control of apoptotic repression in the preimplantation embryo.

During early development, the mammalian embryo is resistant to pro-apoptotic signals because of biochemical properties of the mitochondrion and nucleus. Mitochondria of the bovine two-cell embryo are resistant to depolarization because of low amounts of the proapoptotic protein BAX and high concentrations of the anti-apoptotic protein BCL2. As development proceeds, BAX content increases, BCL2 content declines, and mitochondria becomes capable of pore formation and depolarization in response to pro-apoptotic signals. The nucleus of the two-cell embryo is resistant to degradation by the DNase DFFB because epigenetic modifications, including DNA methylation and histone deacetylation, mask internucleosomal sites for DNA cleavage. Blastomere DNA becomes progressively less methylated during development so that DNA becomes accessible to cleavage by DFFB.

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.

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.