Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkii bouvieri) and rainbow trout (O. mykiss)

Making sense of chromosome polymorphisms in two leptysmine grasshoppers

The touchstone of the 'New Synthesis' was population cytogenetics -rather than genetics - due to the abundant polymorphic inversions in the genus Drosophila. Grasshoppers were not a material of choice because of their conservative karyotypes. However, nowadays seven species of Acrididae were described for polymorphic centric fusions, five of them in South-America. Leptysma argentina and the likely biocontrol of water-hyacinth Cornops aquaticum are semiaquatic Leptysminae (Acrididae: Orthoptera), polymorphic for centric fusions, supernumerary segments and a B-chromosome. We sought to demonstrate the operation of natural selection on them, by detecting: (I) latitudinal clines; (II) regression on environmental variables; (III) deviation from null models, such as linkage equilibrium; (IV) seasonal variation; (V) comparison between age classes and (VI) selection component analyses. All of them were confirmed in L. argentina, just (I) and (II) in C. aquaticum. Furthermore, the relationship between karyotype, phenotype and recombination was confirmed in both species. Karyotype-phenotype relationship may be due to the body enlargement the fusions are associated with, along with a latitudinal transition in voltinism. Karyotype-related recombination reduction in both species may help explain all fusion clines, although there is probably more than one factor at work. No effects were noticed for a supernumerary segment in L. argentina, but it is ubiquitous and certainly non-neutral. C. aquaticum is poised for introduction in South-Africa as a biocontrol of water-hyacinths; the recent discovery of four more segment polymorphisms may imply more chromosomal markers to make sense of its genetic system.

Multiple karyotype differences between populations of the Hoplias malabaricus (Teleostei; Characiformes), a species complex in the gray area of the speciation process

Neotropical fishes exhibit remarkable karyotype diversity, whose evolution is poorly understood. Here, we studied genetic differences in 60 individuals, from 11 localities of one species, the wolf fish Hoplias malabaricus, from populations that include six different "karyomorphs". These differ in Y-X chromosome differentiation, and, in several cases, by fusions with autosomes that have resulted in multiple sex chromosomes. Other differences are also observed in diploid chromosome numbers and morphologies. In an attempt to start understanding how this diversity was generated, we analyzed within- and between-population differences in a genome-wide sequence data set. We detect clear genotype differences between karyomorphs. Even in sympatry, samples with different karyomorphs differ more in sequence than samples from allopatric populations of the same karyomorph, suggesting that they represent populations that are to some degree reproductively isolated. However, sequence divergence between populations with different karyomorphs is remarkably low, suggesting that chromosome rearrangements may have evolved during a brief evolutionary time. We suggest that the karyotypic differences probably evolved in allopatry, in small populations that would have allowed rapid fixation of rearrangements, and that they became sympatric after their differentiation. Further studies are needed to test whether the karyotype differences contribute to reproductive isolation detected between some H. malabaricus karyomorphs.

Transmission ratio distortion detection by neutral genetic markers in the Pura Raza Española horse breed

Transmission Ratio Distortion (TRD) is a genetic phenomenon widely demonstrated in several livestock species, but barely in equine species. The TRD occurs when certain genotypes are over- or under-represented in the offspring of a particular mating and can be caused by a variety of factors during gamete formation or during embryonic development. For this study, 126 394 trios consisting of a stallion, mare, and offspring were genotyped using a panel of 17 neutral microsatellite markers recommended by the International Society for Animal Genetics for paternity tests and individual identification. The number of alleles available for each marker ranges from 13 to 18, been 268 the total number of alleles investigated. The TRDscan v.2.0 software was used with the biallelic procedure to identify regions with distorted segregation ratios. After completing the analysis, a total of 12 alleles (out of 11 microsatellites) were identified with decisive evidence for genotypic TRD; 3 and 9 with additive and heterosis patterns, respectively. In addition, 19 alleles (out of 10 microsatellites) were identified displaying allelic TRD. Among them, 14 and 5 were parent-unspecific and stallion-mare-specific TRD. Out of the TRD regions, 24 genes were identified and annotated, predominantly associated with cholesterol metabolism and homeostasis. These genes are often linked to non-specific symptoms like impaired fertility, stunted growth, and compromised overall health. The results suggest a significant impact on the inheritance of certain genetic traits in horses. Further analysis and validation are needed to better understand the TRD impact before the potential implementation in the horse breeding programme strategies.

