In-depth mechanistic analysis including high-throughput RNA sequencing in the prediction of functional and structural cardiotoxicants using hiPSC cardiomyocytes.

BACKGROUND: Cardiotoxicity remains one of the most reported adverse drug reactions that lead to drug attrition during pre-clinical and clinical drug development. Drug-induced cardiotoxicity may develop as a functional change in cardiac electrophysiology (acute alteration of the mechanical function of the myocardium) and/or as a structural change, resulting in loss of viability and morphological damage to cardiac tissue. RESEARCH DESIGN AND METHODS: Non-clinical models with better predictive value need to be established to improve cardiac safety pharmacology. To this end, high-throughput RNA sequencing (ScreenSeq) was combined with high-content imaging (HCI) and Ca2+ transience (CaT) to analyze compound-treated human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). RESULTS: Analysis of hiPSC-CMs treated with 33 cardiotoxicants and 9 non-cardiotoxicants of mixed therapeutic indications facilitated compound clustering by mechanism of action, scoring of pathway activities related to cardiomyocyte contractility, mitochondrial integrity, metabolic state, diverse stress responses and the prediction of cardiotoxicity risk. The combination of ScreenSeq, HCI and CaT provided a high cardiotoxicity prediction performance with 89% specificity, 91% sensitivity and 90% accuracy. CONCLUSIONS: Overall, this study introduces mechanism-driven risk assessment approach combining structural, functional and molecular high-throughput methods for pre-clinical risk assessment of novel compounds.

Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data.

The Drosophila virilis group is one of the major lineages of Drosophila previously recognised and it has been used as a model for different types of studies. It comprises 13 species whose phylogenetic relationships are not well resolved. In the present study, six nuclear genes (Adh, fused, Gpdh, NonA, CG9631 and CG7219) and the mitochondrial ribosomal RNA genes (12S-16S) have been used to estimate the evolutionary tree of the group using different methods of phylogenetic reconstruction. Different competing evolutionary hypotheses have also been compared using the Approximately Unbiased test to further evaluate the robustness of the inferred trees. Results are, in general, consistent with previous studies in recovering the four major lineages of the group (D. virilis phylad, Drosophila montana subphylad, Drosophila kanekoi subphylad and Drosophila littoralis subphylad), although D. kanekoi, D. littoralis and Drosophila ezoana are here inferred to be more closely related to the D. virilis phylad than to the D. montana subphylad. The age of the crown group, estimated with a Bayesian method that assumes a relaxed molecular clock, is placed in the late Miocene (∼ 10 Mya). The oldest lineages also appeared during this period (∼ 7.5 to ∼ 8.9 Mya), while the ages of the basal nodes of the montana subphylad and the virilis phylad are located in the early Pliocene (∼ 4.9 and ∼ 4.1 Mya). Major cladogenesis events correlate to geological and palaeoclimatic occurrences that most likely affected the freshwater and deciduous forests where these species are found. The inferred biogeographical history of the group, based on the statistical dispersal-vicariance analysis, indicates that the last common ancestor of the group had a Holarctic distribution from which the North American and the Eurasian lineages evolved as a result of a vicariant event.

Multiple loci linked to inversions are associated with eye size variation in species of the Drosophila virilis phylad.

The size and shape of organs is tightly controlled to achieve optimal function. Natural morphological variations often represent functional adaptations to an ever-changing environment. For instance, variation in head morphology is pervasive in insects and the underlying molecular basis is starting to be revealed in the Drosophila genus for species of the melanogaster group. However, it remains unclear whether similar diversifications are governed by similar or different molecular mechanisms over longer timescales. To address this issue, we used species of the virilis phylad because they have been diverging from D. melanogaster for at least 40 million years. Our comprehensive morphological survey revealed remarkable differences in eye size and head shape among these species with D. novamexicana having the smallest eyes and southern D. americana populations having the largest eyes. We show that the genetic architecture underlying eye size variation is complex with multiple associated genetic variants located on most chromosomes. Our genome wide association study (GWAS) strongly suggests that some of the putative causative variants are associated with the presence of inversions. Indeed, northern populations of D. americana share derived inversions with D. novamexicana and they show smaller eyes compared to southern ones. Intriguingly, we observed a significant enrichment of genes involved in eye development on the 4th chromosome after intersecting chromosomal regions associated with phenotypic differences with those showing high differentiation among D. americana populations. We propose that variants associated with chromosomal inversions contribute to both intra- and interspecific variation in eye size among species of the virilis phylad.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies.

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

