Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Genome-wide sequence divergence of satellite DNA could underlie meiotic failure in male hybrids of bighead catfish and North African catfish (Clarias, Clariidae).

Hybrid sterility, a hallmark of postzygotic isolation, arises from parental genome divergence disrupting meiosis. While chromosomal incompatibility is often implicated, the underlying mechanisms remain unclear. This study investigated meiotic behavior and genome-wide divergence in bighead catfish (C. macrocephalus), North African catfish (C. gariepinus), and their sterile male hybrids (important in aquaculture). Repetitive DNA analysis using bioinformatics and cytogenetics revealed significant divergence in satellite DNA (satDNA) families between parental species. Notably, one hybrid exhibited successful meiosis and spermatozoa production, suggesting potential variation in sterility expression. Our findings suggest that genome-wide satDNA divergence, rather than chromosome number differences, likely contributes to meiotic failure and male sterility in these catfish hybrids.

Testing phylogenetic signal with categorical traits and tree uncertainty.

SUMMARY: The phylogenetic signal, frequently used to identify signatures of adaptive evolution or important associations between genes and phenotypes, measures the tendency for recently diverged species to resemble each other more than distantly related species. An example of such a measure is the δ statistic, which uses Shannon entropy to measure the degree of phylogenetic signal between a categorical trait and a phylogeny. In this study, we refined this statistic to account for tree uncertainty, resulting in more accurate assessments of phylogenetic associations. In addition, we provided a more accessible and computationally efficient implementation of the δ statistic that will facilitate its use by the evolutionary community. AVAILABILITY AND IMPLEMENTATION: github.com/diogo-s-ribeiro/delta-statistic.

MicrosatNavigator: exploring nonrandom distribution and lineage-specificity of microsatellite repeat motifs on vertebrate sex chromosomes across 186 whole genomes.

Microsatellites are short tandem DNA repeats, ubiquitous in genomes. They are believed to be under selection pressure, considering their high distribution and abundance beyond chance or random accumulation. However, limited analysis of microsatellites in single taxonomic groups makes it challenging to understand their evolutionary significance across taxonomic boundaries. Despite abundant genomic information, microsatellites have been studied in limited contexts and within a few species, warranting an unbiased examination of their genome-wide distribution in distinct versus closely related-clades. Large-scale comparisons have revealed relevant trends, especially in vertebrates. Here, "MicrosatNavigator", a new tool that allows quick and reliable investigation of perfect microsatellites in DNA sequences, was developed. This tool can identify microsatellites across the entire genome sequences. Using this tool, microsatellite repeat motifs were identified in the genome sequences of 186 vertebrates. A significant positive correlation was noted between the abundance, density, length, and GC bias of microsatellites and specific lineages. The (AC)n motif is the most prevalent in vertebrate genomes, showing distinct patterns in closely related species. Longer microsatellites were observed on sex chromosomes in birds and mammals but not on autosomes. Microsatellites on sex chromosomes of non-fish vertebrates have the lowest GC content, whereas high-GC microsatellites (≥ 50 M% GC) are preferred in bony and cartilaginous fishes. Thus, similar selective forces and mutational processes may constrain GC-rich microsatellites to different clades. These findings should facilitate investigations into the roles of microsatellites in sex chromosome differentiation and provide candidate microsatellites for functional analysis across the vertebrate evolutionary spectrum.

Genomic Signatures of Divergent Ecological Strategies in a Recent Radiation of Neotropical Wild Cats.

Ecological differentiation among diverging species is an important component of the evolutionary process and can be investigated in rapid and recent radiations. Here, we use whole genome sequences of five species from the genus Leopardus, a recently diversified Neotropical lineage with species bearing distinctive morphological, ecological, and behavioral features, to investigate genome-wide diversity, comparative demographic history and signatures of positive selection. Our results show that divergent ecological strategies are reflected in genomic features, for example a generalist species shows historically larger effective population size and higher heterozygosity than habitat specialists. The demographic history of these cats seems to have been jointly driven by climate fluctuations and habitat specialization, with different ecological adaptations leading to distinct trajectories. Finally, a gene involved in vertebrate retinal neurogenesis (POU4F2) was found to be under positive selection in the margay, a cat with notoriously large eyes that are likely associated with its nocturnal and arboreal specializations.

