Differential roles of TGIF family genes in mammalian reproduction

TGIF2 is a potential biomarker for diagnosis and prognosis of glioma

Background

TGFB-induced factor homeobox 2 (TGIF2), a member of the Three-Amino-acid-Loop-Extension (TALE) superfamily, has been implicated in various malignant tumors. However, its prognostic significance in glioma, impact on tumor immune infiltration, and underlying mechanisms in glioma development remain elusive.

Methods

The expression of TGIF2 in various human normal tissues, normal brain tissues, and gliomas was investigated using HPA, TCGA, GTEx, and GEO databases. The study employed several approaches, including Kaplan-Meier analysis, ROC analysis, logistic regression, Cox regression, GO analysis, KEGG analysis, and GSEA, to explore the relationship between TGIF2 expression and clinicopathologic features, prognostic value, and potential biological functions in glioma patients. The impact of TGIF2 on tumor immune infiltration was assessed through Estimate, ssGSEA, and Spearman analysis. Genes coexpressed with TGIF2 were identified, and the protein-protein interaction (PPI) network of these coexpressed genes were constructed using the STRING database and Cytoscape software. Hub genes were identified using CytoHubba plugin, and their clinical predictive value was explored. Furthermore, in vitro experiments were performed by knocking down and knocking out TGIF2 using siRNA and CRISPR/Cas9 gene editing, and the role of TGIF2 in glioma cell invasion and migration was analyzed using transwell assay, scratch wound-healing assay, RT-qPCR, and Western blot.

Results

TGIF2 mRNA was found to be upregulated in 21 cancers, including glioma. High expression of TGIF2 was associated with malignant phenotypes and poor prognosis in glioma patients, indicating its potential as an independent prognostic factor. Furthermore, elevated TGIF2 expression positively correlated with cell cycle regulation, DNA synthesis and repair, extracellular matrix (ECM) components, immune response, and several signaling pathways that promote tumor progression. TGIF2 showed correlations with Th2 cells, macrophages, and various immunoregulatory genes. The hub genes coexpressed with TGIF2 demonstrated significant predictive value. Additionally, in vitro experiments revealed that knockdown and knockout of TGIF2 inhibited glioma cell invasion, migration and suppressed the epithelial-mesenchymal transition (EMT) phenotype.

Conclusion

TGIF2 emerges as a potential biomarker for glioma, possibly linked to tumor immune infiltration and EMT.

Insights into genetic diversity and phenotypic variations in domestic geese through comprehensive population and pan-genome analysis

BACKGROUND

Domestic goose breeds are descended from either the Swan goose (Anser cygnoides) or the Greylag goose (Anser anser), exhibiting variations in body size, reproductive performance, egg production, feather color, and other phenotypic traits. Constructing a pan-genome facilitates a thorough identification of genetic variations, thereby deepening our comprehension of the molecular mechanisms underlying genetic diversity and phenotypic variability.

RESULTS

To comprehensively facilitate population genomic and pan-genomic analyses in geese, we embarked on the task of 659 geese whole genome resequencing data and compiling a database of 155 RNA-seq samples. By constructing the pan-genome for geese, we generated non-reference contigs totaling 612 Mb, unveiling a collection of 2,813 novel genes and pinpointing 15,567 core genes, 1,324 softcore genes, 2,734 shell genes, and 878 cloud genes in goose genomes. Furthermore, we detected an 81.97 Mb genomic region showing signs of genome selection, encompassing the TGFBR2 gene correlated with variations in body weight among geese. Genome-wide association studies utilizing single nucleotide polymorphisms (SNPs) and presence-absence variation revealed significant genomic associations with various goose meat quality, reproductive, and body composition traits. For instance, a gene encoding the SVEP1 protein was linked to carcass oblique length, and a distinct gene-CDS haplotype of the SVEP1 gene exhibited an association with carcass oblique length. Notably, the pan-genome analysis revealed enrichment of variable genes in the "hair follicle maturation" Gene Ontology term, potentially linked to the selection of feather-related traits in geese. A gene presence-absence variation analysis suggested a reduced frequency of genes associated with "regulation of heart contraction" in domesticated geese compared to their wild counterparts. Our study provided novel insights into gene expression features and functions by integrating gene expression patterns across multiple organs and tissues in geese and analyzing population variation.