The Genetic Differentiation of Pyrrhulina (Teleostei, Characiformes) Species is Likely Influenced by Both Geographical Distribution and Chromosomal Rearrangements

Allopatry is generally considered to be one of the main contributors to the remarkable Neotropical biodiversity. However, the role of chromosomal rearrangements including neo-sex chromosomes for genetic diversity is still poorly investigated and understood. Here, we assess the genetic divergence in five Pyrrhulina species using population genomics and combined the results with previously obtained cytogenetic data, highlighting that molecular genetic diversity is consistent with their chromosomal features. The results of a principal coordinate analysis (PCoA) indicated a clear difference among all species while showing a closer relationship of the ones located in the same geographical region. This was also observed in genetic structure analyses that only grouped P. australis and P. marilynae, which were also recovered as sister species in a species tree analysis. We observed a contradictory result for the relationships among the three species from the Amazon basin, as the phylogenetic tree suggested P. obermulleri and P. semifasciata as sister species, while the PCoA showed a high genetic difference between P. semifasciata and all other species. These results suggest a potential role of sex-related chromosomal rearrangements as reproductive barriers between these species.

Construction of Genetic Linkage Maps From a Hybrid Family of Large Yellow Croaker (Larimichthys crocea)

Consensus and sex-specific genetic linkage maps for large yellow croaker (Larimichthys crocea) were constructed using samples from an F₁ family produced by crossing a Daiqu female and a Mindong male. A total of 20,147 single nucleotide polymorphisms (SNPs) by restriction site associated DNA sequencing were assigned to 24 linkage groups (LGs). The total length of the consensus map was 1757.4 centimorgan (cM) with an average marker interval of 0.09 cM. The total length of female and male linkage map was 1533.1 cM and 1279.2 cM, respectively. The average female-to-male map length ratio was 1.2 ± 0.23. Collapsed markers in the genetic maps were re-ordered according to their relative positions in the ASM435267v1 genome assembly to produce integrated genetic linkage maps with 9885 SNPs distributed across the 24 LGs. The recombination pattern of most LGs showed sigmoidal patterns of recombination, with higher recombination in the middle and suppressed recombination at both ends, which corresponds with the presence of sub-telocentric and acrocentric chromosomes in the species. The average recombination rate in the integrated female and male maps was respectively 3.55 cM/Mb and 3.05 cM/Mb. In most LGs, higher recombination rates were found in the integrated female map, compared to the male map, except in LG12, LG16, LG21, LG22, and LG24. Recombination rate profiles within each LG differed between the male and the female, with distinct regions indicating potential recombination hotspots. Separate quantitative trait loci (QTL) and association analyses for growth related traits in 6 months fish were performed, however, no significant QTL was detected. The study indicates that there may be genetic differences between the two strains, which may have implications for the application of DNA-information in the further breeding schemes.

Chromosome reciprocal translocations have accompanied subspecies evolution in bananas

Chromosome rearrangements and the way that they impact genetic differentiation and speciation have long raised questions from evolutionary biologists. They are also a major concern for breeders because of their bearing on chromosome recombination. Banana is a major crop that derives from inter(sub)specific hybridizations between various once geographically isolated Musa species and subspecies. We sequenced 155 accessions, including banana cultivars and representatives of Musa diversity, and genotyped-by-sequencing 1059 individuals from 11 progenies. We precisely characterized six large reciprocal translocations and showed that they emerged in different (sub)species of Musa acuminata, the main contributor to currently cultivated bananas. Most diploid and triploid cultivars analyzed were structurally heterozygous for 1 to 4 M. acuminata translocations, highlighting their complex origin. We showed that all translocations induced a recombination reduction of variable intensity and extent depending on the translocations, involving only the breakpoint regions, a chromosome arm, or an entire chromosome. The translocated chromosomes were found preferentially transmitted in many cases. We explore and discuss the possible mechanisms involved in this preferential transmission and its impact on translocation colonization.

Newly discovered cichlid fish biodiversity threatened by hybridization with non-native species

