Genome analysis of the platypus reveals unique signatures of evolution

Marsupial cathelicidins: characterization, antimicrobial activity and evolution in this unique mammalian lineage

Introduction

Cathelicidins are a family of antimicrobial peptides well-known for their antimicrobial and immunomodulatory functions in eutherian mammals such as humans. However, cathelicidins in marsupials, the other major lineage of mammals, have received little attention despite lineage-specific gene expansions resulting in a large and diverse peptide repertoire.

Methods

We characterized cathelicidins across the marsupial family tree and investigated genomic organisation and evolutionary relationships amongst mammals. Ancestral sequence reconstruction was used to predict ancestral marsupial cathelicidins, which, alongside extant peptides, were synthesized and screened for antimicrobial activity.

Results

We identified 130 cathelicidin genes amongst 14 marsupial species representing 10 families, with gene expansions identified in all species. Cathelicidin genes were encoded in a highly syntenic region of the genome amongst all mammals, although the number of gene clusters differed amongst lineages (eutherians one, marsupials two, and monotremes three). 32 extant and ancestral marsupial cathelicidins displayed rapid, potent, and/or broad-spectrum antibacterial and antifungal activity. Phylogenetic analysis revealed that marsupial and monotreme cathelicidin repertoires may reflect both mammals and birds, as they encode non-classical cathelicidins found only in birds, as well as multiple copies of neutrophil granule protein and classic cathelicidins found only in eutherian mammals.

Conclusion

This study sheds light on the evolutionary history of mammalian cathelicidins and highlights the potential of wildlife for novel bioactive peptide discovery.

Mammalian Species-Specific Resistance to Mammary Cancer

Tumorigenesis in mammals is driven by inherited genetic variants, environmental factors and random errors during normal DNA replication that lead to cancer-causing mutations. These factors initiate uncontrolled cellular proliferation and disrupt the regulation of critical checkpoints. A few mammalian species possess unique protective mechanisms that enable them to resist widespread cancer development and achieve longevity. Tissue-specific tumor protection adds another layer of complexity to this diversity. Breast cancer is a leading cause of human mortality, particularly among females. Driven by the need for new strategies in treatment and prevention, this opinion article explores and supports the idea that herbivores are more resistant to mammary cancer than carnivores and omnivores. This diversity has occurred despite the remarkably similar basic mammary biology. Herbivores' meatless diet cannot explain the differences in cancer resistance, which have accompanied species segregation since the Jurassic era. To investigate the causes of this diversity, the characteristics of tumorigenesis in the human breast-and to a lesser extent in other carnivores-have been compared with data from retrospective analyses of bovine mammary tumor development across various locations over the past century. Well-established genomic, cellular, and systemic triggers of breast cancer exhibit different, or less pronounced tissue-specific activity in the bovine mammary gland, accompanied by novel bovine-specific protective mechanisms. Together, these factors contribute to the near absence of breast cancer in bovines and offer a basis for developing future anticancer strategies.

The Spleen Virome of Australia's Endemic Platypus Is Dominated by Highly Diverse Papillomaviruses

The platypus (Ornithorhynchus anatinus), a unique monotreme, represents a pivotal point in mammalian evolution with its distinctive traits, such as electroreception and venom production. Despite its evolutionary significance, the viral diversity within platypuses remains poorly understood. This study employed next-generation sequencing to investigate the virome of the dead platypuses, uncovering a range of novel and divergent viruses. Among the identified viruses were four complete genomes of papillomaviruses (OaPV1-4) exhibiting substantial divergence from known strains, suggesting a novel genus within the subfamily Secondpapillomavirinae. Additionally, five novel parvoviruses were detected, including two with complete genomes, highlighting the complex viral ecosystem of the platypus. Phylogenetic analysis placed these viruses in unique evolutionary branches, further demonstrating the platypus's evolutionary significance. A circular DNA virus, a tombus-like virus, and a nodamuvirus were also identified, expanding the understanding of viral diversity in monotremes. These findings offer crucial insights into viral evolution in one of the most unique mammalian lineages, emphasising the need for further exploration to assess ecological and pathological impacts on platypus populations.

Chromosome-level echidna genome illuminates evolution of multiple sex chromosome system in monotremes

BACKGROUND

A thorough analysis of genome evolution is fundamental for biodiversity understanding. The iconic monotremes (platypus and echidna) feature extraordinary biology. However, they also exhibit rearrangements in several chromosomes, especially in the sex chromosome chain. Therefore, the lack of a chromosome-level echidna genome has limited insights into genome evolution in monotremes, in particular the multiple sex chromosomes complex.

RESULTS

Here, we present a new long reads-based chromosome-level short-beaked echidna (Tachyglossus aculeatus) genome, which allowed the inference of chromosomal rearrangements in the monotreme ancestor (2n = 64) and each extant species. Analysis of the more complete sex chromosomes uncovered homology between 1 Y chromosome and multiple X chromosomes, suggesting that it is the ancestral X that has undergone reciprocal translocation with ancestral autosomes to form the complex. We also identified dozens of ampliconic genes on the sex chromosomes, with several ancestral ones expressed during male meiosis, suggesting selective constraints in pairing the multiple sex chromosomes.

CONCLUSION

The new echidna genome provides an important basis for further study of the unique biology and conservation of this species.

The evolution of flexibility and function in the Fc domains of IgM, IgY, and IgE

Introduction

Antibody Fc regions harbour the binding sites for receptors that mediate effector functions following antigen engagement by the Fab regions. An extended "hinge" region in IgG allows flexibility between Fab and Fc, but in both the most primitive antibody, IgM, and in the evolutionarily more recent IgE, the hinge is replaced by an additional domain pair in the homodimeric six-domain Fc region. This permits additional flexibility within the Fc region, which has been exploited by nature to modulate antibody effector functions. Thus, in pentameric or hexameric IgM, the Fc regions appear to adopt a planar conformation in solution until antigen binding causes a conformational change and exposes the complement binding sites. In contrast, IgE-Fc principally adopts an acutely bent conformation in solution, but the binding of different receptors is controlled by the degree of bending, and there is allosteric communication between receptor binding sites.

Methods