The evolutionary history of extinct and living lions.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes.

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection.

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.

Biomedical Potential of the Neglected Molluscivorous and Vermivorous Conus Species.

Within the Conidae family, the piscivorous Conus species have been a hotspot target for drug discovery. Here, we assess the relevance of Conus and their other feeding habits, and thus under distinctive evolutionary constraints, to highlight the potential of neglected molluscivorous and vermivorous species in biomedical research and pharmaceutical industry. By singling out the areas with inadequate Conus disquisition, such as the Tamil Nadu Coast and the Andaman Islands, research resources can be expanded and better protected through awareness. In this study, 728 Conus species and 190 species from three other genera (1 from Californiconus, 159 from Conasprella and 30 from Profundiconus) in the Conidae family are assessed. The phylogenetic relationships of the Conidae species are determined and their known feeding habits superimposed. The worm-hunting species appeared first, and later the mollusc- and fish-hunting species were derived independently in the Neogene period (around 23 million years ago). Interestingly, many Conus species in the warm and shallow waters become polyphagous, allowing them to hunt both fish and worms, given the opportunities. Such newly gained trait is multi originated. This is controversial, given the traditional idea that most Conus species are specialized to hunt certain prey categories. However, it shows the functional complexity and great potential of conopeptides from some worm-eating species. Pharmaceutical attempts and relevant omics data have been differentially obtained. Indeed, data from the fish-hunting species receive strong preference over the worm-hunting ones. Expectedly, conopeptides from the fish-hunting species are believed to include the most potential candidates for biomedical research. Our work revisits major findings throughout the Conus evolution and emphasizes the importance of increasing omics surveys complemented with further behavior observation studies. Hence, we claim that Conus species and their feeding habits are equally important, highlighting many places left for Conus exploration worldwide. We also discuss the Conotoxin drug discovery potentials and the urgency of protecting the bioresources of Conus species. In particular, some vermivorous species have demonstrated great potential in malaria therapy, while other conotoxins from several worm- and mollusc-eating species exhibited explicit correlation with SARS-CoV-2. Reclaiming idle data with new perspectives could also promote interdisciplinary studies in both virological and toxicological fields.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Symbiotic Associations in Ascidians: Relevance for Functional Innovation and Bioactive Potential.

Associations between different organisms have been extensively described in terrestrial and marine environments. These associations are involved in roles as diverse as nutrient exchanges, shelter or adaptation to adverse conditions. Ascidians are widely dispersed marine invertebrates associated to invasive behaviours. Studying their microbiomes has interested the scientific community, mainly due to its potential for bioactive compounds production-e.g., ET-73 (trabectedin, Yondelis), an anticancer drug. However, these symbiotic interactions embrace several environmental and biological functions with high ecological relevance, inspiring diverse biotechnological applications. We thoroughly reviewed microbiome studies (microscopic to metagenomic approaches) of around 171 hosts, worldwide dispersed, occurring at different domains of life (Archaea, Bacteria, Eukarya), to illuminate the functions and bioactive potential of associated organisms in ascidians. Associations with Bacteria are the most prevalent, namely with Cyanobacteria, Proteobacteria, Bacteroidetes, Actinobacteria and Planctomycetes phyla. The microbiomes of ascidians belonging to Aplousobranchia order have been the most studied. The integration of worldwide studies characterizing ascidians' microbiome composition revealed several functions including UV protection, bioaccumulation of heavy metals and defense against fouling or predators through production of natural products, chemical signals or competition. The critical assessment and characterization of these communities is extremely valuable to comprehend their biological/ecological role and biotechnological potential.

The genomic context of retrocopies increases their chance of functional relevancy in mammals.

