Human and hamster sera correlate well in identifying antigenic drift among SARS-CoV-2 variants, including JN.1.

Antigenic assessments of SARS-CoV-2 variants inform decisions to update COVID-19 vaccines. Primary infection sera are often used for assessments, but such sera are rare due to population immunity from SARS-CoV-2 infections and COVID-19 vaccinations. Here, we show that neutralization titers and breadth of matched human and hamster pre-Omicron variant primary infection sera correlate well and generate similar antigenic maps. The hamster antigenic map shows modest antigenic drift among XBB sub-lineage variants, with JN.1 and BA.4/BA.5 variants within the XBB cluster, but with fivefold to sixfold antigenic differences between these variants and XBB.1.5. Compared to sera following only ancestral or bivalent COVID-19 vaccinations, or with post-vaccination infections, XBB.1.5 booster sera had the broadest neutralization against XBB sub-lineage variants, although a fivefold titer difference was still observed between JN.1 and XBB.1.5 variants. These findings suggest that antibody coverage of antigenically divergent JN.1 could be improved with a matched vaccine antigen.IMPORTANCEUpdates to COVID-19 vaccine antigens depend on assessing how much vaccine antigens differ antigenically from newer SARS-CoV-2 variants. Human sera from single variant infections are ideal for discriminating antigenic differences among variants, but such primary infection sera are now rare due to high population immunity. It remains unclear whether sera from experimentally infected animals could substitute for human sera for antigenic assessments. This report shows that neutralization titers of variant-matched human and hamster primary infection sera correlate well and recognize variants similarly, indicating that hamster sera can be a proxy for human sera for antigenic assessments. We further show that human sera following an XBB.1.5 booster vaccine broadly neutralized XBB sub-lineage variants but titers were fivefold lower against the more recent JN.1 variant. These findings support updating the current COVID-19 vaccine variant composition and developing a framework for assessing antigenic differences in future variants using hamster primary infection sera.

Clustered cases of human adenovirus types 4, 7, and 14 infections in US Department of Defense Beneficiaries during the 2018-2019 season.

Human adenoviruses (HAdV) are genetically diverse and can infect a number of tissues with severities varied from mild to fatal. HAdV types 3, 4, 7, 11, 14, 21, and 55 were associated with acute respiratory illnesses outbreaks in the United States and in other countries. The risk of outbreaks can be effectively controlled by HAdV vaccination or mitigated by screening and preventive measures. During the influenza season 2018-2019, the DoD Global Respiratory Pathogen Surveillance Program (DoDGRS) received 24 300 respiratory specimens. HAdV samples that produced positive cytopathic effects in viral cultivation were subjected to next-generation sequencing for genome sequence assembly, genome typing, whole genome phylogeny, and sequence comparative analyses. A variety of HAdV types were identified in this study, including HAdV types 1-7, 14, 55, and 56. HAdV types 4, 7, and 14 were found in clustered cases in Colorado, Florida, New York, and South Carolina. Comparative sequence analyses of these isolates revealed the emergence of novel genetic mutations despite the stability of adenovirus genomes. Genomic surveillance of HAdV suggested possible undetected outbreaks and shed light on prevalence, genetic divergence, and viral evolution of HAdV. Continued surveillance will inform risk assessment and countermeasures.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

Transmission of SARS-CoV-2 among recruits in a US Army training environment: a brief report.