CONCLUSION

This accomplishment originates from the discernment of a multitude of selection signals and candidate genes associated with a wide array of traits, thereby markedly enhancing our understanding of the processes underlying domestication and breeding in geese. Moreover, assembling the pan-genome for geese has yielded a comprehensive apprehension of the goose genome, establishing it as an indispensable asset poised to offer innovative viewpoints and make substantial contributions to future geese breeding initiatives.

Molecular characterization and effects of the TGIF1 gene on proliferation and steroidogenesis in yak (Bos grunniens) granulosa cells

TG interaction factor 1 (TGIF1) plays a major role in transcriptional inhibition and suppression of TGF-β signaling, but its functional roles in granulosa cells (GCs) have not been elucidated; in particular, there is no information about the yak (Bos grunniens) TGIF1 gene. Therefore, the objectives of this study were to clone yak TGIF1 and investigate TGIF1 functions in yak GCs. RT‒PCR results showed that the coding region of yak TGIF1 is 759 bp and encodes 252 amino acids. Its nucleotide sequence showed 85.24-99.74% similarity to mouse, human, pig, goat and cattle homologous genes. To explore the functional roles of TGIF1, we studied proliferation, apoptosis, cell cycle progression, steroidogenesis and the expression levels of related genes in yak GCs transfected with small interfering RNA specific to TGIF1. The results showed that TGIF1 knockdown promoted proliferation and cell cycle progression and inhibited apoptosis and estradiol (E2) and progesterone (P4) production in cultured yak GCs. Conversely, TGIF1 overexpression inhibited proliferation and cell cycle progression and stimulated apoptosis and E2 and P4 production. In addition, these functional changes in yak GCs were observed parallel to the expression changes in genes involved in the cell cycle (PCNA, CDK2, CCND1, CCNE1, CDK4 and P53), apoptosis (BCL2, BAX and CASPASE3), and steroidogenesis (CYP11A1, 3β-HSD and StAR). In conclusion, TGIF1 was relatively conserved in the course of animal evolution. TGIF1 inhibited GC viability and stimulated apoptosis and the secretion of E2 and P4 by yak GCs. Our results will help to reveal the mechanism underlying yak follicular development and improve the reproductive efficiency of female yaks.

Metabolic Acidosis Results in Sexually Dimorphic Response in the Heart Tissue

Metabolic acidosis (MA) is a highly prevalent disorder in a significant proportion of the population, resulting from imbalance in blood pH homeostasis. The heart, being an organ with very low regenerative capacity and high metabolic activity, is vulnerable to chronic, although low-grade, MA. To systematically characterize the effect of low-grade MA on the heart, we treated male and female mice with NH4Cl supplementation for 2 weeks and analyzed their blood chemistry and transcriptomic signature of the heart tissue. The reduction of pH and plasma bicarbonate levels without an associated change in anion gap indicated a physiological manifestation of low-grade MA with minimal respiratory compensation. On transcriptomic analysis, we observed changes in cardiac-specific genes with significant gender-based differences due to MA. We found many genes contributing to dilated cardiomyopathy to be altered in males, more than in females, while cardiac contractility and Na/K/ATPase-Src signaling were affected in the opposite way. Our model presents a systems-level understanding of how the cardiovascular tissue is affected by MA. As low-grade MA is a common ailment with many dietary and pharmaceutical interventions, our work presents avenues to limit chronic cardiac damage and disease manifestation, as well as highlighting the sex differences in MA-induced cardiovascular damage.

Small extracellular vesicles derived from dermal fibroblasts promote fibroblast activity and skin development through carrying miR-218 and ITGBL1