Invasive freshwater fishes are known to readily hybridize with indigenous congeneric species, driving loss of unique and irreplaceable genetic resources. Here we reveal that newly discovered (2013-2016) evolutionarily significant populations of Korogwe tilapia (Oreochromis korogwe) from southern Tanzania are threatened by hybridization with the larger invasive Nile tilapia (Oreochromis niloticus). We use a combination of morphology, microsatellite allele frequencies and whole genome sequences to show that O. korogwe from southern lakes (Nambawala, Rutamba and Mitupa) are distinct from geographically disjunct populations in northern Tanzania (Zigi River and Mlingano Dam). We also provide genetic evidence of O. korogwe × niloticus hybrids in three southern lakes and demonstrate heterogeneity in the extent of admixture across the genome. Finally, using the least admixed genomic regions we estimate that the northern and southern O. korogwe populations most plausibly diverged ~140,000 years ago, suggesting that the geographical separation of the northern and southern groups is not a result of a recent translocation, and instead these populations represent independent evolutionarily significant units. We conclude that these newly discovered and phenotypically unique cichlid populations are already threatened by hybridization with an invasive species, and propose that these irreplaceable genetic resources would benefit from conservation interventions.

Genomic Chaos Begets Psychiatric Disorder

The processes that created the primordial genome are inextricably linked to current day vulnerability to developing a psychiatric disorder as summarized in this review article. Chaos and dynamic forces including duplication, transposition, and recombination generated the protogenome. To survive early stages of genome evolution, self-organization emerged to curb chaos. Eventually, the human genome evolved through a delicate balance of chaos/instability and organization/stability. However, recombination coldspots, silencing of transposable elements, and other measures to limit chaos also led to retention of variants that increase risk for disease. Moreover, ongoing dynamics in the genome creates various new mutations that determine liability for psychiatric disorders. Homologous recombination, long-range gene regulation, and gene interactions were all guided by spooky action-at-a-distance, which increased variability in the system. A probabilistic system of life was required to deal with a changing environment. This ensured the generation of outliers in the population, which enhanced the probability that some members would survive unfavorable environmental impacts. Some of the outliers produced through this process in man are ill suited to cope with the complex demands of modern life. Genomic chaos and mental distress from the psychological challenges of modern living will inevitably converge to produce psychiatric disorders in man.

Identification of transcription factors controlling floral morphology in wild Petunia species with contrasting pollination syndromes

Adaptation to different pollinators is an important driver of speciation in the angiosperms. Genetic approaches such as QTL mapping have been successfully used to identify the underlying speciation genes. However, these methods are often limited by widespread suppression of recombination due to divergence between species. While the mutations that caused the interspecific differences in floral color and scent have been elucidated in a variety of plant genera, the genes that are responsible for morphological differences remain mostly unknown. Differences in floral organ length determine the pollination efficiency of hawkmoths and hummingbirds, and therefore the genes that control these differences are potential speciation genes. Identifying such genes is challenging, especially in non-model species and when studying complex traits for which little prior genetic and biochemical knowledge is available. Here we combine transcriptomics with detailed growth analysis to identify candidate transcription factors underlying interspecific variation in the styles of Petunia flowers. Starting from a set of 2284 genes, stepwise filtering for expression in styles, differential expression between species, correlation with growth-related traits, allele-specific expression in interspecific hybrids, and/or high-impact polymorphisms resulted in a set of 43 candidate speciation genes. Validation by virus-induced gene silencing identified two MYB transcription factors, EOBI and EOBII, that were previously shown to regulate floral scent emission, a trait associated with pollination by hawkmoths.

Discriminating between allele- and genotype-specific transmission ratio distortion

Transmission ratio distortion (TRD) is defined as the observed deviation from the expected Mendelian inheritance of alleles from heterozygous parents. This phenomenon is attributed to various biological mechanisms acting on germ cells, embryos or fetuses, or even in early postnatal life. Current statistical approaches typically use two independent parametrizations assuming that TRD relies on allele- or genotype-related mechanisms, although they have never been tested and compared. This study compared allele- and genotype-related TRD models on simulated datasets with 1000 genotyped offspring and real data from 168 sire-dam-offspring beef cattle trios. The analysis of simulated datasets favored the true model of analysis in most cases (>93%), and a low percentage of missidentification occurred under (almost) null dominance (genotype-related model) or similar and moderate-to-low sire- and dam-specific TRD parameters (allele-related model). Moreover, the correlation between simulated and predicted distortion parameters was high (>0.97) under the true model. The comparison of allele- and genotype-related TRD models is an appealing tool to infer the biological source of TRD (i.e. haploid vs. diploid cells) when screening the whole genome. The analysis of beef cattle data corroborated a TRD region previously reported in chromosome 4, although discarding allele-related mechanisms and favoring the genotype-related model as the more reliable one. The results of this study highlight the relevance of implementing and comparing different parametrizations to capture all kinds of TRD, and to compare them using appropriate statistical methods.