BACKGROUND: In 2020, preventive measures were implemented to mitigate the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) among 600-700 recruits arriving weekly at a basic combat training (BCT) facility in the southern United States. Trainees were sorted into companies and platoons (cocoons) at arrival, tested, quarantined for 14 days with daily temperature and respiratory-symptom monitoring and retested before release into larger groups for training where symptomatic testing was conducted. Nonpharmaceutical measures, such as masking, and social distancing, were maintained throughout quarantine and BCT. We assessed for SARS-CoV-2 transmission in the quarantine milieu. METHODS: Nasopharyngeal (NP) swabs were collected at arrival and at the end of quarantine and blood specimens at both timepoints and at the end of BCT. Epidemiological characteristics were analyzed for transmission clusters identified from whole-genome sequencing of NP samples. RESULTS: Among 1403 trainees enrolled from 25 August to 7 October 2020, epidemiological analysis identified three transmission clusters (n = 20 SARS-CoV-2 genomes) during quarantine, which spanned five different cocoons. However, SARS-CoV-2 incidence decreased from 2.7% during quarantine to 1.5% at the end of BCT; prevalence at arrival was 3.3%. CONCLUSIONS: These findings suggest layered SARS-CoV-2 mitigation measures implemented during quarantine minimized the risk of further transmission in BCT.

Genomic surveillance of SARS-CoV-2 in US military compounds in Afghanistan reveals multiple introductions and outbreaks of Alpha and Delta variants.

BACKGROUND: With the emergence and spread of SARS-CoV-2 variants, genomic epidemiology and surveillance have proven invaluable tools for variant tracking. Here, we analyzed SARS-CoV-2 samples collected from personnel located at the US/NATO bases across Afghanistan. RESULTS: Sequencing and phylogenetic analyses revealed at least 16 independent introductions of SARS-CoV-2 into four of these relatively isolated compounds during April and May 2021, including multiple introductions of Alpha and Delta variants. Four of the introductions resulted in sustained spread of the virus within, and in two cases between, the compounds. Three of these outbreaks, one Delta and two Alpha, occurred simultaneously. CONCLUSIONS: Even in rigorously controlled and segregated environments, SARS-CoV-2 introduction and spread may occur frequently.

Origin of a Giant Sex Chromosome.

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

Tbx2a Modulates Switching of RH2 and LWS Opsin Gene Expression.

Sensory systems are tuned by selection to maximize organismal fitness in particular environments. This tuning has implications for intraspecies communication, the maintenance of species boundaries, and speciation. Tuning of color vision largely depends on the sequence of the expressed opsin proteins. To improve tuning of visual sensitivities to shifts in habitat or foraging ecology over the course of development, many organisms change which opsins are expressed. Changes in this developmental sequence (heterochronic shifts) can create differences in visual sensitivity among closely related species. The genetic mechanisms by which these developmental shifts occur are poorly understood. Here, we use quantitative trait locus analyses, genome sequencing, and gene expression studies in African cichlid fishes to identify a role for the transcription factor Tbx2a in driving a switch between long wavelength sensitive (LWS) and Rhodopsin-like (RH2) opsin expression. We identify binding sites for Tbx2a in the LWS promoter and the highly conserved locus control region of RH2 which concurrently promote LWS expression while repressing RH2 expression. We also present evidence that a single change in Tbx2a regulatory sequence has led to a species difference in visual tuning, providing the first mechanistic model for the evolution of rapid switches in sensory tuning. This difference in visual tuning likely has important roles in evolution as it corresponds to differences in diet, microhabitat choice, and male nuptial coloration.

Behavior-dependent cis regulation reveals genes and pathways associated with bower building in cichlid fishes.

Many behaviors are associated with heritable genetic variation [Kendler and Greenspan (2006) Am J Psychiatry 163:1683-1694]. Genetic mapping has revealed genomic regions or, in a few cases, specific genes explaining part of this variation [Bendesky and Bargmann (2011) Nat Rev Gen 12:809-820]. However, the genetic basis of behavioral evolution remains unclear. Here we investigate the evolution of an innate extended phenotype, bower building, among cichlid fishes of Lake Malawi. Males build bowers of two types, pits or castles, to attract females for mating. We performed comparative genome-wide analyses of 20 bower-building species and found that these phenotypes have evolved multiple times with thousands of genetic variants strongly associated with this behavior, suggesting a polygenic architecture. Remarkably, F1 hybrids of a pit-digging and a castle-building species perform sequential construction of first a pit and then a castle bower. Analysis of brain gene expression in these hybrids showed that genes near behavior-associated variants display behavior-dependent allele-specific expression with preferential expression of the pit-digging species allele during pit digging and of the castle-building species allele during castle building. These genes are highly enriched for functions related to neurodevelopment and neural plasticity. Our results suggest that natural behaviors are associated with complex genetic architectures that alter behavior via cis-regulatory differences whose effects on gene expression are specific to the behavior itself.