Skin thickness is closely related to the appearance of human skin, such as sagging and wrinkling, which primarily depends on the level of collagen I synthesized by dermal fibroblasts (DFs). Small extracellular vesicles (SEVs), especially those derived from human DFs (HDFs), are crucial orchestrators in shaping physiological and pathological development of skin. However, the limited supply of human skin prevents the production of a large amount of HDFs-SEVs, and pig skin is used as a model of human skin. In this study, SEVs derived from DFs of Chenghua pigs (CH-SEVs), considered to have superior skin thickness, and Large White pigs (LW-SEVs) were collected to compare their effects on DFs and skin tissue. Our results showed that, compared with LW-SEVs, CH-SEVs more effectively promoted fibroblast proliferation, migration, collagen synthesis and contraction; in addition, in mouse model injected with both SEVs, compared with LW-SEVs, CH-SEVs increased the skin thickness and collagen I content more effectively. Some differentially expressed miRNAs and proteins were found between CH-SEVs and LW-SEVs by small RNA-seq and LC-MS/MS analysis. Interestingly, we identified that CH-SEVs were enriched in miRNA-218 and ITGBL1 protein, which played important roles in promoting fibroblast activity via activation of the downstream TGFβ1-SMAD2/3 pathway in vitro. Furthermore, overexpression of miRNA-218 and ITGBL1 protein increased the thickness and collagen I content of mouse skin in vivo. These results indicate that CH-SEVs can effectively stimulate fibroblast activity and promote skin development and thus have the potential to protect against and repair skin damage.

Species persistence with hybridization in toad-headed lizards driven by divergent selection and low recombination

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

HIV-1 Vpr protein upregulates microRNA-210-5p expression to induce G2 arrest by targeting TGIF2

MicroRNAs (miRNAs) are important molecules that mediate virus-host interactions, mainly by regulating gene expression via gene silencing. Here, we demonstrated that HIV-1 infection upregulated miR-210-5p in HIV-1-inoculated cell lines and in the serum of HIV-1-infected individuals. Luciferase reporter assays and western blotting confirmed that a target protein of miR-210-5p, TGIF2, is regulated by HIV-1 infection. Furthermore, HIV-1 Vpr protein induced miR-210-5p expression. The use of a miR-210-5p inhibitor and TGIF2 overexpression showed that Vpr upregulated miR-210-5p and thereby downregulated TGIF2, which might be one of the mechanisms used by Vpr to induce G2 arrest. Moreover, we identified a transcription factor, NF-κB p50, which upregulated miR-210-5p in response to Vpr protein. In conclusion, we identified a mechanism whereby miR-210-5p, which is induced upon HIV-1 infection, targets TGIF2. This pathway was initiated by Vpr protein activating NF-κB p50, which promoted G2 arrest. These alterations orchestrated by miRNA provide new evidence on how HIV-1 interacts with its host during infection and increase our understanding of the mechanism by which Vpr regulates the cell cycle.

The homeobox gene TGIF1 is required for chicken ovarian cortical development and generation of the juxtacortical medulla

During early embryogenesis in amniotic vertebrates, the gonads differentiate into either ovaries or testes. The first cell lineage to differentiate gives rise to the supporting cells: Sertoli cells in males and pre-granulosa cells in females. These key cell types direct the differentiation of the other cell types in the gonad, including steroidogenic cells. The gonadal surface epithelium and the interstitial cell populations are less well studied, and little is known about their sexual differentiation programs. Here, we show the requirement of the homeobox transcription factor gene TGIF1 for ovarian development in the chicken embryo. TGIF1 is expressed in the two principal ovarian somatic cell populations: the cortex and the pre-granulosa cells of the medulla. TGIF1 expression is associated with an ovarian phenotype in estrogen-mediated sex reversal experiments. Targeted misexpression and gene knockdown indicate that TGIF1 is required, but not sufficient, for proper ovarian cortex formation. In addition, TGIF1 is identified as the first known regulator of juxtacortical medulla development. These findings provide new insights into chicken ovarian differentiation and development, specifically cortical and juxtacortical medulla formation.

Effects of TG interaction factor 1 on synthesis of estradiol and progesterone in granulosa cells of goats through SMAD2/3-SP1 signaling pathway

