Comparative painting reveals strong chromosome homology over 80 million years of bird evolution

Understanding microchromosomal organization and evolution in four representative woodpeckers (Picidae, Piciformes) through BAC-FISH analysis

The genome organization of woodpeckers has several distinctive features e.g., an uncommon accumulation of repetitive sequences, enlarged Z chromosomes, and atypical diploid numbers. Despite the large diversity of species, there is a paucity of detailed cytogenomic studies for this group and we thus aimed to rectify this. Genome organization patterns and hence evolutionary change in the microchromosome formation of four species (Colaptes campestris, Veniliornis spilogaster, Melanerpes candidus, and Picumnus nebulosus) was established through fluorescence in situ hybridization using bacterial artificial chromosomes originally derived from Gallus gallus and Taeniopygia guttata. Findings suggest that P. nebulosus (2n = 110), which was described for the first time, had the most basal karyotype among species of Picidae studied here, and probably arose as a result of fissions of avian ancestral macrochromosomes. We defined a new chromosomal number for V. spilogaster (2n = 88) and demonstrated microchromosomal rearrangements involving C. campestris plus a single, unique hitherto undescribed rearrangement in V. spilogaster. This comprised an inversion after a fusion involving the ancestral microchromosome 12 (homologous to chicken microchromosome 12). We also determined that the low diploid number of M. candidus is related to microchromosome fusions. Woodpeckers thus exhibit significantly rearranged karyotypes compared to the putative ancestral karyotype.

Evolution of bird sex chromosomes: a cytogenomic approach in Palaeognathae species

BACKGROUND

Different patterns of sex chromosome differentiation are seen in Palaeognathae birds, a lineage that includes the ratites (Struthioniformes, Rheiformes, Apterygiformes, Casuariiformes, and the sister group Tinamiformes). While some Tinamiform species have well-differentiated W chromosomes, both Z and W of all the flightless ratites are still morphologically undifferentiated. Here, we conducted a comprehensive analysis of the ZW differentiation in birds using a combination of cytogenetic, genomic, and bioinformatic approaches. The whole set of satDNAs from the emu (Dromaius novaehollandiae) was described and characterized. Furthermore, we examined the in situ locations of these satDNAs alongside several microsatellite repeats and carried out Comparative Genomic Hybridizations in two related species: the greater rhea (Rhea americana) and the tataupa tinamou (Crypturellus tataupa).

RESULTS

From the 24 satDNA families identified (which represent the greatest diversity of satDNAs ever uncovered in any bird species), only three of them were found to accumulate on the emu's sex chromosomes, with no discernible accumulation observed on the W chromosome. The W chromosomes of both the greater rhea and the emu did not exhibit a significant buildup of either C-positive heterochromatin or repetitive DNAs, indicating their large undifferentiation both at morphological and molecular levels. In contrast, the tataupa tinamou has a highly differentiated W chromosome that accumulates several DNA repeats.

CONCLUSION

The findings provide new information on the architecture of the avian genome and an inside look at the starting points of sex chromosome differentiation in birds.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

The Length Polymorphism of the 9th Intron in the Avian CHD1 Gene Allows Sex Determination in Some Species of Palaeognathae

In palaeognathous birds, several PCR-based methods and a range of genes and unknown genomic regions have been studied for the determination of sex. Many of these methods have proven to be unreliable, complex, expensive, and time-consuming. Even the most widely used PCR markers for sex typing in birds, the selected introns of the highly conserved CHD1 gene (primers P2/P8, 1237L/1272H, and 2550F/2718R), have rarely been effective in palaeognathous birds. In this study we used eight species of Palaeognathae to test three PCR markers: CHD1i9 (CHD1 gene intron 9) and NIPBLi16 (NIPBL gene intron 16) that performed properly as Psittaciformes sex differentiation markers, but have not yet been tested in Palaeognathae, as well as the CHD1iA intron (CHD1 gene intron 16), which so far has not been used effectively to sex palaeognathous birds. The results of our research indicate that the CHD1i9 marker effectively differentiates sex in four of the eight species we studied. In Rhea americana, Eudromia elegans, and Tinamus solitarius, the electrophoretic patterns of the amplicons obtained clearly indicate the sex of tested individuals, whereas in Crypturellus tataupa, sexing is possible based on poorly visible female specific bands. Additionally, we present and discuss the results of our in silico investigation on the applicability of CHD1i9 to sex other Palaeognathae that were not tested in this study.

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots

The karyotype of most birds has remained considerably stable during more than 100 million years’ evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

Parrots have undergone substantial karyotype evolution compared to most other birds. Here, Huang et al. analyze chromosome-level genome assemblies for four parrot species and elucidate the complex evolutionary history of parrot chromosomes.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Microchromosomes are building blocks of bird, reptile, and mammal chromosomes

Microchromosomes, once considered unimportant shreds of the chicken genome, are gene-rich elements with a high GC content and few transposable elements. Their origin has been debated for decades. We used cytological and whole-genome sequence comparisons, and chromosome conformation capture, to trace their origin and fate in genomes of reptiles, birds, and mammals. We find that microchromosomes as well as macrochromosomes are highly conserved across birds and share synteny with single small chromosomes of the chordate amphioxus, attesting to their origin as elements of an ancient animal genome. Turtles and squamates (snakes and lizards) share different subsets of ancestral microchromosomes, having independently lost microchromosomes by fusion with other microchromosomes or macrochromosomes. Patterns of fusions were quite different in different lineages. Cytological observations show that microchromosomes in all lineages are spatially separated into a central compartment at interphase and during mitosis and meiosis. This reflects higher interaction between microchromosomes than with macrochromosomes, as observed by chromosome conformation capture, and suggests some functional coherence. In highly rearranged genomes fused microchromosomes retain most ancestral characteristics, but these may erode over evolutionary time; surprisingly, de novo microchromosomes have rapidly adopted high interaction. Some chromosomes of early-branching monotreme mammals align to several bird microchromosomes, suggesting multiple microchromosome fusions in a mammalian ancestor. Subsequently, multiple rearrangements fueled the extraordinary karyotypic diversity of therian mammals. Thus, microchromosomes, far from being aberrant genetic elements, represent fundamental building blocks of amniote chromosomes, and it is mammals, rather than reptiles and birds, that are atypical.