We sought to trace the evolution of Fc conformational diversity from IgM to IgE via the intermediate avian IgY by studying the solution conformations of their Fc regions by small-angle X-ray scattering. We compared four extant proteins: human IgM-Fc homodimer, chicken IgY-Fc, platypus IgE-Fc, and human IgE-Fc. These are examples of proteins that first appeared in the jawed fish [425 million years ago (mya)], tetrapod (310 mya), monotreme (166 mya), and hominid (2.5 mya) clades, respectively.

Results and discussion

We analysed the scattering curves in terms of contributions from a pool of variously bent models chosen by a non-negative linear least-squares algorithm and found that the four proteins form a series in which the proportion of acutely bent material increases: IgM-Fc < IgY-Fc < plIgE-Fc < huIgE-Fc. This follows their order of appearance in evolution. For the huIgM-Fc homodimer, although none are acutely bent, and a significant fraction of the protein is sufficiently bent to expose the C1q-binding site, it predominantly adopts a fully extended conformation. In contrast, huIgE-Fc is found principally to be acutely bent, as expected from earlier studies. IgY-Fc, in this first structural analysis of the complete Fc region, exhibits an ensemble of conformations from acutely bent to fully extended, reflecting IgY's position as an evolutionary intermediate between IgM and IgE.

Compensation of gene dosage on the mammalian X

Changes in gene dosage can have tremendous evolutionary potential (e.g. whole-genome duplications), but without compensatory mechanisms, they can also lead to gene dysregulation and pathologies. Sex chromosomes are a paradigmatic example of naturally occurring gene dosage differences and their compensation. In species with chromosome-based sex determination, individuals within the same population necessarily show 'natural' differences in gene dosage for the sex chromosomes. In this Review, we focus on the mammalian X chromosome and discuss recent new insights into the dosage-compensation mechanisms that evolved along with the emergence of sex chromosomes, namely X-inactivation and X-upregulation. We also discuss the evolution of the genetic loci and molecular players involved, as well as the regulatory diversity and potentially different requirements for dosage compensation across mammalian species.

Body fat and circulating leptin levels in the captive short-beaked echidna (Tachyglossus aculeatus)

It is possible that the reproductive strategy of the short-beaked echidna is related to seasonal changes in fat deposition and energy availability, regulated by seasonal changes in endocrine function. We predicted that circulating leptin levels would be directly proportional to adiposity during most of the year, but that a change in this relationship would occur during the pre-breeding season to allow increased fat deposition. To test this hypothesis, we made use of a captive colony of echidnas to describe and quantify changes in fat distribution and the adipostatic hormone leptin. First we assessed seasonal changes in circulating leptin levels, body mass and adiposity for three male and three female adult echidnas maintained on a standard diet. Second, we explored the relationship between circulating leptin levels and increased caloric intake for an additional five adult female echidnas that were provided with supplemented nutrition. Third we visualised fat distribution in male and female adult echidnas using magnetic resonance imaging (MRI) before and after the breeding season, to determine where fat is deposited in this species. For echidnas maintained on the standard diet, there were no seasonal changes in body mass, body fat or plasma leptin levels. However, female echidnas provided with supplemented nutrition had significantly elevated plasma leptin levels during the breeding season, compared to the pre-and post- breeding periods. MRI showed substantial subcutaneous fat depots extending dorso-laterally from the base of the skull to the base of the tail, in both sexes. Pre-breeding season, both sexes had considerable fat deposition in the pelvic/rump region, whilst the female echidna accumulated most fat in the abdominal region. This study shows that male and female echidnas accumulate body fat in the pelvic/rump and the abdominal regions, respectively and that circulating leptin may promote fattening in female echidnas during the breeding season by means of leptin resistance. However, further research is required to evaluate the precise relationship between seasonal changes in leptin and adiposity.

A family of olfactory receptors uniquely expanded in marsupial and monotreme genomes are expressed by a T cell subset also unique to marsupials and monotremes

Olfactory receptors (OR), expressed on olfactory neurons, mediate the sense of smell. Recently, OR have also been shown to be expressed in non-olfactory tissues, including cells of the immune system. An analysis of single-cell transcriptomes of splenocytes of the grey short-tailed opossum (Monodelphis domestica) found OR are expressed on a subset of T cells, the γμ T cells, that are unique to marsupials and monotremes. A majority of opossum γμ T cells transcriptomes contain OR family 14 transcripts, specifically, from the OR14C subfamily. Amongst the mammals, the OR14 gene family is expanded in the genomes of marsupials and monotremes, and rarer or absent in placental mammals. In summary, here we demonstrate the intriguing correlation that a family of OR genes, abundant in the genomes of marsupials and monotremes, are ectopically expressed in a particular subset of T cells unique to the marsupials and monotremes.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Comparative genomics of monotremes provides insights into the early evolution of mammalian epidermal differentiation genes

The function of the skin as a barrier against the environment depends on the differentiation of epidermal keratinocytes into highly resilient corneocytes that form the outermost skin layer. Many genes encoding structural components of corneocytes are clustered in the epidermal differentiation complex (EDC), which has been described in placental and marsupial mammals as well as non-mammalian tetrapods. Here, we analyzed the genomes of the platypus (Ornithorhynchus anatinus) and the echidna (Tachyglossus aculeatus) to determine the gene composition of the EDC in the basal clade of mammals, the monotremes. We report that mammal-specific subfamilies of EDC genes encoding small proline-rich proteins (SPRRs) and late cornified envelope proteins as well as single-copy EDC genes such as involucrin are conserved in monotremes, suggesting that they have originated in stem mammals. Monotremes have at least one gene homologous to the group of filaggrin (FLG), FLG2 and hornerin (HRNR) in placental mammals, but no clear one-to-one pairwise ortholog of either FLG, FLG2 or HRNR. Caspase-14, a keratinocyte differentiation-associated protease implicated in the processing of filaggrin, is encoded by at least 3 gene copies in the echidna. Our results reveal evolutionarily conserved and clade-specific features of the genetic regulation of epidermal differentiation in monotremes.

Acquisition and Exaptation of Endogenous Retroviruses in Mammalian Placenta