Movement of transposable elements contributes to cichlid diversity.

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Multiple trans QTL and one cis-regulatory deletion are associated with the differential expression of cone opsins in African cichlids.

BACKGROUND: Dissecting the genetic basis of phenotypic diversity is one of the fundamental goals in evolutionary biology. Despite growing evidence for gene expression divergence being responsible for the evolution of complex traits, knowledge about the proximate genetic causes underlying these traits is still limited. African cichlids have diverse visual systems, with different species expressing different combinations of seven cone opsin genes. Using opsin expression variation in African cichlids as a model for gene expression evolution, this study aims to investigate the genetic architecture of opsin expression divergence in this group. RESULTS: Results from a genome-wide linkage mapping on the F2 progeny of an intergeneric cross, between two species with differential opsin expression show that opsins in Lake Malawi cichlids are controlled by multiple quantitative trait loci (QTLs). Most of these QTLs are located in trans to the opsins except for one cis-QTL for SWS1 on LG17. A closer look at this major QTL revealed the presence of a 691 bp deletion in the promoter of the SWS1 opsin (located 751 bp upstream of the start site) that is associated with a decrease in its expression. Phylogenetic footprinting indicates that the region spanning the deletion harbors a microRNA miR-729 and a conserved non-coding element (CNE) that also occurs in zebrafish and other teleosts. This suggests that the deletion might contain ancestrally preserved regulators that have been tuned for SWS1 gene expression in Lake Malawi. While this deletion is not common, it does occur in several other species within the lake. CONCLUSIONS: Differential expression of cichlid opsins is associated with multiple overlapping QTL, with all but one in trans to the opsins they regulate. The one cis-acting factor is a deletion in the promoter of the SWS1 opsin, suggesting that ancestral polymorphic deletions may contribute to cichlid's visual diversity. In addition to expanding our understanding of the molecular landscape of opsin expression in African cichlids, this study sheds light on the molecular mechanisms underlying phenotypic variation in natural populations.

Transcriptome display during tilapia sex determination and differentiation as revealed by RNA-Seq analysis.

BACKGROUND: The factors determining sex in teleosts are diverse. Great efforts have been made to characterize the underlying genetic network in various species. However, only seven master sex-determining genes have been identified in teleosts. While the function of a few genes involved in sex determination and differentiation has been studied, we are far from fully understanding how genes interact to coordinate in this process. RESULTS: To enable systematic insights into fish sexual differentiation, we generated a dynamic co-expression network from tilapia gonadal transcriptomes at 5, 20, 30, 40, 90, and 180 dah (days after hatching), plus 45 and 90 dat (days after treatment) and linked gene expression profiles to both development and sexual differentiation. Transcriptomic profiles of female and male gonads at 5 and 20 dah exhibited high similarities except for a small number of genes that were involved in sex determination, while drastic changes were observed from 90 to 180 dah, with a group of differently expressed genes which were involved in gonadal differentiation and gametogenesis. Weighted gene correlation network analysis identified changes in the expression of Borealin, Gtsf1, tesk1, Zar1, Cdn15, and Rpl that were correlated with the expression of genes previously known to be involved in sex differentiation, such as Foxl2, Cyp19a1a, Gsdf, Dmrt1, and Amh. CONCLUSIONS: Global gonadal gene expression kinetics during sex determination and differentiation have been extensively profiled in tilapia. These findings provide insights into the genetic framework underlying sex determination and sexual differentiation, and expand our current understanding of developmental pathways during teleost sex determination.