Analysis of Sex Chromosome Evolution in the Clade Palaeognathae from Phased Genome Assembly

Birds in the clade Palaeognathae, excluding Tinamiformes, have morphologically conserved karyotypes and less differentiated ZW sex chromosomes compared with those of other birds. In particular, the sex chromosomes of the ostrich and emu have exceptionally large recombining pseudoautosomal regions (PARs), whereas non-PARs are classified into two strata according to the date of their origins: stratum 0 and stratum 1 (S1). However, the construction and analysis of the genome sequences in these regions in the clade Palaeognathae can be challenging because assembling the S1 region is difficult owing to low sequence diversity between gametologs (Z-linked and W-linked sequences). We addressed this issue by applying the Platanus-allee assembler and successfully constructed the haplotype-resolved (phased) assembly for female emu, cassowary, and ostrich using only sequence read data derived from the Illumina platform. Comparative genomic and phylogenetic analyses based on assembled Z-linked and W-linked sequences confirmed that the S1 region of emu and cassowary formed in their common ancestor. Moreover, the interspersed repetitive sequence landscapes in the S1 regions of female emu showed an expansion of younger repetitive elements in the W-linked S1 region, suggesting an interruption in homologous recombination in the S1 region. These results provide novel insights into the trajectory of sex chromosome evolution in the clade Palaeognathae and suggest that the Illumina-based phased assembly method is an effective approach for elucidating the evolutionary process underlying the transition from homomorphic to differentiated sex chromosomes.

Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

It is common that males and females display sexual dimorphisms, which usually result from sex-biased gene expression. Chinese hwamei (Garrulax canorus) is a good model for studying sex-biased gene expression because the song between the sexes is quite different. In this study, we analyze cerebrum and syrinx sex-biased gene expression and evolution using the de novo assembled Chinese hwamei transcriptome. In both the cerebrum and syrinx, our study revealed that most female-biased genes were actively expressed in females only, while most male-biased genes were actively expressed in both sexes. In addition, both male- and female-biased genes were enriched on the putative Z chromosome, suggesting the existence of sexually antagonistic genes and the insufficient dosage compensation of the Z-linked genes. We also identified a 9 Mb sex linkage region on the putative 4A chromosome which enriched more than 20% of female-biased genes. Resultantly, male-biased genes in both tissues had significantly higher Ka/Ks and effective number of codons (ENCs) than unbiased genes, and this suggested that male-biased genes which exhibit accelerated divergence may have resulted from positive selection. Taken together, our results initially revealed the reasons for the differences in singing behavior between males and females of Chinese hwamei.

Long-term stability of sex chromosome gene content allows accurate qPCR-based molecular sexing across birds

Embryos, juveniles, and even adults of many bird species lack pronounced external sexually dimorphic characteristics. Accurate identification of sex is crucial for research (e.g., developmental, population, and evolutionary studies), management of wildlife species, and captive breeding programmes for both conservation and poultry. An accurate molecular sexing method applicable across the entire bird radiation is theoretically possible thanks to the long-term stability of their ZZ/ZW sex chromosomes, but current methods are not applicable in a wide range of bird lineages. Here, we developed a novel molecular sexing method based on the comparison of gene copy number variation by quantitative real-time PCR (qPCR) in conserved Z-specific genes (CHRNA6, DDX4, LPAR1, TMEM161B, VPS13A), i.e. genes linked to Z but absent from W chromosomes. We tested the method across three paleognath and 70 neognath species covering the avian phylogeny. In addition, we designed primers for four Z-specific genes (DOCK8, FUT10, PIGG and PSD3) for qPCR-based molecular sexing in three paleognath species. We have demonstrated that the genes DOCK8, FUT10, PIGG and PSD3 can identify sex in paleognath birds and the genes CHRNA6, DDX4, TMEM161B, and VPS13A can reveal sex in neognath birds. The gene LPAR1 can be used to accurately identify sex in both paleognath and neognath species. Along with outlining a novel method of practical importance for molecular sexing in birds, our study also documents in detail the conservation of sex chromosomes across the avian phylogeny.

Sex determination through X-Y heterogamety in Salix nigra

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

With or Without W? Molecular and Cytogenetic Markers are Not Sufficient for Identification of Environmentally-Induced Sex Reversal in the Bearded Dragon

Transitions from environmental sex determination (ESD) to genotypic sex determination (GSD) require an intermediate step of sex reversal, i.e., the production of individuals with a mismatch between the ancestral genotypic and the phenotypic sex. Among amniotes, the sole well-documented transition in this direction was shown in the laboratory in the central bearded dragon, Pogona vitticeps, where very high incubation temperatures led to the production of females with the male-typical (ZZ) genotype. These sex-reversed females then produced offspring whose sex depended on the incubation temperature. Sex-reversed animals identified by molecular and cytogenetic markers were also reported in the field, and their increasing incidence was speculated as a climate warming-driven transition in sex determination. We show that the molecular and cytogenetic markers normally sex-linked in P. vitticeps are also sex-linked in P. henrylawsoni and P. minor, which points to quite ancient sex chromosomes in this lineage. Nevertheless, we demonstrate, based on a crossing experiment with a male bearded dragon who possesses a mismatch between phenotypic sex and genotype, that the used cytogenetic and molecular markers might not be reliable for the identification of sex reversal. Sex reversal should not be considered as the only mechanism causing a mismatch between genetic sex-linked markers and phenotypic sex, which can emerge also by other processes, here most likely by a rare recombination between regions of sex chromosomes which are normally sex-linked. We warn that sex-linked, even apparently for a long evolutionary time, and sex-specific molecular and cytogenetic markers are not a reliable tool for the identification of sex-reversed individuals in a population and that sex reversal has to be verified by other approaches, particularly by observation of the sex ratio of the progeny.