An old bilbo-like non-LTR retroelement insertion provides insight into the relationship of species of the virilis group.

In Drosophila, at the population and species level, fixation of a TE insertion is an unlikely fate. Of the few reported fixations at the species level most involve non-LTR retroelements. In this work we report the fixation of a non-LTR retroelement in five species (Drosophila littoralis, Drosophila virilis, Drosophila lummei, Drosophila americana and Drosophila novamexicana) of the virilis group of Drosophila. In most species, this TE insertion is being lost through the accumulation of small deletions, but there is also evidence for the accumulation of large deletions. In the americana lineage an insertion of about 900 bp of the non-LTR retroelement is a marker for the Xc inversion. This insertion is, at most, 80 kb away from the basal Xc inversion breakpoint. The presence of a bilbo-like element in D. littoralis but not in D. kanekoi, suggests that D. littoralis is more closely related to species of the virilis phylad than to species of the montana phylad, which is in contrast with the traditional view. Nevertheless, the phylogenetic analyses here performed using a 7 gene dataset suggest that D. littoralis is indeed more closely related to species of the virilis phylad than to species of the montana phylad. The re-evaluation of the phylogenetic relationship of the species of the virilis group, under the assumption of a relaxed molecular clock, results in an estimated age of the bilbo-like element insertion of at least 7.5 Mya.

Origin and Consequences of Chromosomal Inversions in the virilis Group of Drosophila.

In Drosophila, large variations in rearrangement rate have been reported among different lineages and among Muller's elements. Nevertheless, the mechanisms that are involved in the generation of inversions, their increase in frequency, as well as their impact on the genome are not completely understood. This is in part due to the lack of comparative studies on species distantly related to Drosophila melanogaster. Therefore, we sequenced and assembled the genomes of two species of the virilis phylad (Drosophila novamexicana [15010-1031.00] and Drosophila americana [SF12]), which are diverging from D. melanogaster for more than 40 Myr. Based on these data, we identified the precise location of six novel inversion breakpoints. A molecular characterization provided clear evidence that DAIBAM (a miniature inverted-repeat transposable element) was involved in the generation of eight out of the nine inversions identified. In contrast to what has been previously reported for D. melanogaster and close relatives, ectopic recombination is thus the prevalent mechanism of generating inversions in species of the virilis phylad. Using pool-sequencing data for three populations of D. americana, we also show that common polymorphic inversions create a high degree of genetic differentiation between populations for chromosomes X, 4, and 5 over large physical distances. We did not find statistically significant differences in expression levels between D. americana (SF12) and D. novamexicana (15010-1031.00) strains for the three genes surveyed (CG9588, Fig 4, and fab1) flanking three inversion breakpoints.

The Drosophila melanogaster Muc68E Mucin Gene Influences Adult Size, Starvation Tolerance, and Cold Recovery.

Mucins have been implicated in many different biological processes, such as protection from mechanical damage, microorganisms, and toxic molecules, as well as providing a luminal scaffold during development. Nevertheless, it is conceivable that mucins have the potential to modulate food absorption as well, and thus contribute to the definition of several important phenotypic traits. Here we show that the Drosophila melanogaster Muc68E gene is 40- to 60-million-yr old, and is present in Drosophila species of the subgenus Sophophora only. The central repeat region of this gene is fast evolving, and shows evidence for repeated expansions/contractions. This and/or frequent gene conversion events lead to the homogenization of its repeats. The amino acid pattern P[ED][ED][ST][ST][ST] is found in the repeat region of Muc68E proteins from all Drosophila species studied, and can occur multiple times within a single conserved repeat block, and thus may have functional significance. Muc68E is a nonessential gene under laboratory conditions, but Muc68E mutant flies are smaller and lighter than controls at birth. However, at 4 d of age, Muc68E mutants are heavier, recover faster from chill-coma, and are more resistant to starvation than control flies, although they have the same percentage of lipids as controls. Mutant flies have enlarged abdominal size 1 d after chill-coma recovery, which is associated with higher lipid content. These results suggest that Muc68E has a role in metabolism modulation, food absorption, and/or feeding patterns in larvae and adults, and under normal and stress conditions. Such biological function is novel for mucin genes.