Endogenous retroviruses (ERVs) are retrovirus-like sequences that were previously integrated into the host genome. Although most ERVs are inactivated by mutations, deletions, or epigenetic regulation, some remain transcriptionally active and impact host physiology. Several ERV-encoded proteins, such as Syncytins and Suppressyn, contribute to placenta acquisition, a crucial adaptation in mammals that protects the fetus from external threats and other risks while enabling the maternal supply of oxygen, nutrients, and antibodies. In primates, Syncytin-1 and Syncytin-2 facilitate cell-cell fusion for placental formation. Suppressyn is the first ERV-derived protein that inhibits cell fusion by binding to ASCT2, the receptor for Syncytin-1. Furthermore, Syncytin-2 likely inserted into the genome of the common ancestor of Anthropoidea, whereas Syncytin-1 and Suppressyn likely inserted into the ancestor of catarrhines; however, they were inactivated in some lineages, suggesting that multiple exaptation events had occurred. This review discusses the role of ERV-encoded proteins, particularly Syncytins and Suppressyn, in placental development and function, focusing on the integration of ERVs into the host genome and their contribution to the genetic mechanisms underlying placentogenesis. This review provides valuable insights into the molecular and genetic aspects of placentation, potentially shedding light on broader evolutionary and physiological processes in mammals.

Mammalian neurotoxins, Blarina paralytic peptides, cause hyperpolarization of human T-type Ca channel hCav3.2 activation

Among the rare venomous mammals, the short-tailed shrew Blarina brevicauda has been suggested to produce potent neurotoxins in its saliva to effectively capture prey. Several kallikrein-like lethal proteases have been identified, but the active substances of B. brevicauda remained unclear. Here, we report Blarina paralytic peptides (BPPs) 1 and 2 isolated from its submaxillary glands. Synthetic BPP2 showed mealworm paralysis and a hyperpolarization shift (-11 mV) of a human T-type Ca2+ channel (hCav3.2) activation. The amino acid sequences of BPPs were similar to those of synenkephalins, which are precursors of brain opioid peptide hormones that are highly conserved among mammals. However, BPPs rather resembled centipede neurotoxic peptides SLPTXs in terms of disulfide bond connectivity and stereostructure. Our results suggested that the neurotoxin BPPs were the result of convergent evolution as homologs of nontoxic endogenous peptides that are widely conserved in mammals. This finding is of great interest from the viewpoint of the chemical evolution of vertebrate venoms.

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

Return of the forgotten hero: the role of Y chromosome-encoded Zfy in male reproduction

The Y-linked zinc finger gene ZFY is conserved across eutherians and is known to be a critical fertility factor in some species. The initial studies of the mouse homologues, Zfy1 and Zfy2, were performed using mice with spontaneous Y chromosome mutations and Zfy transgenes. These studies revealed that Zfy is involved in multiple processes during spermatogenesis, including removal of germ cells with unpaired chromosomes and control of meiotic sex chromosome inactivation during meiosis I, facilitating the progress of meiosis II, promoting spermiogenesis, and improving assisted reproduction outcomes. Zfy was also identified as a key gene in Y chromosome evolution, protecting this chromosome from extinction by serving as the executioner responsible for meiosis surveillance. Studies with targeted Zfy knock-outs revealed that mice lacking both homologues have severe spermatogenic defects and are infertile. Based on protein structure and in vitro assays, Zfy is expected to drive spermatogenesis as a transcriptional regulator. The combined evidence documents that the presence of at least one Zfy homologue is required for male fertility and that Zfy2 plays a more prominent role. This knowledge reinforces the importance of these factors for mouse spermatogenesis and informs our understanding of the human ZFY variants, which are homologous to the mouse Zfy1 and Zfy2.

The Fast and the Furriest: Investigating the Rate of Selection on Mammalian Toxins

The evolution of venom and the selection pressures that act on toxins have been increasingly researched within toxinology in the last two decades, in part due to the exceptionally high rates of diversifying selection observed in animal toxins. In 2015, Sungar and Moran proposed the 'two-speed' model of toxin evolution linking evolutionary age of a group to the rates of selection acting on toxins but due to a lack of data, mammals were not included as less than 30 species of venomous mammal have been recorded, represented by elusive species which produce small amounts of venom. Due to advances in genomics and transcriptomics, the availability of toxin sequences from venomous mammals has been increasing. Using branch- and site-specific selection models, we present the rates of both episodic and pervasive selection acting upon venomous mammal toxins as a group for the first time. We identified seven toxin groups present within venomous mammals, representing Chiroptera, Eulipotyphla and Monotremata: KLK1, Plasminogen Activator, Desmallipins, PACAP, CRiSP, Kunitz Domain One and Kunitz Domain Two. All but one group (KLK1) was identified by our results to be evolving under both episodic and pervasive diversifying selection with four toxin groups having sites that were implicated in the fitness of the animal by TreeSAAP (Selection on Amino Acid Properties). Our results suggest that venomous mammal ecology, behaviour or genomic evolution are the main drivers of selection, although evolutionary age may still be a factor. Our conclusion from these results indicates that mammalian toxins are following the two-speed model of selection, evolving predominately under diversifying selection, fitting in with other younger venomous taxa like snakes and cone snails-with high amounts of accumulating mutations, leading to more novel adaptions in their toxins.

False gene and chromosome losses in genome assemblies caused by GC content variation and repeats

BACKGROUND

Many short-read genome assemblies have been found to be incomplete and contain mis-assemblies. The Vertebrate Genomes Project has been producing new reference genome assemblies with an emphasis on being as complete and error-free as possible, which requires utilizing long reads, long-range scaffolding data, new assembly algorithms, and manual curation. A more thorough evaluation of the recent references relative to prior assemblies can provide a detailed overview of the types and magnitude of improvements.

RESULTS