Dynamic evolution of transposable elements, demographic history, and gene content of paleognathous birds

Palaeognathae includes ratite and tinamou species that are important for understanding early avian evolution. Here, we analyzed the whole-genome sequences of 15 paleognathous species to infer their demographic histories, which are presently unknown. We found that most species showed a reduction of population size since the beginning of the last glacial period, except for those species distributed in Australasia and in the far south of South America. Different degrees of contraction and expansion of transposable elements (TE) have shaped the paleognathous genome architecture, with a higher transposon removal rate in tinamous than in ratites. One repeat family, AviRTE, likely underwent horizontal transfer from tropical parasites to the ancestor of little and undulated tinamous about 30 million years ago. Our analysis of gene families identified rapid turnover of immune and reproduction-related genes but found no evidence of gene family changes underlying the convergent evolution of flightlessness among ratites. We also found that mitochondrial genes have experienced a faster evolutionary rate in tinamous than in ratites, with the former also showing more degenerated W chromosomes. This result can be explained by the Hill-Robertson interference affecting genetically linked W chromosomes and mitochondria. Overall, we reconstructed the evolutionary history of the Palaeognathae populations, genes, and TEs. Our findings of co-evolution between mitochondria and W chromosomes highlight the key difference in genome evolution between species with ZW sex chromosomes and those with XY sex chromosomes.

Comparative Mapping of the Macrochromosomes of Eight Avian Species Provides Further Insight into Their Phylogenetic Relationships and Avian Karyotype Evolution

Avian genomes typically consist of ~10 pairs of macro- and ~30 pairs of microchromosomes. While inter-chromosomally, a pattern emerges of very little change (with notable exceptions) throughout evolution, intrachromosomal changes remain relatively poorly studied. To rectify this, here we use a pan-avian universally hybridising set of 74 chicken bacterial artificial chromosome (BAC) probes on the macrochromosomes of eight bird species: common blackbird, Atlantic canary, Eurasian woodcock, helmeted guinea fowl, houbara bustard, mallard duck, and rock dove. A combination of molecular cytogenetic, bioinformatics, and mathematical analyses allowed the building of comparative cytogenetic maps, reconstruction of a putative Neognathae ancestor, and assessment of chromosome rearrangement patterns and phylogenetic relationships in the studied neognath lineages. We observe that, as with our previous studies, chicken appears to have the karyotype most similar to the ancestor; however, previous reports of an increased rate of intrachromosomal change in Passeriformes (songbirds) appear not to be the case in our dataset. The use of this universally hybridizing probe set is applicable not only for the re-tracing of avian karyotype evolution but, potentially, for reconstructing genome assemblies.

Repeat Sequence Mapping Shows Different W Chromosome Evolutionary Pathways in Two Caprimulgiformes Families

Although birds belonging to order Caprimulgiformes show extensive karyotype variation, data concerning their genomic organization is still scarce, as most studies have presented only results obtained from conventional staining analyses. Nevertheless, some interesting findings have been observed, such as the W chromosome of the Common Potoo, Nyctibius griseus (2n = 86), which has the same morphology and size of the Z chromosome, a rare feature in Neognathae birds. Hence, we aimed to investigate the process by which the W chromosome of this species was enlarged. For that, we analyzed comparatively the chromosome organization of the Common Potoo and the Scissor-tailed Nightjar, Hydropsalis torquata (2n = 74), which presents the regular differentiated sex chromosomes, by applying C-banding, G-banding and mapping of repetitive DNAs (microsatellite repeats and 18S rDNA). Our results showed an accumulation of constitutive heterochromatin in the W chromosome of both species. However, 9 out of 11 microsatellite sequences hybridized in the large W chromosome in the Common Potoo, while none of them hybridized in the W chromosome of the Scissor-tailed Nightjar. Therefore, we can conclude that the accumulation of microsatellite sequences, and consequent increase in constitutive heterochromatin, was responsible for the enlargement of the W chromosome in the Common Potoo. Based on these results, we conclude that even though these two species belong to the same order, their W chromosomes have gone through different evolutionary histories, with an extra step of accumulation of repetitive sequences in the Common Potoo.

Molecular evolution of the ATP-binding cassette subfamily G member 2 gene subfamily and its paralogs in birds

Background

ATP-binding cassette (ABC) transporters are involved in the active transportation of various endogenous or exogenous substances. Two ABCG2 gene subfamily members have been identified in birds. A detailed comparative study of the ABCG2 and ABCG2-like genes aid our understanding of their evolutionary history at the molecular level and provide a theoretical reference for studying the specific functions of ABCG2 and ABCG2-like genes in birds.

Results

We first identified 77 ABCG2/ABCG2-like gene sequences in the genomes of 41 birds. Further analysis showed that both the nucleic acid and amino acid sequences of ABCG2 and ABCG2-like genes were highly conserved and exhibited high homology in birds. However, significant differences in the N-terminal structure were found between the ABCG2 and ABCG2-like amino acid sequences. A selective pressure analysis showed that the ABCG2 and ABCG2-like genes were affected by purifying selection during the process of bird evolution.

Conclusions

We believe that multiple members of the ABCG2 gene subfamily exist on chromosome 4 in the ancestors of birds. Over the long course of evolution, only the ABCG2 gene was retained on chromosome 4 in birds. The ABCG2-like gene on chromosome 6 might have originated from chromosome replication or fusion. The structural differences between the N terminus of ABCG2 protein and those of ABCG2-like proteins might lead to functional differences between the corresponding genes.

Cross-Species BAC Mapping Highlights Conservation of Chromosome Synteny across Dragon Lizards (Squamata: Agamidae)