Drosophila americana Diapausing Females Show Features Typical of Young Flies.

Diapause is a period of arrested development which is controlled physiologically, preprogrammed environmentally and characterized by metabolic depression that can occur during any stage of insect development. Nevertheless, in the genus Drosophila, diapause is almost always associated with the cessation of ovarian development and reproductive activity in adult females. In this work, we show that, in D. americana (a temperate species of the virilis group), diapause is a genetically determined delay in ovarian development that is triggered by temperature and/or photoperiod. Moreover, we show that in this species diapause incidence increases with latitude, ranging from 13% in the southernmost to 91% in the northernmost range of the distribution. When exposed to diapause inducing conditions, both diapausing and non-diapausing females show a 10% increase in lifespan, that is further increased by 18.6% in diapausing females, although senescence is far from being negligible.ActinD1 expression levels suggest that diapausing females are biologically much younger than their chronological age, and that the fly as a whole, rather than the ovarian developmental one, which is phenotypically more evident, is delayed by diapause. Therefore, diapause candidate genes that show expression levels that are compatible with flies younger than their chronological age may not necessarily play a role in reproductive diapause and in adaptation to seasonally varying environmental conditions [corrected].

Genes belonging to the insulin and ecdysone signaling pathways can contribute to developmental time, lifespan and abdominal size variation in Drosophila americana.

Even within a single genus, such as Drosophila, cases of lineage-specific adaptive evolution have been found. Therefore, the molecular basis of phenotypic variation must be addressed in more than one species group, in order to infer general patterns. In this work, we used D. americana, a species distantly-related to D. melanogaster, to perform an F2 association study for developmental time (DT), chill-coma recovery time (CRT), abdominal size (AS) and lifespan (LS) involving the two strains (H5 and W11) whose genomes have been previously sequenced. Significant associations were found between the 43 large indel markers developed here and DT, AS and LS but not with CRT. Significant correlations are also found between DT and LS, and between AS and LS, that might be explained by variation at genes belonging to the insulin and ecdysone signaling pathways. Since, in this F2 association study a single marker, located close to the Ecdysone receptor (EcR) gene, explained as much as 32.6% of the total variation in DT, we performed a second F2 association study, to determine whether large differences in DT are always due to variation in this genome region. No overlapping signal was observed between the two F2 association studies. Overall, these results illustrate that, in D. americana, pleiotropic genes involved in the highly-conserved insulin and ecdysone signaling pathways are likely responsible for variation observed in ecologically relevant phenotypic traits, although other genes are also involved.

Drosophila americana as a model species for comparative studies on the molecular basis of phenotypic variation.

Understanding the molecular basis of within and between species phenotypic variation is one of the main goals of Biology. In Drosophila, most of the work regarding this issue has been performed in D. melanogaster, but other distantly related species must also be studied to verify the generality of the findings obtained for this species. Here, we make the case for D. americana, a species of the virilis group of Drosophila that has been diverging from the model species, D. melanogaster, for approximately 40 Myr. To determine the suitability of this species for such studies, polymorphism and recombination estimates are presented for D. americana based on the largest nucleotide sequence polymorphism data set so far analyzed (more than 100 data sets) for this species. The polymorphism estimates are also compared with those obtained from the comparison of the genome assembly of two D. americana strains (H5 and W11) here reported. As an example of the general utility of these resources, we perform a preliminary study on the molecular basis of lifespan differences in D. americana. First, we show that there are lifespan differences between D. americana populations from different regions of the distribution range. Then, we perform five F2 association experiments using markers for 21 candidate genes previously identified in D. melanogaster. Significant associations are found between polymorphism at two genes (hep and Lim3) and lifespan. For the F2 association study involving the two sequenced strains (H5 and W11), we identify amino acid differences at Lim3 and Hep that could be responsible for the observed changes in lifespan. For both genes, no large gene expression differences were observed between the two strains.

