Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation.

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

A SARS-CoV-2 targeted siRNA-nanoparticle therapy for COVID-19.

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in humans. Despite several emerging vaccines, there remains no verifiable therapeutic targeted specifically to the virus. Here we present a highly effective small interfering RNA (siRNA) therapeutic against SARS-CoV-2 infection using a novel lipid nanoparticle (LNP) delivery system. Multiple siRNAs targeting highly conserved regions of the SARS-CoV-2 virus were screened, and three candidate siRNAs emerged that effectively inhibit the virus by greater than 90% either alone or in combination with one another. We simultaneously developed and screened two novel LNP formulations for the delivery of these candidate siRNA therapeutics to the lungs, an organ that incurs immense damage during SARS-CoV-2 infection. Encapsulation of siRNAs in these LNPs followed by in vivo injection demonstrated robust repression of virus in the lungs and a pronounced survival advantage to the treated mice. Our LNP-siRNA approaches are scalable and can be administered upon the first sign of SARS-CoV-2 infection in humans. We suggest that an siRNA-LNP therapeutic approach could prove highly useful in treating COVID-19 disease as an adjunctive therapy to current vaccine strategies.

Inhibition of CD38 and supplementation of nicotinamide riboside ameliorate lipopolysaccharide-induced microglial and astrocytic neuroinflammation by increasing NAD.

Neuroinflammation is initiated by activation of the brain's innate immune system in response to an inflammatory challenge. Insufficient control of neuroinflammation leads to enhanced or prolonged pathology in various neurological conditions including multiple sclerosis and Alzheimer's disease. Nicotinamide adenine dinucleotide (NAD+ ) plays critical roles in cellular energy metabolism and calcium homeostasis. Our previous study demonstrated that deletion of CD38, which consumes NAD+ , suppressed cuprizone-induced demyelination, neuroinflammation, and glial activation. However, it is still unknown whether CD38 directly affects neuroinflammation through regulating brain NAD+ level. In this study, we investigated the effect of CD38 deletion and inhibition and supplementation of NAD+ on lipopolysaccharide (LPS)-induced neuroinflammation in mice. Intracerebroventricular injection of LPS significantly increased CD38 expression especially in the hippocampus. Deletion of CD38 decreased LPS-induced inflammatory responses and glial activation. Pre-administration of apigenin, a flavonoid with CD38 inhibitory activity, or nicotinamide riboside (NR), an NAD+ precursor, increased NAD+ level, and significantly suppressed induction of cytokines and chemokines, glial activation and subsequent neurodegeneration after LPS administration. In cell culture, LPS-induced inflammatory responses were suppressed by treatment of primary astrocytes or microglia with apigenin, NAD+ , NR or 78c, the latter a specific CD38 inhibitor. Finally, all these compounds suppressed NF-κB signaling pathway in microglia. These results suggest that CD38-mediated neuroinflammation is linked to NAD+ consumption and that boosting NAD+ by CD38 inhibition and NR supplementation directly suppress neuroinflammation in the brain.

Sex reversal triggers the rapid transition from genetic to temperature-dependent sex.

Sex determination in animals is amazingly plastic. Vertebrates display contrasting strategies ranging from complete genetic control of sex (genotypic sex determination) to environmentally determined sex (for example, temperature-dependent sex determination). Phylogenetic analyses suggest frequent evolutionary transitions between genotypic and temperature-dependent sex determination in environmentally sensitive lineages, including reptiles. These transitions are thought to involve a genotypic system becoming sensitive to temperature, with sex determined by gene-environment interactions. Most mechanistic models of transitions invoke a role for sex reversal. Sex reversal has not yet been demonstrated in nature for any amniote, although it occurs in fish and rarely in amphibians. Here we make the first report of reptile sex reversal in the wild, in the Australian bearded dragon (Pogona vitticeps), and use sex-reversed animals to experimentally induce a rapid transition from genotypic to temperature-dependent sex determination. Controlled mating of normal males to sex-reversed females produces viable and fertile offspring whose phenotypic sex is determined solely by temperature (temperature-dependent sex determination). The W sex chromosome is eliminated from this lineage in the first generation. The instantaneous creation of a lineage of ZZ temperature-sensitive animals reveals a novel, climate-induced pathway for the rapid transition between genetic and temperature-dependent sex determination, and adds to concern about adaptation to rapid global climate change.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction.