Dragon lizards (Squamata: Agamidae) comprise about 520 species in six subfamilies distributed across Asia, Australasia and Africa. Only five species are known to have sex chromosomes. All of them possess ZZ/ZW sex chromosomes, which are microchromosomes in four species from the subfamily Amphibolurinae, but much larger in Phrynocephalus vlangalii from the subfamily Agaminae. In most previous studies of these sex chromosomes, the focus has been on Australian species from the subfamily Amphibolurinae, but only the sex chromosomes of the Australian central bearded dragon (Pogona vitticeps) are well-characterized cytogenetically. To determine the level of synteny of the sex chromosomes of P. vitticeps across agamid subfamilies, we performed cross-species two-colour FISH using two bacterial artificial chromosome (BAC) clones from the pseudo-autosomal regions of P. vitticeps. We mapped these two BACs across representative species from all six subfamilies as well as two species of chameleons, the sister group to agamids. We found that one of these BAC sequences is conserved in macrochromosomes and the other in microchromosomes across the agamid lineages. However, within the Amphibolurinae, there is evidence of multiple chromosomal rearrangements with one of the BACs mapping to the second-largest chromosome pair and to the microchromosomes in multiple species including the sex chromosomes of P. vitticeps. Intriguingly, no hybridization signal was observed in chameleons for either of these BACs, suggesting a likely agamid origin of these sequences. Our study shows lineage-specific evolution of sequences/syntenic blocks and successive rearrangements and reveals a complex history of sequences leading to their association with important biological processes such as the evolution of sex chromosomes and sex determination.

The sex with the reduced sex chromosome dies earlier: a comparison across the tree of life

Many taxa show substantial differences in lifespan between the sexes. However, these differences are not always in the same direction. In mammals, females tend to live longer than males, while in birds, males tend to live longer than females. One possible explanation for these differences in lifespan is the unguarded X hypothesis, which suggests that the reduced or absent chromosome in the heterogametic sex (e.g. the Y chromosome in mammals and the W chromosome in birds) exposes recessive deleterious mutations on the other sex chromosome. While the unguarded X hypothesis is intuitively appealing, it had never been subject to a broad test. We compiled male and female longevity data for 229 species spanning 99 families, 38 orders and eight classes across the tree of life. Consistent with the unguarded X hypothesis, a meta-analysis showed that the homogametic sex, on average, lives 17.6% longer than the heterogametic sex. Surprisingly, we found substantial differences in lifespan dimorphism between female heterogametic species (in which the homogametic sex lives 7.1% longer) and male heterogametic species (in which the homogametic sex lives 20.9% longer). Our findings demonstrate the importance of considering chromosome morphology in addition to sexual selection and environment as potential drivers of sexual dimorphism, and advance our fundamental understanding of the mechanisms that shape an organism's lifespan.

First Report of Sex Chromosomes in Night Lizards (Scincoidea: Xantusiidae)

Squamate reptiles (lizards, snakes, and amphibians) are an outstanding group for studying sex chromosome evolution—they are old, speciose, geographically widespread, and exhibit myriad sex-determining modes. Yet, the vast majority of squamate species lack heteromorphic sex chromosomes. Cataloging the sex chromosome systems of species lacking easily identifiable, heteromorphic sex chromosomes, therefore, is essential before we are to fully understand the evolution of vertebrate sex chromosomes. Here, we use restriction site-associated DNA sequencing (RADseq) to classify the sex chromosome system of the granite night lizard, Xantusia henshawi. RADseq is an effective alternative to traditional cytogenetic methods for determining a species’ sex chromosome system (i.e., XX/XY or ZZ/ZW), particularly in taxa with non-differentiated sex chromosomes. Although many xantusiid lineages have been karyotyped, none possess heteromorphic sex chromosomes. We identified a ZZ/ZW sex chromosome system in X. henshawi—the first such data for this family. Furthermore, we report that the X. henshawi sex chromosome contains fragments of genes found on Gallus gallus chromosomes 7, 12, and 18 (which are homologous to Anolis carolinensis chromosome 2), the first vertebrate sex chromosomes to utilize this linkage group.

Application of DArT seq derived SNP tags for comparative genome analysis in fishes; An alternative pipeline using sequence data from a non-traditional model species, Macquaria ambigua

Bi-allelic Single Nucleotide Polymorphism (SNP) markers are widely used in population genetic studies. In most studies, sequences either side of the SNPs remain unused, although these sequences contain information beyond that used in population genetic studies. In this study, we show how these sequence tags either side of a single nucleotide polymorphism can be used for comparative genome analysis. We used DArTseq (Diversity Array Technology) derived SNP data for a non-model Australian native freshwater fish, Macquaria ambigua, to identify genes linked to SNP associated sequence tags, and to discover homologies with evolutionarily conserved genes and genomic regions. We concatenated 6,776 SNP sequence tags to create a hypothetical genome (representing 0.1–0.3% of the actual genome), which we used to find sequence homologies with 12 model fish species using the Ensembl genome browser with stringent filtering parameters. We identified sequence homologies for 17 evolutionarily conserved genes (cd9b, plk2b, rhot1b, sh3pxd2aa, si:ch211-148f13.1, si:dkey-166d12.2, zgc:66447, atp8a2, clvs2, lyst, mkln1, mnd1, piga, pik3ca, plagl2, rnf6, sec63) along with an ancestral evolutionarily conserved syntenic block (euteleostomi Block_210). Our analysis also revealed repetitive sequences covering approximately 12% of the hypothetical genome where DNA transposon, LTR and non-LTR retrotransposons were most abundant. A hierarchical pattern of the number of sequence homologies with phylogenetically close species validated the approach for repeatability. This new approach of using SNP associated sequence tags for comparative genome analysis may provide insight into the genome evolution of non-model species where whole genome sequences are unavailable.

Repeated sex chromosome evolution in vertebrates supported by expanded avian sex chromosomes