Eukaryote hybrid genomes

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve.

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Later-generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Vive la résistance: genome-wide selection against introduced alleles in invasive hybrid zones

Evolutionary and ecological consequences of hybridization between native and invasive species are notoriously complicated because patterns of selection acting on non-native alleles can vary throughout the genome and across environments. Rapid advances in genomics now make it feasible to assess locus-specific and genome-wide patterns of natural selection acting on invasive introgression within and among natural populations occupying diverse environments. We quantified genome-wide patterns of admixture across multiple independent hybrid zones of native westslope cutthroat trout and invasive rainbow trout, the world's most widely introduced fish, by genotyping 339 individuals from 21 populations using 9380 species-diagnostic loci. A significantly greater proportion of the genome appeared to be under selection favouring native cutthroat trout (rather than rainbow trout), and this pattern was pervasive across the genome (detected on most chromosomes). Furthermore, selection against invasive alleles was consistent across populations and environments, even in those where rainbow trout were predicted to have a selective advantage (warm environments). These data corroborate field studies showing that hybrids between these species have lower fitness than the native taxa, and show that these fitness differences are due to selection favouring many native genes distributed widely throughout the genome.

Evolution of the Banana Genome (Musa acuminata) Is Impacted by Large Chromosomal Translocations

Most banana cultivars are triploid seedless parthenocarpic clones derived from hybridization between Musa acuminata subspecies and sometimes M. balbisiana. M. acuminata subspecies were suggested to differ by a few large chromosomal rearrangements based on chromosome pairing configurations in intersubspecies hybrids. We searched for large chromosomal rearrangements in a seedy M. acuminata ssp. malaccensis banana accession through mate-pair sequencing, BAC-FISH, targeted PCR and marker (DArTseq) segregation in its progeny. We identified a heterozygous reciprocal translocation involving two distal 3 and 10 Mb segments from chromosomes 01 and 04, respectively, and showed that it generated high segregation distortion, reduced recombination and linkage between chromosomes 01 and 04 in its progeny. The two chromosome structures were found to be mutually exclusive in gametes and the rearranged structure was preferentially transmitted to the progeny. The rearranged chromosome structure was frequently found in triploid cultivars but present only in wild malaccensis ssp. accessions, thus suggesting that this rearrangement occurred in M. acuminata ssp. malaccensis. We propose a mechanism for the spread of this rearrangement in Musa diversity and suggest that this rearrangement could have played a role in the emergence of triploid cultivars.

Linkage Map of Lissotriton Newts Provides Insight into the Genetic Basis of Reproductive Isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts.

Night Monkey Hybrids Exhibit De Novo Genomic and Karyotypic Alterations: The First Such Case in Primates

Using molecular chromosomal analyses, we discovered night monkey hybrids produced in captivity from matings between a female Aotus azarae boliviensis (2n = 50) and a male Aotus lemurinus griseimembra (2n = 53). The parents produced seven offspring in total, including one male and six females-a pattern consistent with Haldane's rule. Chromosomal studies were conducted on four of the hybrid offspring. Two of them showed relatively "simple" mixture karyotypes, including different chromosome numbers (2n = 51, 52), which were formed because of a heteromorphic autosome pair in the father (n = 26, 27). The other two hybrid monkeys exhibited de novo genomic and karyotypic alterations. Detailed analysis of the alterations revealed that one individual carried a mixture karyotype of the two parental species and an X chromosome trisomy (53,XXX). The second individual displayed trisomy of chromosome 18 (52,XX,+18) and a reciprocal translocation between autosomes 21 and 23 (52,XX,+18,t(21;23)). Interestingly, the second monkey exhibited mosaicism among blood cells (mos52,XX,+18[87]/52,XX,+18,t(21;23)[85]), but only a single karyotype (52,XX,+18) in skin fibroblast cells. The X- and 18-trisomies were derived from a doubling of the mother's chromosomes in early embryonic cell division, and the reciprocal translocation likely developed in the bone marrow of the offspring, considering that it was observed only in blood cells. Such occurrence of trisomies in hybrid individuals is a unique finding in placental mammals.