BACKGROUND: Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes. RESULTS: We generated 96.7 × 109 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 109 paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females. CONCLUSIONS: The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Waking the sleeping dragon: gene expression profiling reveals adaptive strategies of the hibernating reptile Pogona vitticeps.

BACKGROUND: Hibernation is a physiological state exploited by many animals exposed to prolonged adverse environmental conditions associated with winter. Large changes in metabolism and cellular function occur, with many stress response pathways modulated to tolerate physiological challenges that might otherwise be lethal. Many studies have sought to elucidate the molecular mechanisms of mammalian hibernation, but detailed analyses are lacking in reptiles. Here we examine gene expression in the Australian central bearded dragon (Pogona vitticeps) using mRNA-seq and label-free quantitative mass spectrometry in matched brain, heart and skeletal muscle samples from animals at late hibernation, 2 days post-arousal and 2 months post-arousal. RESULTS: We identified differentially expressed genes in all tissues between hibernation and post-arousal time points; with 4264 differentially expressed genes in brain, 5340 differentially expressed genes in heart, and 5587 differentially expressed genes in skeletal muscle. Furthermore, we identified 2482 differentially expressed genes across all tissues. Proteomic analysis identified 743 proteins (58 differentially expressed) in brain, 535 (57 differentially expressed) in heart, and 337 (36 differentially expressed) in skeletal muscle. Tissue-specific analyses revealed enrichment of protective mechanisms in all tissues, including neuroprotective pathways in brain, cardiac hypertrophic processes in heart, and atrophy protective pathways in skeletal muscle. In all tissues stress response pathways were induced during hibernation, as well as evidence for gene expression regulation at transcription, translation and post-translation. CONCLUSIONS: These results reveal critical stress response pathways and protective mechanisms that allow for maintenance of both tissue-specific function, and survival during hibernation in the central bearded dragon. Furthermore, we provide evidence for multiple levels of gene expression regulation during hibernation, particularly enrichment of miRNA-mediated translational repression machinery; a process that would allow for rapid and energy efficient reactivation of translation from mature mRNA molecules at arousal. This study is the first molecular investigation of its kind in a hibernating reptile, and identifies strategies not yet observed in other hibernators to cope stress associated with this remarkable state of metabolic depression.

Landscape of DNA Methylation on the Marsupial X.

DNA methylation plays a key role in maintaining transcriptional silence on the inactive X chromosome of eutherian mammals. Beyond eutherians, there are limited genome wide data on DNA methylation from other vertebrates. Previous studies of X borne genes in various marsupial models revealed no differential DNA methylation of promoters between the sexes, leading to the conclusion that CpG methylation plays no role in marsupial X-inactivation. Using reduced representation bisulfite sequencing, we generated male and female CpG methylation profiles in four representative vertebrates (mouse, gray short-tailed opossum, platypus, and chicken). A variety of DNA methylation patterns were observed. Platypus and chicken displayed no large-scale differential DNA methylation between the sexes on the autosomes or the sex chromosomes. As expected, a metagene analysis revealed hypermethylation at transcription start sites (TSS) of genes subject to X-inactivation in female mice. This contrasted with the opossum, in which metagene analysis did not detect differential DNA methylation between the sexes at TSSs of genes subject to X-inactivation. However, regions flanking TSSs of these genes were hypomethylated. Our data are the first to demonstrate that, for genes subject to X-inactivation in both eutherian and marsupial mammals, there is a consistent difference between DNA methylation levels at TSSs and immediate flanking regions, which we propose has a silencing effect in both groups.

Transcriptome sequencing of the long-nosed bandicoot (Perameles nasuta) reveals conservation and innovation of immune genes in the marsupial order Peramelemorphia.