Sex chromosomes have evolved from the same autosomes multiple times across vertebrates, suggesting that selection for recombination suppression has acted repeatedly and independently on certain genetic backgrounds. Here, we perform comparative genomics of a bird clade (larks and their sister lineage; Alaudidae and Panuridae) where multiple autosome-sex chromosome fusions appear to have formed expanded sex chromosomes. We detected the largest known avian sex chromosome (195.3 Mbp) and show that it originates from fusions between parts of four avian chromosomes: Z, 3, 4A and 5. Within these four chromosomes, we found evidence of five evolutionary strata where recombination had been suppressed at different time points, and show that stratum age explained the divergence rate of Z-W gametologs. Next, we analysed chromosome content and found that chromosome 3 was significantly enriched for genes with predicted sex-related functions. Finally, we demonstrate extensive homology to sex chromosomes in other vertebrate lineages: chromosomes Z, 3, 4A and 5 have independently evolved into sex chromosomes in fish (Z), turtles (Z, 5), lizards (Z, 4A), mammals (Z, 4A) and frogs (Z, 3, 4A, 5). Our results provide insights into and support for repeated evolution of sex chromosomes in vertebrates.

Conserved sex chromosomes and karyotype evolution in monitor lizards (Varanidae)

Despite their long history with the basal split dating back to the Eocene, all species of monitor lizards (family Varanidae) studied so far share the same chromosome number of 2n = 40. However, there are differences in the morphology of the macrochromosome pairs 5-8. Further, sex determination, which revealed ZZ/ZW sex microchromosomes, was studied only in a few varanid species and only with techniques that did not test their homology. The aim of this study was to (i) test if cryptic interchromosomal rearrangements of larger chromosomal blocks occurred during the karyotype evolution of this group, (ii) contribute to the reconstruction of the varanid ancestral karyotype, and (iii) test homology of sex chromosomes among varanids. We investigated these issues by hybridizing flow sorted chromosome paints from Varanus komodoensis to metaphases of nine species of monitor lizards. The results show that differences in the morphology of the chromosome pairs 5-8 can be attributed to intrachromosomal rearrangements, which led to transitions between acrocentric and metacentric chromosomes in both directions. We also documented the first case of spontaneous triploidy among varanids in Varanus albigularis. The triploid individual was fully grown, which demonstrates that polyploidization is compatible with life in this lineage. We found that the W chromosome differs between species in size and heterochromatin content. The varanid Z chromosome is clearly conserved in all the analyzed species. Varanids, in addition to iguanas, caenophidian snakes, and lacertid lizards, are another squamate group with highly conserved sex chromosomes over a long evolutionary time.

Evaluation of a Chicken 600K SNP genotyping array in non-model species of grouse

The use of single nucleotide polymorphism (SNP) arrays to generate large SNP datasets for comparison purposes have recently become an attractive alternative to other genotyping methods. Although most SNP arrays were originally developed for domestic organisms, they can be effectively applied to wild relatives to obtain large panels of SNPs. In this study, we tested the cross-species application of the Affymetrix 600K Chicken SNP array in five species of North American prairie grouse (Centrocercus and Tympanuchus genera). Two individuals were genotyped per species for a total of ten samples. A high proportion (91%) of the total 580 961 SNPs were genotyped in at least one individual (73–76% SNPs genotyped per species). Principal component analysis with autosomal SNPs separated the two genera, but failed to clearly distinguish species within genera. Gene ontology analysis identified a set of genes related to morphogenesis and development (including genes involved in feather development), which may be primarily responsible for large phenotypic differences between Centrocercus and Tympanuchus grouse. Our study provided evidence for successful cross-species application of the chicken SNP array in grouse which diverged ca. 37 mya from the chicken lineage. As far as we are aware, this is the first reported application of a SNP array in non-passerine birds, and it demonstrates the feasibility of using commercial SNP arrays in research on non-model bird species.

Molecular cytogenetic characterization of repetitive sequences comprising centromeric heterochromatin in three Anseriformes species

The highly repetitive DNA sequence of centromeric heterochromatin is an effective molecular cytogenetic marker for investigating genomic compartmentalization between macrochromosomes and microchromosomes in birds. We isolated four repetitive sequence families of centromeric heterochromatin from three Anseriformes species, viz., domestic duck (Anas platyrhynchos, APL), bean goose (Anser fabalis, AFA), and whooper swan (Cygnus cygnus, CCY), and characterized the sequences by molecular cytogenetic approach. The 190-bp APL-HaeIII and 101-bp AFA-HinfI-S sequences were localized in almost all chromosomes of A. platyrhynchos and A. fabalis, respectively. However, the 192-bp AFA-HinfI-L and 290-bp CCY-ApaI sequences were distributed in almost all microchromosomes of A. fabalis and in approximately 10 microchromosomes of C. cygnus, respectively. APL-HaeIII, AFA-HinfI-L, and CCY-ApaI showed partial sequence homology with the chicken nuclear-membrane-associated (CNM) repeat families, which were localized primarily to the centromeric regions of microchromosomes in Galliformes, suggesting that ancestral sequences of the CNM repeat families are observed in the common ancestors of Anseriformes and Galliformes. These results collectively provide the possibility that homogenization of centromeric heterochromatin occurred between microchromosomes in Anseriformes and Galliformes; however, homogenization between macrochromosomes and microchromosomes also occurred in some centromeric repetitive sequences.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes

BACKGROUND

The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions.

RESULTS

We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest.

CONCLUSIONS

The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

Banding cytogenetics of chimeric hybrids Coturnixcoturnix × Coturnixjaponica and comparative analysis with the domestic fowl