The TG interaction factor 1 (TGIF1) is of the TALE homologue domain protein family and is considered as a transcriptional repressor of SMAD protein that interacts with DNA through a specific consensus-binding site for TG and recruits mSin3A and histone deacetylases to the SMAD complex. In this study, there is the first detailed description of TGIF1 on steroidogenesis in goat granulosa cells. When there is a relatively greater expression of the TGIF1 gene, there is a lesser abundance of CYP11A1, CYP19A1, and StAR mRNA transcript and protein and 3β-HSD mRNA transcript in granulosa cells of goats. Furthermore, there were lesser concentrations of 17β-estradiol (E₂) and progesterone (P₄) in culture medium when there was greater TGIF1 gene expression and there were greater concentrations of these hormones in the culture medium when there was lesser TGIF1 gene expression. There may be functions of TGIF1, therefore, in suppression of SMAD-induced E₂ and P₄ production and in decreasing the phosphorylation of SMAD2/3 in granulosa cells of goats and relative abundance of the SMAD2/3 protein transcription factor, SP1. With suppression of TGIF1 gene expression, there was a reversal of SP1-induced suppression of steroidogenesis-related genes. Results of the present study provide insights about the potential mechanism underlying the regulation of granulosa cell steroidogenesis of goats by TGIF1.

Genome-wide occupancy reveals the localization of H1T2 (H1fnt) to repeat regions and a subset of transcriptionally active chromatin domains in rat spermatids

Background

H1T2/H1FNT is a germ cell-specific linker histone variant expressed during spermiogenesis specifically in round and elongating spermatids. Infertile phenotype of homozygous H1T2 mutant male mice revealed the essential function of H1T2 for the DNA condensation and histone-to-protamine replacement in spermiogenesis. However, the mechanism by which H1T2 imparts the inherent polarity within spermatid nucleus including the additional protein partners and the genomic domains occupied by this linker histone are unknown.

Results

Sequence analysis revealed the presence of Walker motif, SR domains and putative coiled-coil domains in the C-terminal domain of rat H1T2 protein. Genome-wide occupancy analysis using highly specific antibody against the CTD of H1T2 demonstrated the binding of H1T2 to the LINE L1 repeat elements and to a significant percentage of the genic regions (promoter-TSS, exons and introns) of the rat spermatid genome. Immunoprecipitation followed by mass spectrometry analysis revealed the open chromatin architecture of H1T2 occupied chromatin encompassing the H4 acetylation and other histone PTMs characteristic of transcriptionally active chromatin. In addition, the present study has identified the interacting protein partners of H1T2-associated chromatin mainly as nucleo-skeleton components, RNA-binding proteins and chaperones.

Conclusions

Linker histone H1T2 possesses unique domain architecture which can account for the specific functions associated with chromatin remodeling events facilitating the initiation of histone to transition proteins/protamine transition in the polar apical spermatid genome. Our results directly establish the unique function of H1T2 in nuclear shaping associated with spermiogenesis by mediating the interaction between chromatin and nucleo-skeleton, positioning the epigenetically specialized chromatin domains involved in transcription coupled histone replacement initiation towards the apical pole of round/elongating spermatids.

TGIF1 and SF1 polymorphisms are associated with litter size in Small Tail Han sheep

TGF‐β induced factor homeobox 1 (TGIF1) and splicing factor 1 (SF1) are important for mammalian reproduction; however, the effects of these genes on litter size in sheep remain unexplored. In this study, we genotyped 768 ewes from seven sheep breeds at two loci: g.37871539C>T, a synonymous mutation of TGIF1; and g.42314637T>C, a 3′UTR variant of SF1. Our analysis of polymorphism revealed only two genotypes at locus g.37871539C>T in TGIF1, with most sheep populations being moderately polymorphic (0.25 < PIC < 0.5) at this site. In contrast, most breeds exhibited low polymorphism (PIC ≤0.25) at the SF1 locus g.42314637T>C. The association analysis revealed that a synonymous mutation at g.37871539C>T in TGIF1 was highly associated with litter size in Small Tail Han sheep, in which it causes a significant decrease in litter size. Conversely, while the SF1 3′UTR variant g.42314637T>C was also highly associated with litter size in sheep, it causes a significant increase in the number of litter size. Combined, these data provide valuable information regarding candidate genetic markers for sheep breeding programs.