Bandicoots are omnivorous marsupials of the order Peramelemorphia. Conservation concerns and their unique biological characteristics suggest peramelomorphs are worthy research subjects, but knowledge of their genetics and immunology has lagged behind that of other high-profile marsupials. Here, we characterise the transcriptome of the long-nose bandicoot (Perameles nasuta), the first high-throughput data set from any peramelomorph. We investigate the immune gene repertoire of the bandicoot, with a focus on key immune gene families, and compare to previously characterised marsupial and mammalian species. We find that the immune gene complement in bandicoot is often conserved with respect to other marsupials; however, the diversity of expressed transcripts in several key families, such as major histocompatibility complex, T cell receptor µ and natural killer cell receptors, appears greater in the bandicoot than other Australian marsupials, including devil and koala. This transcriptome is an important first step for future studies of bandicoots and the bilby, allowing for population level analysis and construction of bandicoot-specific immunological reagents and assays. Such studies will be critical to understanding the immunology and physiology of Peramelemorphia and to inform the conservation of these unique marsupials.

Characterisation of MHC class I genes in the koala.

Koala (Phascolarctos cinereus) populations are on the decline across the majority of Australia's mainland. Two major diseases threatening the long-term survival of affected koala populations are caused by obligate intracellular pathogens: Chlamydia and koala retrovirus (KoRV). To improve our understanding of the koala immune system, we characterised their major histocompatibility complex (MHC) class I genes, which are centrally involved in presenting foreign peptides derived from intracellular pathogens to cytotoxic T cells. A total of 11 class I genes were identified in the koala genome. Three genes, Phci-UA, UB and UC, showed relatively high genetic variability and were expressed in all 12 examined tissues, whereas the other eight genes had tissue-specific expression and limited polymorphism. Evidence of diversifying selection was detected in Phci-UA and UC, while gene conversion may have played a role in creating new alleles at Phci-UB. We propose that Phci-UA, UB and UC are likely classical MHC genes of koalas, and further research is needed to understand their role in koala chlamydial and KoRV infections.

Widespread convergence in toxin resistance by predictable molecular evolution.

The question about whether evolution is unpredictable and stochastic or intermittently constrained along predictable pathways is the subject of a fundamental debate in biology, in which understanding convergent evolution plays a central role. At the molecular level, documented examples of convergence are rare and limited to occurring within specific taxonomic groups. Here we provide evidence of constrained convergent molecular evolution across the metazoan tree of life. We show that resistance to toxic cardiac glycosides produced by plants and bufonid toads is mediated by similar molecular changes to the sodium-potassium-pump (Na(+)/K(+)-ATPase) in insects, amphibians, reptiles, and mammals. In toad-feeding reptiles, resistance is conferred by two point mutations that have evolved convergently on four occasions, whereas evidence of a molecular reversal back to the susceptible state in varanid lizards migrating to toad-free areas suggests that toxin resistance is maladaptive in the absence of selection. Importantly, resistance in all taxa is mediated by replacements of 2 of the 12 amino acids comprising the Na(+)/K(+)-ATPase H1-H2 extracellular domain that constitutes a core part of the cardiac glycoside binding site. We provide mechanistic insight into the basis of resistance by showing that these alterations perturb the interaction between the cardiac glycoside bufalin and the Na(+)/K(+)-ATPase. Thus, similar selection pressures have resulted in convergent evolution of the same molecular solution across the breadth of the animal kingdom, demonstrating how a scarcity of possible solutions to a selective challenge can lead to highly predictable evolutionary responses.

Differential Gamma Interferon- and Tumor Necrosis Factor Alpha-Driven Cytokine Response Distinguishes Acute Infection of a Metatherian Host with Toxoplasma gondii and Neospora caninum.