Described as "junk" DNA, pseudogenes are dead structures of previously active genes present in genomes. Pseudogenes are categorized into two main classes: processed pseudogenes, formed through retrotransposition, and non-processed pseudogenes, typically originated from gene decay following duplication events. The term "processed pseudogene" has changed to "retrocopy" since they are likely to evolve new functional roles and became a retrogene. Here, we surveyed 38,080 retrocopies from chimpanzee, dog, human, mouse, and rat genomes to assess their potential adaptive value. The retrocopies inserted in the same chromosome of the parental gene have higher chances of remain potentially "active" (absence of premature stop codons and frameshifts) (~26.1%), while those placed into a different chromosome have a twofold decrease chance of continuing potentially "active" (~7.52%). The genomic context of their placement seems associated with their expression. Retrocopies placed in intragenic regions and the same sense of the "host" gene have higher chances of being expressed relative to other genomic contexts. The proximity of retrocopies to their parental gene is associated with a lower decay rate, and their location likely influence their expression. Thus, despite their unclear role, retrocopies are probably involved in adaptive processes. Our results evidence natural selection acting in retrocopies.

Evolutionary genomics of mammalian lung cancer genes reveals signatures of positive selection in APC, RB1 and TP53.

Lung cancer is the type of cancer causing most deaths in humans, with 234,030 new cases of lung cancer diagnosed in the United States in 2018. Recently, Tumor suppressor genes (TSGs) or the control of its pathway became promising drug targets for cancer therapy. A diverse group of TSGs is involved in progression and metastasis of lung cancer. Here, we surveyed nine highly significant mutated genes in 20 mammalian genomes to assess signatures of adaptive evolution using maximum likelihood approaches. We found that three genes (APC, RB1, and TP53) are under strong positive selection, influencing amino acids located in functionally important protein domains, such as three sites in APC found in the APC_N_CC domain, which is responsible for the binding to beta-catenin armadillo repeats that regulate beta-catenin level (beta-catenin is a transcription factor and its misregulation lead to malignant transformation of normal cells). Such sites substitutions mostly increase the stability of the domain. Moreover, substitution of some other sites found in important motifs, such as codon 47 (proline-directed kinase motif) in TP53, modify the phosphorylation activity of TP53 playing a key role in cancer risk. Our findings will open recommendation to drug targeting sites and will foster further research to understand better these proteins function.

Acquisition of social behavior in mammalian lineages is related with duplication events of FPR genes.

Formyl peptide receptors (FPRs) were firstly detected in immune cells where they act as key mediators of leukocyte chemotaxis, promoting the host defense against pathogens. Recently, three paralogs were reported in Homo sapiens (FPR1-3) and seven paralogs in Mus musculus (FPR1, FPRrs1-4, FPRrs6 and FPRrs7), but information from other mammalian lineages is scarce, including ambiguities in the current nomenclature system (e.g. absence of an orthologous relation between human and mouse FPR3). Here, we explored the FPR gene repertoire across 175 mammalian genomes using integrative phylogenetic and synteny analyses to describe the evolutionary history of FPRs in all mammalian orders. FPRs present a well conserved synteny but showed dynamic episodes of duplication events specific to several mammalian orders (Chiroptera, Perissodactyla, Primates and Rodentia), with up to 11 paralogs in some cases. Despite FPRs could be expressed in a panoply of tissues, there is a suggestion that they maintain an exclusive immunological function. However, we observed that species with social behavior have higher repertoire of FPRs in contrast with species with solitary lifestyle. Such evidence suggests a strict relationship between the optimization of the immunological system (by FPR duplication patterns) and the mammalian social behavior.

Neofunctionalization of the UCP1 mediated the non-shivering thermogenesis in the evolution of small-sized placental mammals.

The acquisition of UCP1-mediated non-shivering thermogenesis (NST) was an important event during the evolution of mammals. Here, we assessed the thermogenic neofunctionalization that occurred in the mammalian UCP1, by performing detailed comparative evolutionary genomics analyses (including phylogenetic and selection analyses) of the UCP family members across all major vertebrate classes. Heterogeneously distributed positive selection signatures were found in several UCPs, being preferably located in the mitochondrial matrix domains. Additionally, comparisons with non-mammalian orthologs showed increased evolutionary rates of the mammalian UCP1, not observable in the phylogenetically related UCP2 and UCP3 paralogs. Also, parallel signatures of episodic positive selection (ω > 1) were found in the ancestral branches of both Glires (rodents and lagomorphs) and Afroinsectivores (afrosoricids and macroscelids), underlining the importance of the UCP1 thermogenic activity in these mammalian groups. Finally, we hypothesize that the independent positive selection events that occurred in these two lineages resulted in two UCP1-mediated NST approaches, namely the cold acute response in the Glires and the reproduction success enhancement in the Afroinsectivores.