Pituitary Transcription Factor Mutations Leading to Hypopituitarism

Congenital pituitary hormone deficiency is a disabling condition. It is part of a spectrum of disorders including craniofacial midline developmental defects ranging from holoprosencephaly through septo-optic dysplasia to combined and isolated pituitary hormone deficiency. The first genes discovered in the human disease were based on mouse models of dwarfism due to mutations in transcription factor genes. High-throughput DNA sequencing technologies enabled clinicians and researchers to find novel genetic causes of hypopituitarism for the more than three quarters of patients without a known genetic diagnosis to date. Transcription factor (TF) genes are at the forefront of the functional analysis of novel variants of unknown significance due to the relative ease in in vitro testing in a research lab. Genetic testing in hypopituitarism is of high importance to the individual and their family to predict phenotype composition, disease progression and to avoid life-threatening complications such as secondary adrenal insufficiency.This chapter aims to highlight our current understanding about (1) the contribution of TF genes to pituitary development (2) the diversity of inheritance and phenotype features in combined and select isolated pituitary hormone deficiency and (3) provide an initial assessment on how to approach variants of unknown significance in human hypopituitarism. Our better understanding on how transcription factor gene variants lead to hypopituitarism is a meaningful step to plan advanced therapies to specific genetic changes in the future.

Silencing of TGIF sensitizes MDA-MB-231 human breast cancer cells to cisplatin-induced apoptosis

The present study was designed to explore the sensitivity of MDA-MB-231 cells to cisplatin after silencing the expression of TG-interacting factor (TGIF) protein. Cell viability was measured using an MTT assay. Cell apoptosis was detected by the annexin V and dead cell assay and the Hoechst staining assay. Protein expression was analyzed using western blot analysis. A colony formation assay was also performed. It was observed that cisplatin reduced the expression of TGIF protein in a dose- and time-dependent manner. Silencing TGIF significantly suppressed the cell proliferation and colony formation in MDA-MB-231 cells with the treatment of cisplatin. Results indicated that silencing TGIF could dramatically increase the cisplatin-induced apoptosis rate in MDA-MB-231 cells. The expression of PARP and caspase-3 proteins was correlated with the effect that silencing TGIF enhanced cisplatin sensitivity in MDA-MB-231 cells. The present data showed that silencing TGIF promoted apoptotic sensitivity that was induced by cisplatin in MDA-MB-231 human breast cancer cells and suggested that TGIF might be a therapeutic target for improving the chemotherapy response in triple-negative breast cancer.

Beyond cancer genes: colorectal cancer as robust intrinsic states formed by molecular interactions

Colorectal cancer (CRC) has complex pathological features that defy the linear-additive reasoning prevailing in current biomedicine studies. In pursuing a mechanistic understanding behind such complexity, we constructed a core molecular–cellular interaction network underlying CRC and investigated its nonlinear dynamical properties. The hypothesis and modelling method has been developed previously and tested in various cancer studies. The network dynamics reveal a landscape of several attractive basins corresponding to both normal intestinal phenotype and robust tumour subtypes, identified by their different molecular signatures. Comparison between the modelling results and gene expression profiles from patients collected at the second affiliated hospital of Zhejiang University is presented as validation. The numerical ‘driving’ experiment suggests that CRC pathogenesis may depend on pathways involved in gastrointestinal track development and molecules associated with mesenchymal lineage differentiation, such as Stat5, BMP, retinoic acid signalling pathways, Runx and Hox transcription families. We show that the multi-faceted response to immune stimulation and therapies, as well as different carcinogenesis and metastasis routes, can be straightforwardly understood and analysed under such a framework.

Y-chromosomal genes affecting male fertility: A review