Toxoplasma gondii and Neospora caninum (both Apicomplexa) are closely related cyst-forming coccidian parasites that differ significantly in their host ranges and ability to cause disease. Unlike eutherian mammals, Australian marsupials (metatherian mammals) have long been thought to be highly susceptible to toxoplasmosis and neosporosis because of their historical isolation from the parasites. In this study, the carnivorous fat-tailed dunnart (Sminthopsis crassicaudata) was used as a disease model to investigate the immune response and susceptibility to infection of an Australian marsupial to T. gondii and N. caninum The disease outcome was more severe in N. caninum-infected dunnarts than in T. gondii-infected dunnarts, as shown by the severity of clinical and histopathological features of disease and higher tissue parasite burdens in the tissues evaluated. Transcriptome sequencing (RNA-seq) of spleens from infected dunnarts and mitogen-stimulated dunnart splenocytes was used to define the cytokine repertoires. Changes in mRNA expression during the time course of infection were measured using quantitative reverse transcription-PCR (qRT-PCR) for key Th1 (gamma interferon [IFN-γ] and tumor necrosis factor alpha [TNF-α]), Th2 (interleukin 4 [IL-4] and IL-6), and Th17 (IL-17A) cytokines. The results show qualitative differences in cytokine responses by the fat-tailed dunnart to infection with N. caninum and T. gondii Dunnarts infected with T. gondii were capable of mounting a more effective Th1 immune response than those infected with N. caninum, indicating the role of the immune response in the outcome scenarios of parasite infection in this marsupial mammal.

Amplification of microsatellite repeat motifs is associated with the evolutionary differentiation and heterochromatinization of sex chromosomes in Sauropsida.

The sex chromosomes in Sauropsida (reptiles and birds) have evolved independently many times. They show astonishing diversity in morphology ranging from cryptic to highly differentiated sex chromosomes with male (XX/XY) and female heterogamety (ZZ/ZW). Comparing such diverse sex chromosome systems thus provides unparalleled opportunities to capture evolution of morphologically differentiated sex chromosomes in action. Here, we describe chromosomal mapping of 18 microsatellite repeat motifs in eight species of Sauropsida. More than two microsatellite repeat motifs were amplified on the sex-specific chromosome, W or Y, in five species (Bassiana duperreyi, Aprasia parapulchella, Notechis scutatus, Chelodina longicollis, and Gallus gallus) of which the sex-specific chromosomes were heteromorphic and heterochromatic. Motifs (AAGG)n and (ATCC)n were amplified on the W chromosome of Pogona vitticeps and the Y chromosome of Emydura macquarii, respectively. By contrast, no motifs were amplified on the W chromosome of Christinus marmoratus, which is not much differentiated from the Z chromosome. Taken together with previously published studies, our results suggest that the amplification of microsatellite repeats is tightly associated with the differentiation and heterochromatinization of sex-specific chromosomes in sauropsids as well as in other taxa. Although some motifs were common between the sex-specific chromosomes of multiple species, no correlation was observed between this commonality and the species phylogeny. Furthermore, comparative analysis of sex chromosome homology and chromosomal distribution of microsatellite repeats between two closely related chelid turtles, C. longicollis and E. macquarii, identified different ancestry and differentiation history. These suggest multiple evolutions of sex chromosomes in the Sauropsida.

Characterization of the antimicrobial peptide family defensins in the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii).

Defensins comprise a family of cysteine-rich antimicrobial peptides with important roles in innate and adaptive immune defense in vertebrates. We characterized alpha and beta defensin genes in three Australian marsupials: the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii) and identified 48, 34, and 39 defensins, respectively. One hundred and twelve have the classical antimicrobial peptides characteristics required for pathogen membrane targeting, including cationic charge (between 1+ and 15+) and a high proportion of hydrophobic residues (>30%). Phylogenetic analysis shows that gene duplication has driven unique and species-specific expansions of devil, koala, and tammar wallaby beta defensins and devil alpha defensins. Defensin genes are arranged in three genomic clusters in marsupials, whereas further duplications and translocations have occurred in eutherians resulting in four and five gene clusters in mice and humans, respectively. Marsupial defensins are generally under purifying selection, particularly residues essential for defensin structural stability. Certain hydrophobic or positively charged sites, predominantly found in the defensin loop, are positively selected, which may have functional significance in defensin-target interaction and membrane insertion.

Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome.