Comparative eye and liver differentially expressed genes reveal monochromatic vision and cancer resistance in the shortfin mako shark (Isurus oxyrinchus).

The shortfin mako, Isurus oxyrinchus is an oceanic pelagic shark found worldwide in tropical and subtropical waters. However, the understanding of its biology at molecular level is still incipient. We sequenced the messenger RNA isolated from eye and liver tissues. De novo transcriptome yielded a total of 705,940 transcripts. A total of 3774 genes were differentially expressed (DEGs), with 1612 in the eye and 2162 in the liver. Most DEGs in the eye were related to structural and signaling functions, including nonocular and ocular opsin genes, whereas nine out of ten most overexpressed genes in the liver were related to tumor suppression, wound healing, and human diseases. Furthermore, DEGs findings provide insights on the monochromatic shark vision and a repertory of cancer-related genes, which may be insightful to elucidate shark resistance to cancer. Therefore, our results provide valuable sequence resources for future functional and population studies.

Proteogenomic Characterization of the Cement and Adhesive Gland of the Pelagic Gooseneck Barnacle Lepas anatifera.

We focus on the stalked goose barnacle L. anatifera adhesive system, an opportunistic less selective species for the substrate, found attached to a variety of floating objects at seas. Adhesion is an adaptative character in barnacles, ensuring adequate positioning in the habitat for feeding and reproduction. The protein composition of the cement multicomplex and adhesive gland was quantitatively studied using shotgun proteomic analysis. Overall, 11,795 peptide sequences were identified in the gland and 2206 in the cement, clustered in 1689 and 217 proteinGroups, respectively. Cement specific adhesive proteins (CPs), proteases, protease inhibitors, cuticular and structural proteins, chemical cues, and many unannotated proteins were found, among others. In the cement, CPs were the most abundant (80.5%), being the bulk proteins CP100k and -52k the most expressed of all, and CP43k-like the most expressed interfacial protein. Unannotated proteins comprised 4.7% of the cement proteome, ranking several of them among the most highly expressed. Eight of these proteins showed similar physicochemical properties and amino acid composition to known CPs and classified through Principal Components Analysis (PCA) as new CPs. The importance of PCA on the identification of unannotated non-conserved adhesive proteins, whose selective pressure is on their relative amino acid abundance, was demonstrated.

Measuring phylogenetic signal between categorical traits and phylogenies.

MOTIVATION: Determining whether a trait and phylogeny share some degree of phylogenetic signal is a flagship goal in evolutionary biology. Signatures of phylogenetic signal can assist the resolution of a broad range of evolutionary questions regarding the tempo and mode of phenotypic evolution. However, despite the considerable number of strategies to measure it, few and limited approaches exist for categorical traits. Here, we used the concept of Shannon entropy and propose the δ statistic for evaluating the degree of phylogenetic signal between a phylogeny and categorical traits. RESULTS: We validated δ as a measure of phylogenetic signal: the higher the δ-value the higher the degree of phylogenetic signal between a given tree and a trait. Based on simulated data we proposed a threshold-based classification test to pinpoint cases of phylogenetic signal. The assessment of the test's specificity and sensitivity suggested that the δ approach should only be applied to 20 or more species. We have further tested the performance of δ in scenarios of branch length and topology uncertainty, unbiased and biased trait evolution and trait saturation. Our results showed that δ may be applied in a wide range of phylogenetic contexts. Finally, we investigated our method in 14 360 mammalian gene trees and found that olfactory receptor genes are significantly associated with the mammalian activity patterns, a result that is congruent with expectations and experiments from the literature. Our application shows that δ can successfully detect molecular signatures of phenotypic evolution. We conclude that δ represents a useful measure of phylogenetic signal since many phenotypes can only be measured in categories. AVAILABILITY AND IMPLEMENTATION: https://github.com/mrborges23/delta_statistic. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.