Novel Sex Chromosomes in 3 Cichlid Fishes from Lake Tanganyika.

African cichlids are well known for their adaptive radiations, but it is now apparent that they also harbor an extraordinary diversity of sex chromosome systems. In this study, we sequenced pools of males and females from species in 3 different genera of cichlids from Lake Tanganyika. We then searched for regions that were differentiated following the patterns expected for sex chromosomes. We report 2 novel sex chromosomes systems, an XY system on LG19 in Tropheus sp. "black" and a ZW system on LG7 in Hemibates stenosoma. We also identify a ZW system on LG5 in Cyprichromis leptosoma that may be convergent with a system previously described in Lake Malawi cichlids. Our data also identify candidate single nucleotide polymorphisms for the blue/yellow tail color polymorphism observed among male C. leptosoma.

A high quality assembly of the Nile Tilapia (Oreochromis niloticus) genome reveals the structure of two sex determination regions.

BACKGROUND: Tilapias are the second most farmed fishes in the world and a sustainable source of food. Like many other fish, tilapias are sexually dimorphic and sex is a commercially important trait in these fish. In this study, we developed a significantly improved assembly of the tilapia genome using the latest genome sequencing methods and show how it improves the characterization of two sex determination regions in two tilapia species. RESULTS: A homozygous clonal XX female Nile tilapia (Oreochromis niloticus) was sequenced to 44X coverage using Pacific Biosciences (PacBio) SMRT sequencing. Dozens of candidate de novo assemblies were generated and an optimal assembly (contig NG50 of 3.3Mbp) was selected using principal component analysis of likelihood scores calculated from several paired-end sequencing libraries. Comparison of the new assembly to the previous O. niloticus genome assembly reveals that recently duplicated portions of the genome are now well represented. The overall number of genes in the new assembly increased by 27.3%, including a 67% increase in pseudogenes. The new tilapia genome assembly correctly represents two recent vasa gene duplication events that have been verified with BAC sequencing. At total of 146Mbp of additional transposable element sequence are now assembled, a large proportion of which are recent insertions. Large centromeric satellite repeats are assembled and annotated in cichlid fish for the first time. Finally, the new assembly identifies the long-range structure of both a ~9Mbp XY sex determination region on LG1 in O. niloticus, and a ~50Mbp WZ sex determination region on LG3 in the related species O. aureus. CONCLUSIONS: This study highlights the use of long read sequencing to correctly assemble recent duplications and to characterize repeat-filled regions of the genome. The study serves as an example of the need for high quality genome assemblies and provides a framework for identifying sex determining genes in tilapia and related fish species.

Dynamic Sequence Evolution of a Sex-Associated B Chromosome in Lake Malawi Cichlid Fish.

B chromosomes are extra chromosomes found in many species of plants, animals, and fungi. B chromosomes often manipulate common cellular processes to increase their frequency, sometimes to the detriment of organismal fitness. Here, we characterize B chromosomes in several species of Lake Malawi cichlid fish. Whole genome sequencing of Metriaclima zebra "Boadzulu" individuals revealed blocks of sequence with unusually high sequence coverage, indicative of increased copy number of those sequences. These regions of high sequence coverage were found only in females. SNPs unique to the high copy number sequences permitted the design of specific amplification primers. These primers amplified fragments only in Metriaclima lombardoi individuals that carried a cytologically identified B chromosome (B-carriers), indicating these extra copies are located on the B chromosome. These same primers were used to identify B-carrying individuals in additional species from Lake Malawi. Across 7 species, a total of 43 B-carriers were identified among 323 females. B-carriers were exclusively female; no B chromosomes were observed in the 317 males surveyed from these species. Quantitative analysis of the copy number variation of B-specific sequence blocks suggests that B-carriers possess a single B chromosome, consistent with previous karyotyping of M. lombardoi A single B chromosome in B-carriers is consistent with 2 potential drive mechanisms: one involving nondisjunction and preferential segregation in a mitotic division prior to the germ-line, and the other involving preferential segregation during meiosis I.