BACKGROUND: Squamates (lizards and snakes) are a speciose lineage of reptiles displaying considerable karyotypic diversity, particularly among lizards. Understanding the evolution of this diversity requires comparison of genome organisation between species. Although the genomes of several squamate species have now been sequenced, only the green anole lizard has any sequence anchored to chromosomes. There is only limited gene mapping data available for five other squamates. This makes it difficult to reconstruct the events that have led to extant squamate karyotypic diversity. The purpose of this study was to anchor the recently sequenced central bearded dragon (Pogona vitticeps) genome to chromosomes to trace the evolution of squamate chromosomes. Assigning sequence to sex chromosomes was of particular interest for identifying candidate sex determining genes. RESULTS: By using two different approaches to map conserved blocks of genes, we were able to anchor approximately 42 % of the dragon genome sequence to chromosomes. We constructed detailed comparative maps between dragon, anole and chicken genomes, and where possible, made broader comparisons across Squamata using cytogenetic mapping information for five other species. We show that squamate macrochromosomes are relatively well conserved between species, supporting findings from previous molecular cytogenetic studies. Macrochromosome diversity between members of the Toxicofera clade has been generated by intrachromosomal, and a small number of interchromosomal, rearrangements. We reconstructed the ancestral squamate macrochromosomes by drawing upon comparative cytogenetic mapping data from seven squamate species and propose the events leading to the arrangements observed in representative species. In addition, we assigned over 8 Mbp of sequence containing 219 genes to the Z chromosome, providing a list of genes to begin testing as candidate sex determining genes. CONCLUSIONS: Anchoring of the dragon genome has provided substantial insight into the evolution of squamate genomes, enabling us to reconstruct ancestral macrochromosome arrangements at key positions in the squamate phylogeny, demonstrating that fusions between macrochromosomes or fusions of macrochromosomes and microchromosomes, have played an important role during the evolution of squamate genomes. Assigning sequence to the sex chromosomes has identified NR5A1 as a promising candidate sex determining gene in the dragon.

Development of a SNP-based assay for measuring genetic diversity in the Tasmanian devil insurance population.

BACKGROUND: The Tasmanian devil (Sarcophilus harrisii) has undergone a recent, drastic population decline due to the highly contagious devil facial tumor disease. The tumor is one of only two naturally occurring transmissible cancers and is almost inevitably fatal. In 2006 a disease-free insurance population was established to ensure that the Tasmanian devil is protected from extinction. The insurance program is dependent upon preserving as much wild genetic diversity as possible to maximize the success of subsequent reintroductions to the wild. Accurate genotypic data is vital to the success of the program to ensure that loss of genetic diversity does not occur in captivity. Until recently, microsatellite markers have been used to study devil population genetics, however as genetic diversity is low in the devil and potentially decreasing in the captive population, a more sensitive genotyping assay is required. METHODS: Utilising the devil reference genome and whole genome re-sequencing data, we have identified polymorphic regions for use in a custom genotyping assay. These regions were amplified using PCR and sequenced on the Illumina MiSeq platform to refine a set a markers to genotype the Tasmanian devil insurance population. RESULTS: We have developed a set of single nucleotide polymorphic (SNP) markers, assayed by amplicon sequencing, that provide a high-throughput method for monitoring genetic diversity and assessing familial relationships among devils. To date we have used a total of 267 unique SNPs within both putatively neutral and functional loci to genotype 305 individuals in the Tasmanian devil insurance population. We have used these data to assess genetic diversity in the population as well as resolve the parentage of 21 offspring. CONCLUSIONS: Our molecular data has been incorporated with studbook management practices to provide more accurate pedigree information and to inform breeding recommendations. The assay will continue to be used to monitor the genetic diversity of the insurance population of Tasmanian devils with the aim of reducing inbreeding and maximizing success of reintroductions to the wild.

Evolutionary history of novel genes on the tammar wallaby Y chromosome: Implications for sex chromosome evolution.