The Common quail Coturnixcoturnix Linnaeus, 1758 is a wild migratory bird which is distributed in Eurasia and North Africa, everywhere with an accelerating decline in population size. This species is protected by the Bonn and Berne conventions (1979) and by annex II/1 of the Birds Directive (2009). In Algeria, its breeding took place at the hunting centre in the west of the country. Breeding errors caused uncontrolled crosses between the Common quail and Japanese quail Coturnixjaponica Temminck & Schlegel, 1849. In order to help to preserve the natural genetic heritage of the Common quail and to lift the ambiguity among the populations of quail raised in Algeria, it seemed essential to begin to describe the chromosomes of this species in the country since no cytogenetic study has been reported to date. Fibroblast cultures from embryo and adult animal were initiated. Double synchronization with excess thymidine allowed us to obtain high resolution chromosomes blocked at prometaphase stage. The karyotype and the idiogram in GTG morphological banding (G-bands obtained with trypsin and Giemsa) corresponding to larger chromosomes 1-12 and ZW pair were thus established. The diploid set of chromosomes was estimated as 2N=78. Cytogenetic analysis of expected hybrid animals revealed the presence of a genetic introgression and cellular chimerism. This technique is effective in distinguishing the two quail taxa. Furthermore, the comparative chromosomal analysis of the two quails and domestic chicken Gallusgallusdomesticus Linnaeus, 1758 has been conducted. Differences in morphology and/or GTG band motifs were observed on 1, 2, 4, 7, 8 and W chromosomes. Neocentromere occurrence was suggested for Common quail chromosome 1 and Chicken chromosomes 4 and W. Double pericentric inversion was observed on the Common quail chromosome 2 while pericentric inversion hypothesis was proposed for Chicken chromosome 8. A deletion on the short arm of the Common quail chromosome 7 was also found. These results suggest that Common quail would be a chromosomally intermediate species between Chicken and Japanese quail. The appearance of only a few intrachromosomal rearrangements that occurred during evolution suggests that the organization of the genome is highly conserved between these three galliform species.

A Genetic Map of Ostrich Z Chromosome and the Role of Inversions in Avian Sex Chromosome Evolution

Recombination arrest is a necessary step for the evolution of distinct sex chromosomes. Structural changes, such as inversions, may represent the mechanistic basis for recombination suppression and comparisons of the structural organization of chromosomes as given by chromosome-level assemblies offer the possibility to infer inversions across species at some detail. In birds, deduction of the process of sex chromosome evolution has been hampered by the lack of a validated chromosome-level assembly from a representative of one of the two basal clades of modern birds, Paleognathae. We therefore developed a high-density genetic linkage map of the ostrich Z chromosome and used this to correct an existing assembly, including correction of a large chimeric superscaffold and the order and orientation of other superscaffolds. We identified the pseudoautosomal region as a 52 Mb segment (≈60% of the Z chromosome) where recombination occurred in both sexes. By comparing the order and location of genes on the ostrich Z chromosome with that of six bird species from the other major clade of birds (Neognathae), and of reptilian outgroup species, 25 Z-linked inversions were inferred in the avian lineages. We defined Z chromosome organization in an early avian ancestor and identified inversions spanning the candidate sex-determining DMRT1 gene in this ancestor, which could potentially have triggered the onset of avian sex chromosome evolution. We conclude that avian sex chromosome evolution has been characterized by a complex process of probably both Z-linked and W-linked inversions (and/or other processes). This study illustrates the need for validated chromosome-level assemblies for inference of genome evolution.

Sex Chromosomes of the Iconic Moth Abraxas grossulariata (Lepidoptera, Geometridae) and Its Congener A. sylvata

The magpie moth, Abraxas grossulariata, is an iconic species in which female heterogamety was discovered at the beginning of the 20th century. However, the sex chromosomes of this species have not yet been cytologically identified. We describe the sex chromosomes of A. grossulariata and its congener, A. sylvata. Although these species split only around 9.5 million years ago, and both species have the expected WZ/ZZ chromosomal system of sex determination and their sex chromosomes share the major ribosomal DNA (rDNA) representing the nucleolar organizer region (NOR), we found major differences between their karyotypes, including between their sex chromosomes. The species differ in chromosome number, which is 2n = 56 in A. grossularita and 2n = 58 in A. sylvata. In addition, A. grossularita autosomes exhibit massive autosomal blocks of heterochromatin, which is a very rare phenomenon in Lepidoptera, whereas the autosomes of A. sylvata are completely devoid of distinct heterochromatin. Their W chromosomes differ greatly. Although they are largely composed of female-specific DNA sequences, as shown by comparative genomic hybridization, cross-species W-chromosome painting revealed considerable sequence differences between them. The results suggest a relatively rapid molecular divergence of Abraxas W chromosomes by the independent spreading of female-specific repetitive sequences.

Karyotype Evolution in Birds: From Conventional Staining to Chromosome Painting

In the last few decades, there have been great efforts to reconstruct the phylogeny of Neoaves based mainly on DNA sequencing. Despite the importance of karyotype data in phylogenetic studies, especially with the advent of fluorescence in situ hybridization (FISH) techniques using different types of probes, the use of chromosomal data to clarify phylogenetic proposals is still minimal. Additionally, comparative chromosome painting in birds is restricted to a few orders, while in mammals, for example, virtually all orders have already been analyzed using this method. Most reports are based on comparisons using Gallus gallus probes, and only a small number of species have been analyzed with more informative sets of probes, such as those from Leucopternis albicollis and Gyps fulvus, which show ancestral macrochromosomes rearranged in alternative patterns. Despite this, it is appropriate to review the available cytogenetic information and possible phylogenetic conclusions. In this report, the authors gather both classical and molecular cytogenetic data and describe some interesting and unique characteristics of karyotype evolution in birds.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2^m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2^m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2^m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2^m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

Weird Animals, Sex, and Genome Evolution