Here we evaluate new vertebrate genome references relative to the previous assemblies for the same species and, in two cases, the same individuals, including a mammal (platypus), two birds (zebra finch, Anna's hummingbird), and a fish (climbing perch). We find that up to 11% of genomic sequence is entirely missing in the previous assemblies. In the Vertebrate Genomes Project zebra finch assembly, we identify eight new GC- and repeat-rich micro-chromosomes with high gene density. The impact of missing sequences is biased towards GC-rich 5'-proximal promoters and 5' exon regions of protein-coding genes and long non-coding RNAs. Between 26 and 60% of genes include structural or sequence errors that could lead to misunderstanding of their function when using the previous genome assemblies.

CONCLUSIONS

Our findings reveal novel regulatory landscapes and protein coding sequences that have been greatly underestimated in previous assemblies and are now present in the Vertebrate Genomes Project reference genomes.

Widespread false gene gains caused by duplication errors in genome assemblies

Background

False duplications in genome assemblies lead to false biological conclusions. We quantified false duplications in popularly used previous genome assemblies for platypus, zebra finch, and Anna’s Hummingbird, and their new counterparts of the same species generated by the Vertebrate Genomes Project, of which the Vertebrate Genomes Project pipeline attempted to eliminate false duplications through haplotype phasing and purging. These assemblies are among the first generated by the Vertebrate Genomes Project where there was a prior chromosomal level reference assembly to compare with.

Results

Whole genome alignments revealed that 4 to 16% of the sequences are falsely duplicated in the previous assemblies, impacting hundreds to thousands of genes. These lead to overestimated gene family expansions. The main source of the false duplications is heterotype duplications, where the haplotype sequences were relatively more divergent than other parts of the genome leading the assembly algorithms to classify them as separate genes or genomic regions. A minor source is sequencing errors. Ancient ATP nucleotide binding gene families have a higher prevalence of false duplications compared to other gene families. Although present in a smaller proportion, we observe false duplications remaining in the Vertebrate Genomes Project assemblies that can be identified and purged.

Conclusions

This study highlights the need for more advanced assembly methods that better separate haplotypes and sequence errors, and the need for cautious analyses on gene gains.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02764-1.

Megasatellite formation and evolution in vertebrate genes

Since formation of the first proto-eukaryotes, gene repertoire and genome complexity have significantly increased. Among genetic elements responsible for this increase are tandem repeats. Here we describe a genome-wide analysis of large tandem repeats, called megasatellites, in 58 vertebrate genomes. Two bursts occurred, one after the radiation between Agnatha and Gnathostomata fishes and the second one in therian mammals. Megasatellites are enriched in subtelomeric regions and frequently encoded in genes involved in transcription regulation, intracellular trafficking, and cell membrane metabolism, reminiscent of what is observed in fungus genomes. The presence of many introns within young megasatellites suggests that an exon-intron DNA segment is first duplicated and amplified before accumulation of mutations in intronic parts partially erases the megasatellite in such a way that it becomes detectable only in exons. Our results suggest that megasatellite formation and evolution is a dynamic and still ongoing process in vertebrate genomes.

Monotreme-Specific Conserved Putative Proteins Derived from Retroviral Reverse Transcriptase

Endogenous retroviruses (ERVs) have played an essential role in the evolution of mammals. ERV-derived genes are reported in the therians, many of which are involved in placental development; however, the contribution of the ERV-derived genes in monotremes, which are oviparous mammals, remains to be uncovered. Here, we conducted a comprehensive search for possible ERV-derived genes in platypus and echidna genomes and identified three reverse transcriptase-like genes named RTOM1, RTOM2, and RTOM3 clustered in the GRIP2 intron. Comparative genomic analyses revealed that RTOM1, RTOM2, and RTOM3 are strongly conserved and are under purifying selection between these species. These could be generated by tandem duplications before the divergence of platypus and echidna. All RTOM transcripts were specifically expressed in the testis, possibly suggesting their physiological importance. This is the first study reporting monotreme-specific de novo gene candidates derived from ERVs, which provides new insights into the unique evolution of monotremes.

Molecular evolution of vision-related genes may contribute to marsupial photic niche adaptations

Vision plays an essential role in the life of many animals. While most mammals are night-active (nocturnal), many have adapted to novel light environments. This includes diurnal (day-active) and crepuscular (twilight-active) species. Here, we used integrative approaches to investigate the molecular evolution of 112 vision-related genes across 19 genomes representing most marsupial orders. We found that four genes (GUCA1B, GUCY2F, RGR, and SWS2) involved in retinal phototransduction likely became functionally redundant in the ancestor of marsupials, a group of largely obligate nocturnal mammals. We also show evidence of rapid evolution and positive selection of bright-light vision genes in the common ancestor of Macropus (kangaroos, wallaroos, and wallabies). Macropus-specific amino acid substitutions in opsin genes (LWS and SWS1), in particular, may be an adaptation for crepuscular vision in this genus via opsin spectral sensitivity tuning. Our study set the stage for functional genetics studies and provides a stepping stone to future research efforts that fully capture the visual repertoire of marsupials.

Activation of Xist by an evolutionarily conserved function of KDM5C demethylase

XX female and XY male therian mammals equalize X-linked gene expression through the mitotically-stable transcriptional inactivation of one of the two X chromosomes in female somatic cells. Here, we describe an essential function of the X-linked homolog of an ancestral X-Y gene pair, Kdm5c-Kdm5d, in the expression of Xist lncRNA, which is required for stable X-inactivation. Ablation of Kdm5c function in females results in a significant reduction in Xist RNA expression. Kdm5c encodes a demethylase that enhances Xist expression by converting histone H3K4me2/3 modifications into H3K4me1. Ectopic expression of mouse and human KDM5C, but not the Y-linked homolog KDM5D, induces Xist in male mouse embryonic stem cells (mESCs). Similarly, marsupial (opossum) Kdm5c but not Kdm5d also upregulates Xist in male mESCs, despite marsupials lacking Xist, suggesting that the KDM5C function that activates Xist in eutherians is strongly conserved and predates the divergence of eutherian and metatherian mammals. In support, prototherian (platypus) Kdm5c also induces Xist in male mESCs. Together, our data suggest that eutherian mammals co-opted the ancestral demethylase KDM5C during sex chromosome evolution to upregulate Xist for the female-specific induction of X-inactivation.

The Reaction Specificity of Mammalian ALOX15 Orthologs is Changed During Late Primate Evolution and These Alterations Might Offer Evolutionary Advantages for Hominidae

Arachidonic acid lipoxygenases (ALOXs) have been implicated in the immune response of mammals. The reaction specificity of these enzymes is decisive for their biological functions and ALOX classification is based on this enzyme property. Comparing the amino acid sequences and the functional properties of selected mammalian ALOX15 orthologs we previously hypothesized that the reaction specificity of these enzymes can be predicted based on their amino acid sequences (Triad Concept) and that mammals, which are ranked in evolution below gibbons, express arachidonic acid 12-lipoxygenating ALOX15 orthologs. In contrast, Hominidae involving the great apes and humans possess 15-lipoxygenating enzymes (Evolutionary Hypothesis). These two hypotheses were based on sequence data of some 60 mammalian ALOX15 orthologs and about half of them were functionally characterized. Here, we compared the ALOX15 sequences of 152 mammals representing all major mammalian subclades expressed 44 novel ALOX15 orthologs and performed extensive mutagenesis studies of their triad determinants. We found that ALOX15 genes are absent in extant Prototheria but that corresponding enzymes frequently occur in Metatheria and Eutheria. More than 90% of them catalyze arachidonic acid 12-lipoxygenation and the Triad Concept is applicable to all of them. Mammals ranked in evolution above gibbons express arachidonic acid 15-lipoxygenating ALOX15 orthologs but enzymes with similar specificity are only present in less than 5% of mammals ranked below gibbons. This data suggests that ALOX15 orthologs have been introduced during Prototheria-Metatheria transition and put the Triad Concept and the Evolutionary Hypothesis on a much broader and more reliable experimental basis.

Why are you talking with snakes? To get new evolutionary insights in cardiac electrophysiology!

Brette, Le Guennec, and Thireau discuss recent findings on evolutionary cardiac electrophysiology.

The Evolution of Viviparity in Vertebrates

In the vertebrate tree of life, viviparity or live birth has independently evolved many times, resulting in a rich diversity of reproductive strategies. Viviparity is believed to be a mode of reproduction that evolved from the ancestral condition of oviparity or egg laying, where most of the fetal development occurs outside the body. Today, there is not a simple model of parity transition to explain this species-specific divergence in modes of reproduction. Most evidence points to a gradual series of evolutionary adaptations that account for this phenomenon of reproduction, elegantly displayed by various viviparous squamates that exhibit placentae formed by the appositions of maternal and embryonic tissues, which share significant homology with the tissues that form the placenta in therian mammals. In an era where the genomes of many vertebrate species are becoming available, studies are now exploring the molecular basis of this transition from oviparity to viviparity, and in some rare instances its possible reversibility, such as the Australian three-toed skink (Saiphos equalis). In contrast to the parity diversity in squamates, mammals are viviparous with the notable exception of the egg-laying monotremes. Advancing computational tools coupled with increasing genome availability across species that utilize different reproductive strategies promise to reveal the molecular underpinnings of the ancestral transition of oviparity to viviparity. As a result, the dramatic changes in reproductive physiology and anatomy that accompany these parity changes can be reinterpreted. This chapter will briefly explore the vertebrate modes of reproduction using a phylogenetic framework and where possible highlight the role of potential candidate genes that may help explain the polygenic origins of live birth.

Population analysis of retrotransposons in giraffe genomes supports RTE decline and widespread LINE1 activity in Giraffidae

BACKGROUND

The majority of structural variation in genomes is caused by insertions of transposable elements (TEs). In mammalian genomes, the main TE fraction is made up of autonomous and non-autonomous non-LTR retrotransposons commonly known as LINEs and SINEs (Long and Short Interspersed Nuclear Elements). Here we present one of the first population-level analysis of TE insertions in a non-model organism, the giraffe. Giraffes are ruminant artiodactyls, one of the few mammalian groups with genomes that are colonized by putatively active LINEs of two different clades of non-LTR retrotransposons, namely the LINE1 and RTE/BovB LINEs as well as their associated SINEs. We analyzed TE insertions of both types, and their associated SINEs in three giraffe genome assemblies, as well as across a population level sampling of 48 individuals covering all extant giraffe species.

RESULTS

The comparative genome screen identified 139,525 recent LINE1 and RTE insertions in the sampled giraffe population. The analysis revealed a drastically reduced RTE activity in giraffes, whereas LINE1 is still actively propagating in the genomes of extant (sub)-species. In concert with the extremely low activity of the giraffe RTE, we also found that RTE-dependent SINEs, namely Bov-tA and Bov-A2, have been virtually immobile in the last 2 million years. Despite the high current activity of the giraffe LINE1, we did not find evidence for the presence of currently active LINE1-dependent SINEs. TE insertion heterozygosity rates differ among the different (sub)-species, likely due to divergent population histories.

CONCLUSIONS

The horizontally transferred RTE/BovB and its derived SINEs appear to be close to inactivation and subsequent extinction in the genomes of extant giraffe species. This is the first time that the decline of a TE family has been meticulously analyzed from a population genetics perspective. Our study shows how detailed information about past and present TE activity can be obtained by analyzing large-scale population-level genomic data sets.

Nine Levels of Explanation : A Proposed Expansion of Tinbergen's Four-Level Framework for Understanding the Causes of Behavior

Tinbergen's classic "On Aims and Methods of Ethology" (Zeitschrift für Tierpsychologie, 20, 1963) proposed four levels of explanation of behavior, which he thought would soon apply to humans. This paper discusses the need for multilevel explanation; Huxley and Mayr's prior models, and others that followed; Tinbergen's differences with Lorenz on "the innate"; and Mayr's ultimate/proximate distinction. It synthesizes these approaches with nine levels of explanation in three categories: phylogeny, natural selection, and genomics (ultimate causes); maturation, sensitive period effects, and routine environmental effects (intermediate causes); and hormonal/metabolic processes, neural circuitry, and eliciting stimuli (proximate causes), as a respectful extension of Tinbergen's levels. The proposed classification supports and builds on Tinbergen's multilevel model and Mayr's ultimate/proximate continuum, adding intermediate causes in accord with Tinbergen's emphasis on ontogeny. It requires no modification of Standard Evolutionary Theory or The Modern Synthesis, but shows that much that critics claim was missing was in fact part of Neo-Darwinian theory (so named by J. Mark Baldwin in The American Naturalist in 1896) all along, notably reciprocal causation in ontogeny, niche construction, cultural evolution, and multilevel selection. Updates of classical examples in ethology are offered at each of the nine levels, including the neuroethological and genomic findings Tinbergen foresaw. Finally, human examples are supplied at each level, fulfilling his hope of human applications as part of the biology of behavior. This broad ethological framework empowers us to explain human behavior-eventually completely-and vindicates the idea of human nature, and of humans as a part of nature.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Origin and Isoform Specific Functions of Exchange Proteins Directly Activated by cAMP: A Phylogenetic Analysis

Exchange proteins directly activated by cAMP (EPAC1 and EPAC2) are one of the several families of cellular effectors of the prototypical second messenger cAMP. To understand the origin and molecular evolution of EPAC proteins, we performed a comprehensive phylogenetic analysis of EPAC1 and EPAC2. Our study demonstrates that unlike its cousin PKA, EPAC proteins are only present in multicellular Metazoa. Within the EPAC family, EPAC1 is only associated with chordates, while EPAC2 spans the entire animal kingdom. Despite a much more contemporary origin, EPAC1 proteins show much more sequence diversity among species, suggesting that EPAC1 has undergone more selection and evolved faster than EPAC2. Phylogenetic analyses of the individual cAMP binding domain (CBD) and guanine nucleotide exchange (GEF) domain of EPACs, two most conserved regions between the two isoforms, further reveal that EPAC1 and EPAC2 are closely clustered together within both the larger cyclic nucleotide receptor and RAPGEF families. These results support the notion that EPAC1 and EPAC2 share a common ancestor resulting from a fusion between the CBD of PKA and the GEF from RAPGEF1. On the other hand, the two terminal extremities and the RAS-association (RA) domains show the most sequence diversity between the two isoforms. Sequence diversities within these regions contribute significantly to the isoform-specific functions of EPACs. Importantly, unique isoform-specific sequence motifs within the RA domain have been identified.

PEG10 viral aspartic protease domain is essential for the maintenance of fetal capillary structure in the mouse placenta

The therian-specific gene paternally expressed 10 (Peg10) plays an essential role in placenta formation: Peg10 knockout mice exhibit early embryonic lethality as a result of severe placental defects. The PEG10 protein exhibits homology with long terminal repeat (LTR) retrotransposon GAG and POL proteins; therefore, we generated mice harboring a mutation in the highly conserved viral aspartic protease motif in the POL-like region of PEG10 because this motif is essential for the life cycle of LTR retrotransposons/retroviruses. Intriguingly, frequent perinatal lethality, not early embryonic lethality, was observed with fetal and placental growth retardation starting mid-gestation. In the mutant placentas, severe defects were observed in the fetal vasculature, where PEG10 is expressed in the three trophoblast cell layers that surround fetal capillary endothelial cells. Thus, Peg10 has essential roles, not only in early placenta formation, but also in placental vasculature maintenance from mid- to late-gestation. This implies that along the feto-maternal placenta interface an interaction occurs between two retrovirus-derived genes, Peg10 and retrotransposon Gag like 1 (Rtl1, also called Peg11), that is essential for the maintenance of fetal capillary endothelial cells.

Summary: Disruption of the highly conserved viral aspartic protease domain in PEG10 causes placental abnormality leading to perinatal lethality in mice.

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Evolution of the Avian Eggshell Biomineralization Protein Toolkit - New Insights From Multi-Omics

The avian eggshell is a remarkable biomineral, which is essential for avian reproduction; its properties permit embryonic development in the desiccating terrestrial environment, and moreover, are critically important to preserve unfertilized egg quality for human consumption. This calcium carbonate (CaCO₃) bioceramic is made of 95% calcite and 3.5% organic matrix; it protects the egg contents against microbial penetration and mechanical damage, allows gaseous exchange, and provides calcium for development of the embryonic skeleton. In vertebrates, eggshell occurs in the Sauropsida and in a lesser extent in Mammalia taxa; avian eggshell calcification is one of the fastest known CaCO₃ biomineralization processes, and results in a material with excellent mechanical properties. Thus, its study has triggered a strong interest from the researcher community. The investigation of eggshell biomineralization in birds over the past decades has led to detailed characterization of its protein and mineral constituents. Recently, our understanding of this process has been significantly improved using high-throughput technologies (i.e., proteomics, transcriptomics, genomics, and bioinformatics). Presently, more or less complete eggshell proteomes are available for nine birds, and therefore, key proteins that comprise the eggshell biomineralization toolkit are beginning to be identified. In this article, we review current knowledge on organic matrix components from calcified eggshell. We use these data to analyze the evolution of selected matrix proteins and underline their role in the biological toolkit required for eggshell calcification in avian species. Amongst the panel of eggshell-associated proteins, key functional domains are present such as calcium-binding, vesicle-binding and protein-binding. These technical advances, combined with progress in mineral ultrastructure analyses, have opened the way for new hypotheses of mineral nucleation and crystal growth in formation of the avian eggshell, including transfer of amorphous CaCO₃ in vesicles from uterine cells to the eggshell mineralization site. The enrichment of multi-omics datasets for bird species is critical to understand the evolutionary context for development of CaCO₃ biomineralization in metazoans, leading to the acquisition of the robust eggshell in birds (and formerly dinosaurs).

Platypus and echidna genomes reveal mammalian biology and evolution

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

New Insights into the Mammalian Egg Zona Pellucida

Mammalian oocytes are surrounded by an extracellular coat called the zona pellucida (ZP), which, from an evolutionary point of view, is the most ancient of the coats that envelope vertebrate oocytes and conceptuses. This matrix separates the oocyte from cumulus cells and is responsible for species-specific recognition between gametes, preventing polyspermy and protecting the preimplantation embryo. The ZP is a dynamic structure that shows different properties before and after fertilization. Until very recently, mammalian ZP was believed to be composed of only three glycoproteins, ZP1, ZP2 and ZP3, as first described in mouse. However, studies have revealed that this composition is not necessarily applicable to other mammals. Such differences can be explained by an analysis of the molecular evolution of the ZP gene family, during which ZP genes have suffered pseudogenization and duplication events that have resulted in differing models of ZP protein composition. The many discoveries made in recent years related to ZP composition and evolution suggest that a compilation would be useful. Moreover, this review analyses ZP biosynthesis, the role of each ZP protein in different mammalian species and how these proteins may interact among themselves and with other proteins present in the oviductal lumen.

Origin and Expansion of the Serine Protease Repertoire in the Myelomonocyte Lineage

The deepest evolutionary branches of the trypsin/chymotrypsin family of serine proteases are represented by the digestive enzymes of the gastrointestinal tract and the multi-domain proteases of the blood coagulation and complement system. Similar to the very old digestive system, highly diverse cleavage specificities emerged in various cell lineages of the immune defense system during vertebrate evolution. The four neutrophil serine proteases (NSPs) expressed in the myelomonocyte lineage, neutrophil elastase, proteinase 3, cathepsin G, and neutrophil serine protease 4, collectively display a broad repertoire of (S1) specificities. The origin of NSPs can be traced back to a circulating liver-derived trypsin-like protease, the complement factor D ancestor, whose activity is tightly controlled by substrate-induced activation and TNFα-induced locally upregulated protein secretion. However, the present-day descendants are produced and converted to mature enzymes in precursor cells of the bone marrow and are safely sequestered in granules of circulating neutrophils. The potential site and duration of action of these cell-associated serine proteases are tightly controlled by the recruitment and activation of neutrophils, by stimulus-dependent regulated secretion of the granules, and by various soluble inhibitors in plasma, interstitial fluids, and in the inflammatory exudate. An extraordinary dynamic range and acceleration of immediate defense responses have been achieved by exploiting the high structural plasticity of the trypsin fold.

Draft Genome Assemblies and Annotations of Agrypnia vestita Walker, and Hesperophylax magnus Banks Reveal Substantial Repetitive Element Expansion in Tube Case-Making Caddisflies (Insecta: Trichoptera)

Trichoptera (caddisflies) play an essential role in freshwater ecosystems; for instance, larvae process organic material from the water and are food for a variety of predators. Knowledge on the genomic diversity of caddisflies can facilitate comparative and phylogenetic studies thereby allowing scientists to better understand the evolutionary history of caddisflies. Although Trichoptera are the most diverse aquatic insect order, they remain poorly represented in terms of genomic resources. To date, all long-read based genomes have been sequenced from individuals in the retreat-making suborder, Annulipalpia, leaving ∼275 Ma of evolution without high-quality genomic resources. Here, we report the first long-read based de novo genome assemblies of two tube case-making Trichoptera from the suborder Integripalpia, Agrypnia vestita Walker and Hesperophylax magnus Banks. We find that these tube case-making caddisflies have genome sizes that are at least 3-fold larger than those of currently sequenced annulipalpian genomes and that this pattern is at least partly driven by major expansion of repetitive elements. In H. magnus, long interspersed nuclear elements alone exceed the entire genome size of some annulipalpian counterparts suggesting that caddisflies have high potential as a model for understanding genome size evolution in diverse insect lineages.

Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution

Members of the APOBEC family of cytidine deaminases show antiviral activities in mammalian cells through lethal editing in the genomes of small DNA viruses, herpesviruses and retroviruses, and potentially those of RNA viruses such as coronaviruses. Consistent with the latter, APOBEC-like directional C→U transitions of genomic plus-strand RNA are greatly overrepresented in SARS-CoV-2 genome sequences of variants emerging during the COVID-19 pandemic. A C→U mutational process may leave evolutionary imprints on coronavirus genomes, including extensive homoplasy from editing and reversion at targeted sites and the occurrence of driven amino acid sequence changes in viral proteins. If sustained over longer periods, this process may account for the previously reported marked global depletion of C and excess of U bases in human seasonal coronavirus genomes. This review synthesizes the current knowledge on APOBEC evolution and function and the evidence of their role in APOBEC-mediated genome editing of SARS-CoV-2 and other coronaviruses.

-

SARS-CoV-2 sequence variants contain an overabundance of C- > U transitions

-

C- > U transitions are the hallmark of the activity of APOBEC cytosine deaminases

-

Further work is needed to determine APOBEC's role in coronavirus evolution

An Investigation of ZZ/ZW and XX/XY Sex Determination Systems in North African Catfish (Clarias gariepinus, )

An investigation of sex-specific loci may provide important insights into fish sex determination strategies. This may be useful for biotechnological purposes, for example, to produce all-male or all-female fish for commercial breeding. The North African catfish species, Clarias gariepinus, has been widely adopted for aquaculture because its superior growth and disease resistance render the species suitable for hybridization with other catfish to improve the productivity and quality of fish meat. This species has either a ZZ/ZW or XX/XY sex determination system. Here, we investigate and characterize these systems using high-throughput genome complexity reduction sequencing as Diversity Arrays Technology. This approach was effective in identifying moderately sex-linked loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 30 perfectly sexed individuals of C. gariepinus. However, SNPs based markers were not found in this study. In total, 41 loci met the criteria for being moderately male-linked (with male vs. female ratios 80:20 and 70:30), while 25 loci were found to be moderately linked to female sex. No strictly male- or female-linked loci were detected. Seven moderately male-linked loci were partially homologous to some classes of transposable elements and three moderately male-linked loci were partially homologous to functional genes. Our data showed that the male heterogametic XX/XY sex determination system should co-exist with the ZZ/ZW system in C. gariepinus. Our finding of the co-existence of XX/XY and ZZ/ZW systems can be applied to benefit commercial breeding of this species in Thailand. This approach using moderately sex-linked loci provides a solid baseline for revealing sex determination mechanisms and identify potential sex determination regions in catfish, allowing further investigation of genetic improvements in breeding programs.

Evolution of reproductive strategies in the species-rich land snail subfamily Phaedusinae (Stylommatophora: Clausiliidae)

Most of the present knowledge on animal reproductive mode evolution, and possible factors driving transitions between oviparity and viviparity is based on studies on vertebrates. The species rich door snail (Clausiliidae) subfamily Phaedusinae represents a suitable and unique model for further examining parity evolution, as three different strategies, oviparity, viviparity, and the intermediate mode of embryo-retention, occur in this group. The present study reconstructs the evolution of reproductive strategies in Phaedusinae based on time-calibrated molecular phylogenetics, reproductive mode examinations and ancestral state reconstruction. Our phylogenetic analysis employing multiple mitochondrial and nuclear markers identified a well-supported clade (including the tribes Phaedusini and Serrulinini) that contains species exhibiting various reproductive strategies. This clade evolved from an oviparous most recent common ancestor according to our reconstruction. All non-oviparous taxa are confined to a highly supported subclade, coinciding with the tribe Phaedusini. Both oviparity and viviparity occur frequently in different lineages of this subclade that are not closely related. During Phaedusini diversification, multiple transitions in reproductive strategy must have taken place, which could have been promoted by a high fitness of embryo-retaining species. The evolutionary success of this group might result from the maintenance of various strategies.

Molecular insights into the divergence and diversity of post-testicular maturation strategies

Competition to achieve paternity has coerced the development of a multitude of male reproductive strategies. In one of the most well-studied examples, the spermatozoa of all mammalian species must undergo a series of physiological changes as they transit the male (epididymal maturation) and female (capacitation) reproductive tracts prior to realizing their potential to fertilize an ovum. However, the origin and adaptive advantage afforded by these intricate processes of post-testicular sperm maturation remain to be fully elucidated. Here, we review literature pertaining to the nature and the physiological role of epididymal maturation and subsequent capacitation in comparative vertebrate taxa including representative species from the avian, reptilian, and mammalian lineages. Such insights are discussed in terms of the framework they provide for helping to understand the evolutionary significance of post-testicular sperm maturation.

Comprehensive genomic analysis reveals dynamic evolution of endogenous retroviruses that code for retroviral-like protein domains

Background

Endogenous retroviruses (ERVs) are remnants of ancient retroviral infections of mammalian germline cells. A large proportion of ERVs lose their open reading frames (ORFs), while others retain them and become exapted by the host species. However, it remains unclear what proportion of ERVs possess ORFs (ERV-ORFs), become transcribed, and serve as candidates for co-opted genes.

Results

We investigated characteristics of 176,401 ERV-ORFs containing retroviral-like protein domains (gag, pro, pol, and env) in 19 mammalian genomes. The fractions of ERVs possessing ORFs were overall small (~ 0.15%) although they varied depending on domain types as well as species. The observed divergence of ERV-ORF from their consensus sequences showed bimodal distributions, suggesting that a large proportion of ERV-ORFs either recently, or anciently, inserted themselves into mammalian genomes. Alternatively, very few ERVs lacking ORFs were found to exhibit similar divergence patterns. To identify candidates for ERV-derived genes, we estimated the ratio of non-synonymous to synonymous substitution rates (dN/dS) for ERV-ORFs in human and non-human mammalian pairs, and found that approximately 42% of the ERV-ORFs showed dN/dS < 1. Further, using functional genomics data including transcriptome sequencing, we determined that approximately 9.7% of these selected ERV-ORFs exhibited transcriptional potential.

Conclusions

These results suggest that purifying selection operates on a certain portion of ERV-ORFs, some of which may correspond to uncharacterized functional genes hidden within mammalian genomes. Together, our analyses suggest that more ERV-ORFs may be co-opted in a host-species specific manner than we currently know, which are likely to have contributed to mammalian evolution and diversification.

Transcriptomic Analysis of Four Cerianthid (Cnidaria, Ceriantharia) Venoms

Tube anemones, or cerianthids, are a phylogenetically informative group of cnidarians with complex life histories, including a pelagic larval stage and tube-dwelling adult stage, both known to utilize venom in stinging-cell rich tentacles. Cnidarians are an entirely venomous group that utilize their proteinaceous-dominated toxins to capture prey and defend against predators, in addition to several other ecological functions, including intraspecific interactions. At present there are no studies describing the venom for any species within cerianthids. Given their unique development, ecology, and distinct phylogenetic-placement within Cnidaria, our objective is to evaluate the venom-like gene diversity of four species of cerianthids from newly collected transcriptomic data. We identified 525 venom-like genes between all four species. The venom-gene profile for each species was dominated by enzymatic protein and peptide families, which is consistent with previous findings in other cnidarian venoms. However, we found few toxins that are typical of sea anemones and corals, and furthermore, three of the four species express toxin-like genes closely related to potent pore-forming toxins in box jellyfish. Our study is the first to provide a survey of the putative venom composition of cerianthids and contributes to our general understanding of the diversity of cnidarian toxins.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Marsupial and monotreme milk-a review of its nutrient and immune properties

All mammals are characterized by the ability of females to produce milk. Marsupial (metatherian) and monotreme (prototherian) young are born in a highly altricial state and rely on their mother’s milk for the first part of their life. Here we review the role and importance of milk in marsupial and monotreme development. Milk is the primary source of sustenance for young marsupials and monotremes and its composition varies at different stages of development. We applied nutritional geometry techniques to a limited number of species with values available to analyze changes in macronutrient composition of milk at different stages. Macronutrient energy composition of marsupial milk varies between species and changes concentration during the course of lactation. As well as nourishment, marsupial and monotreme milk supplies growth and immune factors. Neonates are unable to mount a specific immune response shortly after birth and therefore rely on immunoglobulins, immunological cells and other immunologically important molecules transferred through milk. Milk is also essential to the development of the maternal-young bond and is achieved through feedback systems and odor preferences in eutherian mammals. However, we have much to learn about the role of milk in marsupial and monotreme mother-young bonding. Further research is warranted in gaining a better understanding of the role of milk as a source of nutrition, developmental factors and immunity, in a broader range of marsupial species, and monotremes.

Macroevolutionary Analyses Suggest That Environmental Factors, Not Venom Apparatus, Play Key Role in Terebridae Marine Snail Diversification

How species diversification occurs remains an unanswered question in predatory marine invertebrates, such as sea snails of the family Terebridae. However, the anatomical disparity found throughput the Terebridae provides a unique perspective for investigating diversification patterns in venomous predators. In this study, a new dated molecular phylogeny of the Terebridae is used as a framework for investigating diversification of the family through time, and for testing the putative role of intrinsic and extrinsic traits, such as shell size, larval ecology, bathymetric distribution, and anatomical features of the venom apparatus, as drivers of terebrid species diversification. Macroevolutionary analysis revealed that when diversification rates do not vary across Terebridae clades, the whole family has been increasing its global diversification rate since 25 Ma. We recovered evidence for a concurrent increase in diversification of depth ranges, while shell size appeared to have undergone a fast divergence early in terebrid evolutionary history. Our data also confirm that planktotrophy is the ancestral larval ecology in terebrids, and evolutionary modeling highlighted that shell size is linked to larval ecology of the Terebridae, with species with long-living pelagic larvae tending to be larger and have a broader size range than lecithotrophic species. Although we recovered patterns of size and depth trait diversification through time and across clades, the presence or absence of a venom gland (VG) did not appear to have impacted Terebridae diversification. Terebrids have lost their venom apparatus several times and we confirm that the loss of a VG happened in phylogenetically clustered terminal taxa and that reversal is extremely unlikely. Our findings suggest that environmental factors, and not venom, have had more influence on terebrid evolution.