We report here the isolation and sequencing of 10 Y-specific tammar wallaby (Macropus eugenii) BAC clones, revealing five hitherto undescribed tammar wallaby Y genes (in addition to the five genes already described) and several pseudogenes. Some genes on the wallaby Y display testis-specific expression, but most have low widespread expression. All have partners on the tammar X, along with homologs on the human X. Nonsynonymous and synonymous substitution ratios for nine of the tammar XY gene pairs indicate that they are each under purifying selection. All 10 were also identified as being on the Y in Tasmanian devil (Sarcophilus harrisii; a distantly related Australian marsupial); however, seven have been lost from the human Y. Maximum likelihood phylogenetic analyses of the wallaby YX genes, with respective homologs from other vertebrate representatives, revealed that three marsupial Y genes (HCFC1X/Y, MECP2X/Y, and HUWE1X/Y) were members of the ancestral therian pseudoautosomal region (PAR) at the time of the marsupial/eutherian split; three XY pairs (SOX3/SRY, RBMX/Y, and ATRX/Y) were isolated from each other before the marsupial/eutherian split, and the remaining three (RPL10X/Y, PHF6X/Y, and UBA1/UBE1Y) have a more complex evolutionary history. Thus, the small marsupial Y chromosome is surprisingly rich in ancient genes that are retained in at least Australian marsupials and evolved from testis-brain expressed genes on the X.

Comparative transcriptomics reveals striking similarities between the bovine and feline isolates of Tritrichomonas foetus: consequences for in silico drug-target identification.

BACKGROUND: Few, if any, protozoan parasites are reported to exhibit extreme organ tropism like the flagellate Tritrichomonas foetus. In cattle, T. foetus infects the reproductive system causing abortion, whereas the infection in cats results in chronic large bowel diarrhoea. In the absence of a T. foetus genome, we utilized a de novo approach to assemble the transcriptome of the bovine and feline genotype to identify host-specific adaptations and virulence factors specific to each genotype. Furthermore, a subset of orthologs was used to characterize putative druggable targets and expose complications of in silico drug target mining in species with indefinite host-ranges. RESULTS: Illumina RNA-seq reads were assembled into two representative bovine and feline transcriptomes containing 42,363 and 36,559 contigs, respectively. Coding and non-coding regions of the genome libraries revealed striking similarities, with 24,620 shared homolog pairs reduced down to 7,547 coding orthologs between the two genotypes. The transcriptomes were near identical in functional category distribution; with no indication of selective pressure acting on orthologs despite differences in parasite origins/host. Orthologs formed a large proportion of highly expressed transcripts in both genotypes (bovine genotype: 76%, feline genotype: 56%). Mining the libraries for protease virulence factors revealed the cysteine proteases (CP) to be the most common. In total, 483 and 445 bovine and feline T. foetus transcripts were identified as putative proteases based on MEROPS database, with 9 hits to putative protease inhibitors. In bovine T. foetus, CP8 is the preferentially transcribed CP while in the feline genotype, transcription of CP7 showed higher abundance. In silico druggability analysis of the two genotypes revealed that when host sequences are taken into account, drug targets are genotype-specific. CONCLUSION: Gene discovery analysis based on RNA-seq data analysis revealed prominent similarities between the bovine and feline T. foetus, suggesting recent adaptation to their respective host/niche. T. foetus represents a unique case of a mammalian protozoan expanding its parasitic grasp across distantly related host lineages. Consequences of the host-range for in silico drug targeting are exposed here, demonstrating that targets of the parasite in one host are not necessarily ideal for the same parasite in another host.

Molecular evolution of Dmrt1 accompanies change of sex-determining mechanisms in reptilia.

In reptiles, sex-determining mechanisms have evolved repeatedly and reversibly between genotypic and temperature-dependent sex determination. The gene Dmrt1 directs male determination in chicken (and presumably other birds), and regulates sex differentiation in animals as distantly related as fruit flies, nematodes and humans. Here, we show a consistent molecular difference in Dmrt1 between reptiles with genotypic and temperature-dependent sex determination. Among 34 non-avian reptiles, a convergently evolved pair of amino acids encoded by sequence within exon 2 near the DM-binding domain of Dmrt1 distinguishes species with either type of sex determination. We suggest that this amino acid shift accompanied the evolution of genotypic sex determination from an ancestral condition of temperature-dependent sex determination at least three times among reptiles, as evident in turtles, birds and squamates. This novel hypothesis describes the evolution of sex-determining mechanisms as turnover events accompanied by one or two small mutations.

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia.

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

State-transition modeling of blood transcriptome predicts disease evolution and treatment response in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention, before phenotypic changes become detectable.