The Drosophila melanogaster methuselah gene: a novel gene with ancient functions.

The Drosophila melanogaster G protein-coupled receptor gene, methuselah (mth), has been described as a novel gene that is less than 10 million years old. Nevertheless, it shows a highly specific expression pattern in embryos, larvae, and adults, and has been implicated in larval development, stress resistance, and in the setting of adult lifespan, among others. Although mth belongs to a gene subfamily with 16 members in D. melanogaster, there is no evidence for functional redundancy in this subfamily. Therefore, it is surprising that a novel gene influences so many traits. Here, we explore the alternative hypothesis that mth is an old gene. Under this hypothesis, in species distantly related to D. melanogaster, there should be a gene with features similar to those of mth. By performing detailed phylogenetic, synteny, protein structure, and gene expression analyses we show that the D. virilis GJ12490 gene is the orthologous of mth in species distantly related to D. melanogaster. We also show that, in D. americana (a species of the virilis group of Drosophila), a common amino acid polymorphism at the GJ12490 orthologous gene is significantly associated with developmental time, size, and lifespan differences. Our results imply that GJ12490 orthologous genes are candidates for developmental time and lifespan differences in Drosophila in general.

Drosophila genes that affect meiosis duration are among the meiosis related genes that are more often found duplicated.

Using a phylogenetic approach, the examination of 33 meiosis/meiosis-related genes in 12 Drosophila species, revealed nine independent gene duplications, involving the genes cav, mre11, meiS332, polo and mtrm. Evidence is provided that at least eight out of the nine gene duplicates are functional. Therefore, the rate at which Drosophila meiosis/meiosis-related genes are duplicated and retained is estimated to be 0.0012 per gene per million years, a value that is similar to the average for all Drosophila genes. It should be noted that by using a phylogenetic approach the confounding effect of concerted evolution, that is known to lead to overestimation of the duplication and retention rate, is avoided. This is an important issue, since even in our moderate size sample, evidence for long-term concerted evolution (lasting for more than 30 million years) was found for the meiS332 gene pair in species of the Drosophila subgenus. Most striking, in contrast to theoretical expectations, is the finding that genes that encode proteins that must follow a close stoichiometric balance, such as polo, mtrm and meiS332 have been found duplicated. The duplicated genes may be examples of gene neofunctionalization. It is speculated that meiosis duration may be a trait that is under selection in Drosophila and that it has different optimal values in different species.

A comparative study of the short term cold resistance response in distantly related Drosophila species: the role of regucalcin and frost.

The molecular basis of short term cold resistance (indexed as chill-coma recovery time) has been mostly addressed in D. melanogaster, where candidate genes (Dca (also known as smp-30) and Frost (Fst)) have been identified. Nevertheless, in Drosophila, the ability to tolerate short term exposure to low temperatures evolved several times independently. Therefore, it is unclear whether variation in the same candidate genes is also responsible for short term cold resistance in distantly related Drosophila species. It should be noted that Dca is a candidate gene for cold resistance in the Sophophora subgenus only, since there is no orthologous gene copy in the Drosophila subgenus. Here we show that, in D. americana (Drosophila subgenus), there is a north-south gradient for a variant at the 5' non-coding region of regucalcin (a Dca-like gene; in D. melanogaster the proteins encoded by the two genes share 71.9% amino acid identities) but in our D. americana F2 association experiment there is no association between this polymorphism and chill-coma recovery times. Moreover, we found no convincing evidence that this gene is up-regulated after cold shock in both D. americana and D. melanogaster. Size variation in the Fst PEST domain (putatively involved in rapid protein degradation) is observed when comparing distantly related Drosophila species, and is associated with short term cold resistance differences in D. americana. Nevertheless, this effect is likely through body size variation. Moreover, we show that, even at two hours after cold shock, when up-regulation of this gene is maximal in D. melanogaster (about 48 fold expression change), in D. americana this gene is only moderately up-regulated (about 3 fold expression change). Our work thus shows that there are important differences regarding the molecular basis of cold resistance in distantly related Drosophila species.