A SNP Based Linkage Map of the Arctic Charr (Salvelinus alpinus) Genome Provides Insights into the Diploidization Process After Whole Genome Duplication

Diploidization, which follows whole genome duplication events, does not occur evenly across the genome. In salmonid fishes, certain pairs of homeologous chromosomes preserve tetraploid loci in higher frequencies toward the telomeres due to residual tetrasomic inheritance. Research suggests this occurs only in homeologous pairs where one chromosome arm has undergone a fusion event. We present a linkage map for Arctic charr (Salvelinus alpinus), a salmonid species with relatively fewer chromosome fusions. Genotype by sequencing identified 19,418 SNPs, and a linkage map consisting of 4508 markers was constructed from a subset of high quality SNPs and microsatellite markers that were used to anchor the new map to previous versions. Both male- and female-specific linkage maps contained the expected number of 39 linkage groups. The chromosome type associated with each linkage group was determined, and 10 stable metacentric chromosomes were identified, along with a chromosome polymorphism involving the sex chromosome AC04. Two instances of a weak form of pseudolinkage were detected in the telomeric regions of homeologous chromosome arms in both female and male linkage maps. Chromosome arm homologies within the Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss) genomes were determined. Paralogous sequence variants (PSVs) were identified, and their comparative BLASTn hit locations showed that duplicate markers exist in higher numbers on seven pairs of homeologous arms, previously identified as preserving tetrasomy in salmonid species. Homeologous arm pairs where neither arm has been part of a fusion event in Arctic charr had fewer PSVs, suggesting faster diploidization rates in these regions.

Salmonid Chromosome Evolution as Revealed by a Novel Method for Comparing RADseq Linkage Maps

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ∼65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esox lucius for identifying post-duplicated homeologs. To facilitate this process, we developed MapComp to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MapComp, available at: https://github.com/enormandeau/mapcomp/

Bayesian analysis of parent-specific transmission ratio distortion in seven Spanish beef cattle breeds

Transmission ratio distortion (TRD) is the departure from the expected Mendelian ratio in offspring, a poorly investigated biological phenomenon in livestock species. Given the current availability of specific parametric methods for the analysis of segregation data, this study focused on the screening of TRD in 602 402 single nucleotide polymorphisms covering all autosomal chromosomes in seven Spanish beef cattle breeds. On average, 0.13% (n = 786) and 0.01% (n = 29) of genetic markers evidenced sire- or dam-specific TRD respectively. There were no single nucleotide polymorphisms accounting for both sire- and dam-specific TRD at the same time, and only one marker (rs43147474) accounted for (sire-specific) TRD in all seven breeds. It must be noted that rs43147474 is located in the fourth intronic region of the GTP-binding protein 10 gene, and this locus has been previously linked to the maintenance of mitochondria and nucleolar architectures. Alternatively, other candidate genes surround this hot-spot for sire-specific TRD in the cattle genome, and they are related to embryonic and postnatal lethality as well as prostate cancer, among others. This research characterized the distribution of TRD in the bovine genome, highlighting heterogeneous results when comparing across breeds.

Temporal Genetic Variance and Propagule-Driven Genetic Structure Characterize Naturalized Rainbow Trout (Oncorhynchus mykiss) from a Patagonian Lake Impacted by Trout Farming

Knowledge about the genetic underpinnings of invasions—a theme addressed by invasion genetics as a discipline—is still scarce amid well documented ecological impacts of non-native species on ecosystems of Patagonia in South America. One of the most invasive species in Patagonia’s freshwater systems and elsewhere is rainbow trout (Oncorhynchus mykiss). This species was introduced to Chile during the early twentieth century for stocking and promoting recreational fishing; during the late twentieth century was reintroduced for farming purposes and is now naturalized. We used population- and individual-based inference from single nucleotide polymorphisms (SNPs) to illuminate three objectives related to the establishment and naturalization of Rainbow Trout in Lake Llanquihue. This lake has been intensively used for trout farming during the last three decades. Our results emanate from samples collected from five inlet streams over two seasons, winter and spring. First, we found that significant intra- population (temporal) genetic variance was greater than inter-population (spatial) genetic variance, downplaying the importance of spatial divergence during the process of naturalization. Allele frequency differences between cohorts, consistent with variation in fish length between spring and winter collections, might explain temporal genetic differences. Second, individual-based Bayesian clustering suggested that genetic structure within Lake Llanquihue was largely driven by putative farm propagules found at one single stream during spring, but not in winter. This suggests that farm broodstock might migrate upstream to breed during spring at that particular stream. It is unclear whether interbreeding has occurred between “pure” naturalized and farm trout in this and other streams. Third, estimates of the annual number of breeders (N_b) were below 73 in half of the collections, suggestive of genetically small and recently founded populations that might experience substantial genetic drift. Our results reinforce the notion that naturalized trout originated recently from a small yet genetically diverse source and that farm propagules might have played a significant role in the invasion of Rainbow Trout within a single lake with intensive trout farming. Our results also argue for proficient mitigation measures that include management of escapes and strategies to minimize unintentional releases from farm facilities.