Comparative analysis of a sex chromosome from the blackchin tilapia, Sarotherodon melanotheron.

BACKGROUND: Inversions and other structural polymorphisms often reduce the rate of recombination between sex chromosomes, making it impossible to fine map sex-determination loci using traditional genetic mapping techniques. Here we compare distantly related species of tilapia that each segregate an XY system of sex-determination on linkage group 1. We use whole genome sequencing to identify shared sex-patterned polymorphisms, which are candidates for the ancestral sex-determination mutation. RESULTS: We found that Sarotherodon melanotheron segregates an XY system on LG1 in the same region identified in Oreochromis niloticus. Both species have higher densities of sex-patterned SNPs, as well as elevated number of ancestral copy number variants in this region when compared to the rest of the genome, but the pattern of differentiation along LG1 differs between species. The number of sex-patterned SNPs shared by the two species is small, but larger than expected by chance, suggesting that a novel Y-chromosome arose just before the divergence of the two species. We identified a shared sex-patterned SNP that alters a Gata4 binding site near Wilms tumor protein that might be responsible for sex-determination. CONCLUSIONS: Shared sex-patterned SNPs, insertions and deletions suggest an ancestral sex-determination system that is common to both S. melanotheron and O. niloticus. Functional analyses are needed to evaluate shared SNPs near candidate genes that might play a role in sex-determination of these species. Interspecific variation in the sex chromosomes of tilapia species provides an excellent model system for understanding the evolution of vertebrate sex chromosomes.

Integrated analysis of miRNA and mRNA expression profiles in tilapia gonads at an early stage of sex differentiation.

BACKGROUND: MicroRNAs (miRNAs) represent a second regulatory network that has important effects on gene expression and protein translation during biological process. However, the possible role of miRNAs in the early stages of fish sex differentiation is not well understood. In this study, we carried an integrated analysis of miRNA and mRNA expression profiles to explore their possibly regulatory patterns at the critical stage of sex differentiation in tilapia. RESULTS: We identified 279 pre-miRNA genes in tilapia genome, which were highly conserved in other fish species. Based on small RNA library sequencing, we identified 635 mature miRNAs in tilapia gonads, in which 62 and 49 miRNAs showed higher expression in XX and XY gonads, respectively. The predicted targets of these sex-biased miRNAs (e.g., miR-9, miR-21, miR-30a, miR-96, miR-200b, miR-212 and miR-7977) included genes encoding key enzymes in steroidogenic pathways (Cyp11a1, Hsd3b, Cyp19a1a, Hsd11b) and key molecules involved in vertebrate sex differentiation (Foxl2, Amh, Star1, Sf1, Dmrt1, and Gsdf). These genes also showed sex-biased expression in tilapia gonads at 5 dah. Some miRNAs (e.g., miR-96 and miR-737) targeted multiple genes involved in steroid synthesis, suggesting a complex miRNA regulatory network during early sex differentiation in this fish. CONCLUSIONS: The sequence and expression patterns of most miRNAs in tilapia are conserved in fishes, indicating the basic functions of vertebrate miRNAs might share a common evolutionary origin. This comprehensive analysis of miRNA and mRNA at the early stage of molecular sex differentiation in tilapia XX and XY gonads lead to the discovery of differentially expressed miRNAs and their putative targets, which will facilitate studies of the regulatory network of molecular sex determination and differentiation in fishes.

An improved genome reference for the African cichlid, Metriaclima zebra.