The mammalian sex-chromosomes (X and Y) have evolved from autosomes and are involved in sex determination and reproductive traits. The Y-chromosome is the smallest chromosome that consists of 2-3% of the haploid genome and may contain between 70 and 200 genes. The Y-chromosome plays major role in male fertility and is suitable to study the evolutionary relics, speciation, and male infertility and/or subfertility due to its unique features such as long non-recombining region, abundance of repetitive sequences, and holandric inheritance pattern. During evolution, many holandric genes were deleted. The current review discusses the mammalian holandric genes and their functions. The commonly encountered infertility and/or subfertility problems due to point or gross mutation (deletion) of the Y-chromosomal genes have also been discussed. For example, loss or microdeletion of sex-determining region, Y-linked gene results in XY males that exhibit female characteristics, deletion of RNA binding motif, Y-encoded in azoospermic factor b region results in the arrest of spermatogenesis at meiosis. The holandric genes have been covered for associating the mutations with male factor infertility.

The X factor: X chromosome dosage compensation in the evolutionarily divergent monotremes and marsupials

Marsupials and monotremes represent evolutionarily divergent lineages from the majority of extant mammals which are eutherian, or placental, mammals. Monotremes possess multiple X and Y chromosomes that appear to have arisen independently of eutherian and marsupial sex chromosomes. Dosage compensation of X-linked genes occurs in monotremes on a gene-by-gene basis, rather than through chromosome-wide silencing, as is the case in eutherians and marsupials. Specifically, studies in the platypus have shown that for any given X-linked gene, a specific proportion of nuclei within a cell population will silence one locus, with the percentage of cells undergoing inactivation at that locus being highly gene-specific. Hence, it is perhaps not surprising that the expression level of X-linked genes in female platypus is almost double that in males. This is in contrast to the situation in marsupials where one of the two X chromosomes is inactivated in females by the long non-coding RNA RSX, a functional analogue of the eutherian XIST. However, marsupial X chromosome inactivation differs from that seen in eutherians in that it is exclusively the paternal X chromosome that is silenced. In addition, marsupials appear to have globally upregulated X-linked gene expression in both sexes, thus balancing their expression levels with those of the autosomes, a process initially proposed by Ohno in 1967 as being a fundamental component of the X chromosome dosage compensation mechanism but which may not have evolved in eutherians.

Induction of Apoptosis and Growth Suppression by Homeobox Gene TGIFLX in Prostate Cancer Cell Line Lncap

BACKGROUND

TGIFLX, a Homoproteins cluster member located on the X chromosome, has a critical role in male reproduction and prostate development. Previous studies have shown the erratic expression of TGIFLX gene in a large proportion of prostate tumors. However TGIFLX function in prostate development remains unknown. The purpose of this study was to evaluate the consequences of TGIFLX expression on prostate cancer cell lines (LNCaP).

METHOD

Inducible Tet-On gene expression system was used with a regulatory capability by doxycycline induction. In this system, stable LNCaP cells with TGIFLX tet-on plasmid were able to induce TGIFLX expression by doxycycline treatment. TGIFLX gene expression was confirmed by RT-PCR.

RESULTS

Induction of gene expression caused cell proliferation decrement and apoptosis increment in LNCaP TGIFLX cells compared with control cells (P<0.01). Also, by using PEGFPN1 plasmid protein in this study localization was shown in nucleus. The gene was cloned in the plasmid and transfected to LNcap cells with plasmid PEGFPN1 TGIFLX and the plasmid was PEGFPN1. The TGIFLX expression was confirmed by RT-PCR and fluorescent microscopy.

CONCLUSION

TGIFLX expression demonstrated a tumor suppressor characterization in a prostatic cancer cell line with low grade of tumorigenicity (LNCaP). More cell lines with different level of tumorogenicity need to be investigated for further clarification of the TGIFLX gene function.

Genome sequence of an Australian kangaroo, Macropus eugenii, provides insight into the evolution of mammalian reproduction and development

BACKGROUND

We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development.

RESULTS

The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements.

CONCLUSIONS

Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.