High-density interspecific genetic linkage mapping provides insights into genomic incompatibility between channel catfish and blue catfish

Catfish is the leading aquaculture species in the United States. The interspecific hybrid catfish produced by mating female channel catfish with male blue catfish outperform both of their parent species in a number of traits. However, mass production of the hybrids has been difficult because of reproductive isolation. Investigations of genome structure and organization of the hybrids provide insights into the genetic basis for maintenance of species divergence in the face of gene flow, thereby helping develop strategies for introgression and efficient production of the hybrids for aquaculture. In this study, we constructed a high-density genetic linkage map using the hybrid catfish system with the catfish 250K SNP array. A total of 26,238 SNPs were mapped to 29 linkage groups, with 12,776 unique marker positions. The linkage map spans approximately 3240 cM with an average intermarker distance of 0.25 cM. A fraction of markers (986 of 12,776) exhibited significant deviation from the expected Mendelian ratio of segregation, and they were clustered in major genomic blocks across 15 LGs, most notably LG9 and LG15. The distorted markers exhibited significant bias for maternal alleles among the backcross progenies, suggesting strong selection against the blue catfish alleles. The clustering of distorted markers within genomic blocks should lend insights into speciation as marked by incompatibilities between the two species. Such findings should also have profound implications for understanding the genomic evolution of closely related species as well as the introgression of hybrid production programs in aquaculture.

Co-localization of growth QTL with differentially expressed candidate genes in rainbow trout

We tested whether genes differentially expressed between large and small rainbow trout co-localized with familial QTL regions for body size. Eleven chromosomes, known from previous work to house QTL for weight and length in rainbow trout, were examined for QTL in half-sibling families produced in September (1 XY male and 1 XX neomale) and December (1 XY male). In previous studies, we identified 108 candidate genes for growth expressed in the liver and white muscle in a subset of the fish used in this study. These gene sequences were BLASTN aligned against the rainbow trout and stickleback genomes to determine their location (rainbow trout) and inferred location based on synteny with the stickleback genome. Across the progeny of all three males used in the study, 63.9% of the genes with differential expression appear to co-localize with the QTL regions on 6 of the 11 chromosomes tested in these males. Genes that co-localized with QTL in the mixed-sex offspring of the two XY males primarily showed up-regulation in the muscle of large fish and were related to muscle growth, metabolism, and the stress response.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high-density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.

Repetitive Sequence and Sex Chromosome Evolution in Vertebrates

Sex chromosomes are the most dynamic entity in any genome having unique morphology, gene content, and evolution. They have evolved multiple times and independently throughout vertebrate evolution. One of the major genomic changes that pertain to sex chromosomes involves the amplification of common repeats. It is hypothesized that such amplification of repeats facilitates the suppression of recombination, leading to the evolution of heteromorphic sex chromosomes through genetic degradation of Y or W chromosomes. Although contrasting evidence is available, it is clear that amplification of simple repetitive sequences played a major role in the evolution of Y and W chromosomes in vertebrates. In this review, we present a brief overview of the repetitive DNA classes that accumulated during sex chromosome evolution, mainly focusing on vertebrates, and discuss their possible role and potential function in this process.

Reduced recombination patterns in Robertsonian hybrids between chromosomal races of the house mouse: chiasma analyses