Making my career in Australia exposed me to the tyranny of distance, but it gave me opportunities to study our unique native fauna. Distantly related animal species present genetic variation that we can use to explore the most fundamental biological structures and processes. I have compared chromosomes and genomes of kangaroos and platypus, tiger snakes and emus, devils (Tasmanian) and dragons (lizards). I particularly love the challenges posed by sex chromosomes, which, apart from determining sex, provide stunning examples of epigenetic control and break all the evolutionary rules that we currently understand. Here I describe some of those amazing animals and the insights on genome structure, function, and evolution they have afforded us. I also describe my sometimes-random walk in science and the factors and people who influenced my direction. Being a woman in science is still not easy, and I hope others will find encouragement and empathy in my story.

Spot the difference: Solving the puzzle of hidden pictures in the lizard genome for identification of regeneration factors

All living things share some common life processes, such as growth and reproduction, and have the ability to respond to their environment. However, each type of organism has its own specialized way of managing biological events. Genetic sequences determine phenotypic and physiological traits. Based on genetic information, comparative genomics has been used to delineate the differences and similarities between various genomes, and significant progress has been made in understanding regenerative biology by comparing the genomes of a variety of lower animal models of regeneration, such as planaria, zebra fish, and newts. However, the genome of lizards has been relatively ignored until recently, even though lizards have been studied as an excellent amniote model of tissue regeneration. Very recently, whole genome sequences of lizards have been uncovered, and several attempts have been made to find regeneration factors based on genetic information. In this article, recent advances in comparative analysis of the lizard genome are introduced, and their biological implications and putative applications for regenerative medicine and stem cell biology are discussed. [BMB Reports 2016; 49(5): 249-254]

Highly conserved Z and molecularly diverged W chromosomes in the fish genus Triportheus (Characiformes, Triportheidae)

The main objectives of this study were to test: (1) whether the W-chromosome differentiation matches to species' evolutionary divergence (phylogenetic concordance) and (2) whether sex chromosomes share a common ancestor within a congeneric group. The monophyletic genus Triportheus (Characiformes, Triportheidae) was the model group for this study. All species in this genus so far analyzed have ZW sex chromosome system, where the Z is always the largest chromosome of the karyotype, whereas the W chromosome is highly variable ranging from almost homomorphic to highly heteromorphic. We applied conventional and molecular cytogenetic approaches including C-banding, ribosomal DNA mapping, comparative genomic hybridization (CGH) and cross-species whole chromosome painting (WCP) to test our questions. We developed Z- and W-chromosome paints from T. auritus for cross-species WCP and performed CGH in a representative species (T. signatus) to decipher level of homologies and rates of differentiation of W chromosomes. Our study revealed that the ZW sex chromosome system had a common origin, showing highly conserved Z chromosomes and remarkably divergent W chromosomes. Notably, the W chromosomes have evolved to different shapes and sequence contents within ~15-25 Myr of divergence time. Such differentiation highlights a dynamic process of W-chromosome evolution within congeneric species of Triportheus.

Sixteen kiwi (Apteryx spp) transcriptomes provide a wealth of genetic markers and insight into sex chromosome evolution in birds

Background

Kiwi represent the most basal extant avian lineage (paleognaths) and exhibit biological attributes that are unusual or extreme among living birds, such as large egg size, strong olfaction, nocturnality, flightlessness and long lifespan. Despite intense interest in their evolution and their threatened status, genomic resources for kiwi were virtually non-existent until the recent publication of a single genome. Here we present the most comprehensive kiwi transcriptomes to date, obtained via Illumina sequencing of whole blood and de novo assembly of mRNA sequences of eight individuals from each of the two rarest kiwi species, little spotted kiwi (LSK; Apteryx owenii) and rowi (A. rowi).

Results

Sequences obtained were orthologous with a wide diversity of functional genes despite the sequencing of a single tissue type. Individual and composite assemblies contain more than 7900 unique protein coding transcripts in each of LSK and rowi that show strong homology with chicken (Gallus gallus), including those associated with growth, development, disease resistance, reproduction and behavior. The assemblies also contain 66,909 SNPs that distinguish between LSK and rowi, 12,384 SNPs among LSK (associated with 3088 genes), and 29,313 SNPs among rowi (associated with 4953 genes). We found 3084 transcripts differentially expressed between LSK and rowi and 150 transcripts differentially expressed between the sexes. Of the latter, 83 could be mapped to chicken chromosomes with 95% syntenic with chromosome Z.

Conclusions

Our study has simultaneously sequenced multiple species, sexes, and individual kiwi at thousands of genes, and thus represents a significant leap forward in genomic resources available for kiwi. The expression pattern we observed among chromosome Z related genes in kiwi is similar to that observed in ostriches and emu, suggesting a common and ancestral pattern of sex chromosome homomorphy, recombination, and gene dosage among living paleognaths. The transcriptome assemblies described here will provide a rich resource for polymorphic marker development and studies of adaptation of these highly unusual and endangered birds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2714-2) contains supplementary material, which is available to authorized users.

Evolution of vertebrate sex chromosomes and dosage compensation

Differentiated sex chromosomes in mammals and other vertebrates evolved independently but in strikingly similar ways. Vertebrates with differentiated sex chromosomes share the problems of the unequal expression of the genes borne on sex chromosomes, both between the sexes and with respect to autosomes. Dosage compensation of genes on sex chromosomes is surprisingly variable - and can even be absent - in different vertebrate groups. Systems that compensate for different gene dosages include a wide range of global, regional and gene-by-gene processes that differ in their extent and their molecular mechanisms. However, many elements of these control systems are similar across distant phylogenetic divisions and show parallels to other gene silencing systems. These dosage systems cannot be identical by descent but were probably constructed from elements of ancient silencing mechanisms that are ubiquitous among vertebrates and shared throughout eukaryotes.

The X Chromosome of Hemipteran Insects: Conservation, Dosage Compensation and Sex-Biased Expression

Insects of the order Hemiptera (true bugs) use a wide range of mechanisms of sex determination, including genetic sex determination, paternal genome elimination, and haplodiploidy. Genetic sex determination, the prevalent mode, is generally controlled by a pair of XY sex chromosomes or by an XX/X0 system, but different configurations that include additional sex chromosomes are also present. Although this diversity of sex determining systems has been extensively studied at the cytogenetic level, only the X chromosome of the model pea aphid Acyrthosiphon pisum has been analyzed at the genomic level, and little is known about X chromosome biology in the rest of the order.