BACKGROUND: Problems associated with using draft genome assemblies are well documented and have become more pronounced with the use of short read data for de novo genome assembly. We set out to improve the draft genome assembly of the African cichlid fish, Metriaclima zebra, using a set of Pacific Biosciences SMRT sequencing reads corresponding to 16.5× coverage of the genome. Here we characterize the improvements that these long reads allowed us to make to the state-of-the-art draft genome previously assembled from short read data. RESULTS: Our new assembly closed 68% of the existing gaps and added 90.6 Mbp of new non-gap sequence to the existing draft assembly of M. zebra. Comparison of the new assembly to the sequence of several bacterial artificial chromosome clones confirmed the accuracy of the new assembly. The closure of sequence gaps revealed thousands of new exons, allowing significant improvement in gene models. We corrected one known misassembly, and identified and fixed other likely misassemblies. 63.5 Mbp (70%) of the new sequence was classified as repetitive and the new sequence allowed for the assembly of many more transposable elements. CONCLUSIONS: Our improvements to the M. zebra draft genome suggest that a reasonable investment in long reads could greatly improve many comparable vertebrate draft genome assemblies.

Interspecific variation in Rx1 expression controls opsin expression and causes visual system diversity in African cichlid fishes.

The mechanisms underlying natural phenotypic diversity are key to understanding evolution and speciation. Cichlid fishes are among the most speciose vertebrates and an ideal model for identifying genes controlling species differences. Cichlids have diverse visual sensitivities that result from species expressing subsets of seven cichlid cone opsin genes. We previously identified a quantitative trait locus (QTL) that tunes visual sensitivity by varying SWS2A (short wavelength sensitive 2A) opsin expression in a genetic cross between two Lake Malawi cichlid species. Here, we identify Rx1 (retinal and anterior neural fold homeobox) as the causative gene for the QTL using fine mapping and RNAseq in retinal transcriptomes. Rx1 is differentially expressed between the parental species and correlated with SWS2A expression in the F2 progeny. Expression of Rx1 and SWS2A is also correlated in a panel of 16 Lake Malawi cichlid species. Association mapping in this panel identified a 413-bp deletion located 2.5-kb upstream of the Rx1 translation start site that is correlated with decreased Rx1 expression. This deletion explains 62% of the variance in SWS2A expression across 53 cichlid species in 29 genera. The deletion occurs in both the sand and rock-dwelling cichlid clades, suggesting that it is an ancestral polymorphism. Our finding supports the hypothesis that mixing and matching of ancestral polymorphisms can explain the diversity of present day cichlid phenotypes.

Origin and evolution of B chromosomes in the cichlid fish Astatotilapia latifasciata based on integrated genomic analyses.

Approximately 15% of eukaryotes contain supernumerary B chromosomes. When present, B chromosomes frequently represent as much as 5% of the genome. Despite thousands of reports describing the distribution of supernumeraries in various taxa, a comprehensive theory for the origin, maintenance, and evolution of B chromosomes has not emerged. Here, we sequence the complete genomes of individual cichlid fish (Astatotilapia latifasciata) with and without B chromosomes, as well as microdissected B chromosomes, to identify DNA sequences on the B. B sequences were further analyzed through quantitative polymerase chain reaction and in situ hybridization. We find that the B chromosome contains thousands of sequences duplicated from essentially every chromosome in the ancestral karyotype. Although most genes on the B chromosome are fragmented, a few are largely intact, and we detect evidence that at least three of them are transcriptionally active. We propose a model in which the B chromosome originated early in the evolutionary history of Lake Victoria cichlids from a small fragment of one autosome. DNA sequences originating from several autosomes, including protein-coding genes and transposable elements, subsequently inserted into this proto-B. We propose that intact B chromosome genes involved with microtubule organization, kinetochore structure, recombination and progression through the cell cycle may play a role in driving the transmission of the B chromosome. Furthermore, our work suggests that karyotyping is an essential step prior to genome sequencing to avoid problems in genome assembly and analytical biases created by the presence of high copy number sequences on the B chromosome.