The recombination suppression models of chromosomal speciation posit that chromosomal rearrangements act as partial barriers to gene flow allowing these regions to accumulate genetic incompatibilities, thus contributing to the divergence of populations. Empirical and theoretical studies exploring the requirements of these models have mostly focused on the role of inversions. Here, the recombination landscape of heterozygosity for Robertsonian (Rb) fusions is investigated in the house mouse. Laboratory-bred F1 males and females between highly differentiated races from Tunisia (Rb: 2n=22, Standard, St: 2n=40) were produced in which all Rb fusions are present as trivalents in meiosis. Recombination patterns were determined by the analysis of chiasmata and compared with previous data on the Tunisian parental mice. A comparative analysis was performed on wild-caught male mice spanning the hybrid zone between two Italian races (2n=40, 2n=22). The results showed that the chiasma characteristics of both male and female Tunisian F1 and Italian hybrids clearly differed from those of Rb and St mice. Not only was the mean chiasma number (CN) intermediate between those of the parental mice in both geographic samples, but the distribution of chiasmata along the chromosomal arms of the F1 showed a distinct mosaic pattern. In short, the proximal region in the F1 exhibited a reduced CN similar to that observed in homozygous Rb, whereas distal regions more closely matched those in St mice. These results suggest that Rb rearrangements (homozygous or heterozygous) reduce recombination in the proximal regions of the chromosomes supporting their potential role in recombination-mediated speciation models.

Comparative mapping between Coho Salmon (Oncorhynchus kisutch) and three other salmonids suggests a role for chromosomal rearrangements in the retention of duplicated regions following a whole genome duplication event

Whole genome duplication has been implicated in evolutionary innovation and rapid diversification. In salmonid fishes, however, whole genome duplication significantly pre-dates major transitions across the family, and re-diploidization has been a gradual process between genomes that have remained essentially collinear. Nevertheless, pairs of duplicated chromosome arms have diverged at different rates from each other, suggesting that the retention of duplicated regions through occasional pairing between homeologous chromosomes may have played an evolutionary role across species pairs. Extensive chromosomal arm rearrangements have been a key mechanism involved in re-dipliodization of the salmonid genome; therefore, we investigated their influence on degree of differentiation between homeologs across salmon species. We derived a linkage map for coho salmon and performed comparative mapping across syntenic arms within the genus Oncorhynchus, and with the genus Salmo, to determine the phylogenetic relationship between chromosome arrangements and the retention of undifferentiated duplicated regions. A 6596.7 cM female coho salmon map, comprising 30 linkage groups with 7415 and 1266 nonduplicated and duplicated loci, respectively, revealed uneven distribution of duplicated loci along and between chromosome arms. These duplicated regions were conserved across syntenic arms across Oncorhynchus species and were identified in metacentric chromosomes likely formed ancestrally to the divergence of Oncorhynchus from Salmo. These findings support previous studies in which observed pairings involved at least one metacentric chromosome. Re-diploidization in salmon may have been prevented or retarded by the formation of metacentric chromosomes after the whole genome duplication event and may explain lineage-specific innovations in salmon species if functional genes are found in these regions.

Comparative genome mapping between Chinook salmon (Oncorhynchus tshawytscha) and rainbow trout (O. mykiss) based on homologous microsatellite loci

Comparative genome mapping can rapidly facilitate the transfer of DNA sequence information from a well-characterized species to one that is less described. Chromosome arm numbers are conserved between members of the teleost family Salmonidae, order Salmoniformes, permitting rapid alignment of large syntenic blocks of DNA between members of the group. However, extensive Robertsonian rearrangements after an ancestral whole-genome duplication event has resulted in different chromosome numbers across Salmonid taxa. In anticipation of the rapid application of genomic data across members of the Pacific salmon genus Oncorhynchus, we mapped the genome of Chinook salmon (O. tshawytscha) by using 361 microsatellite loci and compared linkage groups to those already derived for a well-characterized species rainbow trout (O. mykiss). The Chinook salmon female map length was 1526 cM, the male map 733 cM, and the consensus map between the two sexes was 2206 cM. The average female to male recombination ratio was 5.43 (range 1-42.8 across all pairwise marker comparisons). We detected 34 linkage groups that corresponded with all chromosome arms mapped with homologous loci in rainbow trout and inferred that 16 represented metacentric chromosomes and 18 represented acrocentric chromosomes. Up to 13 chromosomes were conserved between the two species, suggesting that their structure precedes the divergence between Chinook salmon and rainbow trout. However, marker order differed in one of these linkage groups. The remaining linkage group structures reflected independent Robertsonian chromosomal arrangements, possibly after divergence. The putative linkage group homologies presented here are expected to facilitate future DNA sequencing efforts in Chinook salmon.