In this study, we take advantage of published DNA- and RNA-seq data from three additional Hemiptera species to perform a comparative analysis of the gene content and expression of the X chromosome throughout this clade. We find that, despite showing evidence of dosage compensation, the X chromosomes of these species show female-biased expression, and a deficit of male-biased genes, in direct contrast to the pea aphid X. We further detect an excess of shared gene content between these very distant species, suggesting that despite the diversity of sex determining systems, the same chromosomal element is used as the X throughout a large portion of the order.

Intrachromosomal rearrangements in two representatives of the genus Saltator (Thraupidae, Passeriformes) and the occurrence of heteromorphic Z chromosomes

Saltator is a genus within family Thraupidae, the second largest family of Passeriformes, with more than 370 species found exclusively in the New World. Despite this, only a few species have had their karyotypes analyzed, most of them only with conventional staining. The diploid number is close to 80, and chromosome morphology is similar to the usual avian karyotype. Recent studies using cross-species chromosome painting have shown that, although the chromosomal morphology and number are similar to many species of birds, Passeriformes exhibit a complex pattern of paracentric and pericentric inversions in the chromosome homologous to GGA1q in two different suborders, Oscines and Suboscines. Hence, considering the importance and species richness of Thraupidae, this study aims to analyze two species of genus Saltator, the golden-billed saltator (S. aurantiirostris) and the green-winged saltator (S. similis) by means of classical cytogenetics and cross-species chromosome painting using Gallus gallus and Leucopternis albicollis probes, and also 5S and 18S rDNA and telomeric sequences. The results show that the karyotypes of these species are similar to other species of Passeriformes. Interestingly, the Z chromosome appears heteromorphic in S. similis, varying in morphology from acrocentric to metacentric. 5S and 18S probes hybridize to one pair of microchromosomes each, and telomeric sequences produce signals only in the terminal regions of chromosomes. FISH results are very similar to the Passeriformes already analyzed by means of molecular cytogenetics (Turdus species and Elaenia spectabilis). However, the paracentric and pericentric inversions observed in Saltator are different from those detected in these species, an observation that helps to explain the probable sequence of rearrangements. As these rearrangements are found in both suborders of Passeriformes (Oscines and Suboscines), we propose that the fission of GGA1 and inversions in GGA1q have occurred very early after the radiation of this order.

Gene turnover in the avian globin gene families and evolutionary changes in hemoglobin isoform expression

The apparent stasis in the evolution of avian chromosomes suggests that birds may have experienced relatively low rates of gene gain and loss in multigene families. To investigate this possibility and to explore the phenotypic consequences of variation in gene copy number, we examined evolutionary changes in the families of genes that encode the α- and β-type subunits of hemoglobin (Hb), the tetrameric α2β2 protein responsible for blood-O2 transport. A comparative genomic analysis of 52 bird species revealed that the size and membership composition of the α- and β-globin gene families have remained remarkably constant during approximately 100 My of avian evolution. Most interspecific variation in gene content is attributable to multiple independent inactivations of the α(D)-globin gene, which encodes the α-chain subunit of a functionally distinct Hb isoform (HbD) that is expressed in both embryonic and definitive erythrocytes. Due to consistent differences in O2-binding properties between HbD and the major adult-expressed Hb isoform, HbA (which incorporates products of the α(A)-globin gene), recurrent losses of α(D)-globin contribute to among-species variation in blood-O2 affinity. Analysis of HbA/HbD expression levels in the red blood cells of 122 bird species revealed high variability among lineages and strong phylogenetic signal. In comparison with the homologous gene clusters in mammals, the low retention rate for lineage-specific gene duplicates in the avian globin gene clusters suggests that the developmental regulation of Hb synthesis in birds may be more highly conserved, with orthologous genes having similar stage-specific expression profiles and similar functional properties in disparate taxa.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.

Temporal genomic evolution of bird sex chromosomes

BACKGROUND

Sex chromosomes exhibit many unusual patterns in sequence and gene expression relative to autosomes. Birds have evolved a female heterogametic sex system (male ZZ, female ZW), through stepwise suppression of recombination between chrZ and chrW. To address the broad patterns and complex driving forces of Z chromosome evolution, we analyze here 45 newly available bird genomes and four species' transcriptomes, over their course of recombination loss between the sex chromosomes.

RESULTS

We show Z chromosomes in general have a significantly higher substitution rate in introns and synonymous protein-coding sites than autosomes, driven by the male-to-female mutation bias ('male-driven evolution' effect). Our genome-wide estimate reveals that the degree of such a bias ranges from 1.6 to 3.8 among different species. G + C content of third codon positions exhibits the same trend of gradual changes with that of introns, between chrZ and autosomes or regions with increasing ages of becoming Z-linked, therefore codon usage bias in birds is probably driven by the mutational bias. On the other hand, Z chromosomes also evolve significantly faster at nonsynonymous sites relative to autosomes ('fast-Z' evolution). And species with a lower level of intronic heterozygosities tend to evolve even faster on the Z chromosome. Further analysis of fast-evolving genes' enriched functional categories and sex-biased expression patterns support that, fast-Z evolution in birds is mainly driven by genetic drift. Finally, we show in species except for chicken, gene expression becomes more male-biased within Z-linked regions that have became hemizygous in females for a longer time, suggesting a lack of global dosage compensation in birds, and the reported regional dosage compensation in chicken has only evolved very recently.

CONCLUSIONS

In conclusion, we uncover that the sequence and expression patterns of Z chromosome genes covary with their ages of becoming Z-linked. In contrast to the mammalian X chromosomes, such patterns are mainly driven by mutational bias and genetic drift in birds, due to the opposite sex-biased inheritance of Z vs. X.