A high-density linkage map for Astyanax mexicanus using genotyping-by-sequencing technology

Quantitative trait loci concentrate in specific regions of the Mexican cavefish genome and reveal key candidate genes for cave-associated evolution

A major goal of modern biology is connecting phenotype with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here, we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. These analyses revealed that QTL clusters in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics associated with increased opportunities for mutation, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb, which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific regions of the genome may play significant roles in driving adaptation to the cave environment in A. mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Extraordinary model systems for regeneration

Regeneration is the remarkable phenomenon through which an organism can regrow lost or damaged parts with fully functional replacements, including complex anatomical structures, such as limbs. In 2019, Development launched its 'Model systems for regeneration' collection, a series of articles introducing some of the most popular model organisms for studying regeneration in vivo. To expand this topic further, this Perspective conveys the voices of five expert biologists from the field of regenerative biology, each of whom showcases some less well-known, but equally extraordinary, species for studying regeneration.

The spatiotemporal and genetic architecture of extraoral taste buds in Astyanax cavefish

Intense environmental pressures can yield both regressive and constructive traits through complex evolutionary mechanisms. Although regression is well-studied, the biological bases of constructive features are less well understood. Cave-dwelling Astyanax fish harbor prolific extraoral taste buds on their heads, which are absent in conspecific surface-dwellers. Here, we present novel ontogenetic data demonstrating extraoral taste buds appear gradually and late in life history. This appearance is similar but non-identical in different cavefish populations, where patterning has evolved to permit taste bud re-specification across the endoderm-ectoderm germ layer boundary. Quantitative genetic analyses revealed that spatially distinct taste buds on the head are primarily mediated by two different cave-dominant loci. While the precise function of this late expansion on to the head is unknown, the appearance of extraoral taste buds coincides with a dietary shift from live-foods to bat guano, suggesting an adaptive mechanism to detect nutrition in food-starved caves. This work provides fundamental insight to a constructive evolutionary feature, arising late in life history, promising a new window into unresolved features of vertebrate sensory organ development.

Evolution in Sinocyclocheilus cavefish is marked by rate shifts, reversals, and origin of novel traits

Background

Natural model systems are indispensable for exploring adaptations in response to environmental pressures. Sinocyclocheilus of China, the most diverse cavefish clade in the world (75 species), provide unique opportunities to understand recurrent evolution of stereotypic traits (such as eye loss and sensory expansion) in the context of a deep and diverse phylogenetic group. However, they remain poorly understood in terms of their morphological evolution. Therefore, we explore key patterns of morphological evolution, habitat utilization and geographic distribution in these fishes.

Results

We constructed phylogenies and categorized 49 species based on eye-related condition (Blind, Micro-eyed, and Normal-eyed), habitat types (Troglobitic—cave-restricted; Troglophilic—cave-associated; Surface—outside caves) and existence of horns. Geometric-morphometric analyses show Normal-eyed morphs with fusiform shapes segregating from Blind/Micro-eyed deeper bodied morphs along the first principal-component axis; second axis accounts for shape complexity related to horns. The body shapes showed a significant association with eye-related condition and horn, but not habitat types. Ancestral reconstructions suggest at least three independent origins of Blind morphs, each with different levels of modification in relation to their ancestral Normal-eyed morphs; Sinocyclocheilus are also pre-adapted for cave dwelling. Our geophylogeny shows an east-to-west diversification spanning Pliocene and Pleistocene, with early-diversifying Troglobitic species dominating subterranean habitats of karstic plains whereas predominantly Surface forms inhabit hills to the west. Evolutionary rates analyses suggest that lineages leading to Blind morphs were characterized by significant rate shifts, such as a slowdown in body size evolution and a 5–20 fold increase in rate of eye regression, possibly explained by limited resource availability. Body size and eye size have undergone reversals, but not horns, a trait entailing considerable time to form.

Conclusions

Sinocyclocheilus occupied cave habitats in response to drying associated with aridification of China during late Miocene and the Pliocene. The prominent cave-adaptations (eye-regression, horn-evolution) occur in clades associated with the extensive subterranean cave system in Guangxi and Guizhou provinces. Integration of morphology, phylogeny, rate analyses, molecular-dating and distribution show not only several remarkable patterns of evolution, but also interesting exceptions to these patterns signifying the diversification of Sinocyclocheilus as an invaluable model system to explore evolutionary novelty.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01776-y.

A chromosome-level genome of Astyanax mexicanus surface fish for comparing population-specific genetic differences contributing to trait evolution

Identifying the genetic factors that underlie complex traits is central to understanding the mechanistic underpinnings of evolution. Cave-dwelling Astyanax mexicanus populations are well adapted to subterranean life and many populations appear to have evolved troglomorphic traits independently, while the surface-dwelling populations can be used as a proxy for the ancestral form. Here we present a high-resolution, chromosome-level surface fish genome, enabling the first genome-wide comparison between surface fish and cavefish populations. Using this resource, we performed quantitative trait locus (QTL) mapping analyses and found new candidate genes for eye loss such as dusp26. We used CRISPR gene editing in A. mexicanus to confirm the essential role of a gene within an eye size QTL, rx3, in eye formation. We also generated the first genome-wide evaluation of deletion variability across cavefish populations to gain insight into this potential source of cave adaptation. The surface fish genome reference now provides a more complete resource for comparative, functional and genetic studies of drastic trait differences within a species.

An Asymmetric Genetic Signal Associated with Mechanosensory Expansion in Cave-Adapted Fish

A key challenge in contemporary biology is connecting genotypic variation to phenotypic diversity. Quantitative genetics provides a powerful technique for identifying regions of the genome that covary with phenotypic variation. Here, we present a quantitative trait loci (QTL) analysis of a natural freshwater fish system, Astyanax mexicanus, that harbors two morphs corresponding to a cave and surface fish. Following their divergence ~500 Kya, cavefish have adapted to the extreme pressures of the subterranean biome. As a consequence, cavefish have lost numerous features, but evolved gains for a variety of constructive features including behavior. Prior work found that sensory tissues (neuromasts) present in the "eye orbit" region of the skull associate with sensitivity to vibrations in water. This augmented sensation is believed to facilitate foraging behavior in the complete darkness of a cave, and may impact on evolved lateral swimming preference. To this point, however, it has remained unclear how morphological variation integrates with behavioral variation through heritable factors. Using a QTL approach, we discovered the genetic architecture of neuromasts present in the eye orbit region, demonstrating that this feature is under genetic control. Interestingly, linked loci were asymmetric-signals were detected using only data collected from the right, but not left, side of the face. This finding may explain enhanced sensitivity and/or feedback of water movements mediating a lateral swimming preference. The locus we discovered based on neuromast position maps near established QTL for eye size and a facial bone morphology, raising the intriguing possibility that eye loss, sensory expansion, and the cranial skeleton may be integrated for evolving adaptive behaviors. Thus, this work will further our understanding of the functional consequence of key loci that influence the evolutionary origin of changes impacting morphology, behavior, and adaptation.

Heart Regeneration in the Mexican Cavefish

Although Astyanax mexicanus surface fish regenerate their hearts after injury, their Pachón cave-dwelling counterparts cannot and, instead, form a permanent fibrotic scar, similar to the human heart. Myocardial proliferation peaks at similar levels in both surface fish and Pachón 1 week after injury. However, in Pachón, this peak coincides with a strong scarring and immune response, and ultimately, cavefish cardiomyocytes fail to replace the scar. We identified lrrc10 to be upregulated in surface fish compared with Pachón after injury. Similar to cavefish, knockout of lrrc10 in zebrafish impairs heart regeneration without affecting wound cardiomyocyte proliferation. Furthermore, using quantitative trait locus (QTL) analysis, we have linked the degree of heart regeneration to three loci in the genome, identifying candidate genes fundamental to the difference between scarring and regeneration. Our study provides evidence that successful heart regeneration entails a delicate interplay between cardiomyocyte proliferation and scarring.

•

Astyanax mexicanus surface fish regenerate their hearts, but Pachón cavefish cannot

•

Successful regeneration is a delicate interplay between proliferation and scarring

•

lrrc10 is required for heart regeneration

•

Using QTL analysis, we have identified three genomic loci linked to heart regeneration

High-Density Genetic Map Construction and QTL Mapping of Leaf and Needling Traits in Ziziphus jujuba Mill

The Chinese jujube (Ziziphus jujuba Mill., 2n = 2x = 24), one of the most popular fruit trees in Asia, is widely cultivated and utilized in China, where it is traditionally consumed as both a fresh and dried food resource. A high-density genetic map can provide the necessary framework for quantitative trait loci (QTL) analyses and map-based gene cloning and molecular breeding. In this study, we constructed a new high-density genetic linkage map via a genotyping-by-sequencing approach. For the consensus linkage map, a total of 3,792 markers spanning 2,167.5 cM were mapped onto 12 linkage groups, with an average marker interval distance of 0.358 cM. The genetic map anchored 301 Mb (85.7%) of scaffolds from the sequenced Z. jujuba "Junzao" genome. Based on this genetic map, 30 potential QTLs were detected, including 27 QTLs for leaf traits and 3 QTLs for needling length. This high-density genetic map and the identified QTLs for relevant agronomic traits lay the groundwork for functional genetic mapping, map-based cloning, and marker-assisted selection in jujube.

An Adult Brain Atlas Reveals Broad Neuroanatomical Changes in Independently Evolved Populations of Mexican Cavefish

A shift in environmental conditions impacts the evolution of complex developmental and behavioral traits. The Mexican cavefish, Astyanax mexicanus, is a powerful model for examining the evolution of development, physiology, and behavior because multiple cavefish populations can be compared to an extant, ancestral-like surface population of the same species. Many behaviors have diverged in cave populations of A. mexicanus, and previous studies have shown that cavefish have a loss of sleep, reduced stress, an absence of social behaviors, and hyperphagia. Despite these findings, surprisingly little is known about the changes in neuroanatomy that underlie these behavioral phenotypes. Here, we use serial sectioning to generate brain atlases of surface fish and three independent cavefish populations. Volumetric reconstruction of serial-sectioned brains confirms convergent evolution on reduced optic tectum volume in all cavefish populations tested. In addition, we quantified volumes of specific neuroanatomical loci within several brain regions that have previously been implicated in behavioral regulation, including the hypothalamus, thalamus, and habenula. These analyses reveal an enlargement of the hypothalamus in all cavefish populations relative to surface fish, as well as subnuclei-specific differences within the thalamus and prethalamus. Taken together, these analyses support the notion that changes in environmental conditions are accompanied by neuroanatomical changes in brain structures associated with behavior. This atlas provides a resource for comparative neuroanatomy of additional brain regions and the opportunity to associate brain anatomy with evolved changes in behavior.

Stable transgenesis in Astyanax mexicanus using the Tol2 transposase system

BACKGROUND

Astyanax mexicanus is a well-established fish model system for evolutionary and developmental biology research. These fish exist as surface forms that inhabit rivers and 30 different populations of cavefish. Despite important progress in the deployment of new technologies, deep mechanistic insights into the genetic basis of evolution, development, and behavior have been limited by a lack of transgenic lines commonly used in genetic model systems.

RESULTS

Here, we expand the toolkit of transgenesis by characterizing two novel stable transgenic lines that were generated using the highly efficient Tol2 system, commonly used to generate transgenic zebrafish. A stable transgenic line consisting of the zebrafish ubiquitin promoter expresses enhanced green fluorescent protein ubiquitously throughout development in a surface population of Astyanax. To define specific cell-types, a Cntnap2-mCherry construct labels lateral line mechanosensory neurons in zebrafish. Strikingly, both constructs appear to label the predicted cell types, suggesting many genetic tools and defined promoter regions in zebrafish are directly transferrable to cavefish.

CONCLUSION

The lines provide proof-of-principle for the application of Tol2 transgenic technology in A. mexicanus. Expansion on these initial transgenic lines will provide a platform to address broadly important problems in the quest to bridge the genotype-phenotype gap.

Construction of a high-density genetic map: genotyping by sequencing (GBS) to map purple seed coat color (Psc) in hulless barley

Background

Colored hulless barley are more suitable in food processing compared to normal (yellow) varieties because it is rich in bioactive compounds and produces higher extraction pearling fractions. Therefore, seed coat color is an important agronomic trait for the breeding and study of hulless barley.

Results

Genotyping-by-sequencing single-nucleotide polymorphism (GBS-SNP) analysis of a doubled haploid (DH) mapping population (Nierumuzha × Kunlun10) was conducted to map the purple seed coat color genes (Psc). A high-density genetic map of hulless barley was constructed, which contains 3662 efficient SNP markers with 1129 bin markers. Seven linkage groups were resolved, which had a total length of 645.56 cM. Chromosome length ranged from 60.21 cM to 127.21 cM, with average marker density of 0.57 cM. A total of five loci accounting for 3.79% to 23.86% of the observed phenotypic variation for Psc were detected using this high-density map. Five structural candidate genes (F3’M, HID, UF3GT, UFGT and 5MAT) and one regulatory factor (Ant1) related to flavonoid or anthocyanin biosynthesis were identified..

Conclusions

Five structural candidate genes and one regulatory factor related to flavonoid or anthocyanin biosynthesis have been identified using a high-density genetic map of hulless barley. This study lays the foundation for map-based cloning of Psc but provides a valuable tool for studying marker-trait associations and its application to marker-assisted breeding of hulless barley.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0072-6) contains supplementary material, which is available to authorized users.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Genetic analysis reveals candidate genes for activity QTL in the blind Mexican tetra, Astyanax mexicanus

Animal models provide useful tools for exploring the genetic basis of morphological, physiological and behavioral phenotypes. Cave-adapted species are particularly powerful models for a broad array of phenotypic changes with evolutionary, developmental and clinical relevance. Here, we explored the genetic underpinnings of previously characterized differences in locomotor activity patterns between the surface-dwelling and Pachón cave-dwelling populations of Astyanax mexicanus. We identified multiple novel QTL underlying patterns in overall levels of activity (velocity), as well as spatial tank use (time spent near the top or bottom of the tank). Further, we demonstrated that different regions of the genome mediate distinct patterns in velocity and tank usage. We interrogated eight genomic intervals underlying these activity QTL distributed across six linkage groups. In addition, we employed transcriptomic data and draft genomic resources to generate and evaluate a list of 36 potential candidate genes. Interestingly, our data support the candidacy of a number of genes, but do not suggest that differences in the patterns of behavior observed here are the result of alterations to certain candidate genes described in other species (e.g., teleost multiple tissue opsins, melanopsins or members of the core circadian clockwork). This study expands our knowledge of the genetic architecture underlying activity differences in surface and cavefish. Future studies will help define the role of specific genes in shaping complex behavioral phenotypes in Astyanax and other vertebrate taxa.

Uncovering the Genetic Architecture of Seed Weight and Size in Intermediate Wheatgrass through Linkage and Association Mapping

Intermediate wheatgrass [IWG; (Host) Barkworth & D.R. Dewey subsp. ] is being developed as a new perennial grain crop that has a large allohexaploid genome similar to that of wheat ( L.). Breeding for increased seed weight is one of the primary goals for improving grain yield of IWG. As a new crop, however, the genetic architecture of seed weight and size has not been characterized, and selective breeding of IWG may be more intricate than wheat because of its self-incompatible mating system and perennial growth habit. Here, seed weight, seed area size, seed width, and seed length were evaluated across multiple years, in a heterogeneous breeding population comprised of 1126 genets and two clonally replicated biparental populations comprised of 172 and 265 genets. Among 10,171 DNA markers discovered using genotyping-by-sequencing (GBS) in the breeding population, 4731 markers were present in a consensus genetic map previously constructed using seven full-sib populations. Thirty-three quantitative trait loci (QTL) associated with seed weight and size were identified using association mapping (AM), of which 23 were verified using linkage mapping in the biparental populations. About 37.6% of seed weight variation in the breeding population was explained by 15 QTL, 12 of which also contributed to either seed length or seed width. When performing either phenotypic selection or genomic selection for seed weight, we observed the frequency of favorable QTL alleles were increased to >46%. Thus, by combining AM and genomic selection, we can effectively select the favorable QTL alleles for seed weight and size in IWG breeding populations.

A High-Density Genetic Linkage Map and QTL Fine Mapping for Body Weight in Crucian Carp (Carassius auratus) Using 2b-RAD Sequencing

A high-resolution genetic linkage map is essential for a wide range of genetics and genomics studies such as comparative genomics analysis and QTL fine mapping. Crucian carp (Carassius auratus) is widely distributed in Eurasia, and is an important aquaculture fish worldwide. In this study, a high-density genetic linkage map was constructed for crucian carp using 2b-RAD technology. The consensus map contains 8487 SNP markers, assigning to 50 linkage groups (LGs) and spanning 3762.88 cM, with an average marker interval of 0.44 cM and genome coverage of 98.8%. The female map had 4410 SNPs, and spanned 3500.42 cM (0.79 cM/marker), while the male map had 4625 SNPs and spanned 3346.33 cM (0.72 cM/marker). The average recombination ratio of female to male was 2.13:1, and significant male-biased recombination suppressions were observed in LG47 and LG49. Comparative genomics analysis revealed a clear 2:1 syntenic relationship between crucian carp LGs and chromosomes of zebrafish and grass carp, and a 1:1 correspondence, but extensive chromosomal rearrangement, between crucian carp and common carp, providing evidence that crucian carp has experienced a fourth round of whole genome duplication (4R-WGD). Eight chromosome-wide QTL for body weight at 2 months after hatch were detected on five LGs, explaining 10.1-13.2% of the phenotypic variations. Potential candidate growth-related genes, such as an EGF-like domain and TGF-β, were identified within the QTL intervals. This high-density genetic map and QTL analysis supplies a basis for genome evolutionary studies in cyprinid fishes, genome assembly, and QTL fine mapping for complex traits in crucian carp.

Evolution and conservation of Characidium sex chromosomes

Fish species exhibit substantial variation in the degree of genetic differentiation between sex chromosome pairs, and therefore offer the opportunity to study the full range of sex chromosome evolution. We used restriction-site associated DNA sequencing (RAD-seq) to study the sex chromosomes of Characidium gomesi, a species with conspicuous heteromorphic ZW/ZZ sex chromosomes. We screened 9863 single-nucleotide polymorphisms (SNPs), corresponding to ~1 marker/100 kb distributed across the genome for sex-linked variation. With this data set, we identified 26 female-specific RAD loci, putatively located on the W chromosome, as well as 148 sex-associated SNPs showing significant differentiation (average F_ST=0.144) between males and females, and therefore in regions of more recent divergence between the Z and W chromosomes. In addition, we detected 25 RAD loci showing extreme heterozygote deficiency in females but which were in Hardy-Weinberg equilibrium in males, consistent with degeneration of the W chromosome and therefore female hemizygosity. We validated seven female-specific and two sex-associated markers in a larger sample of C. gomesi, of which three localised to the W chromosome, thereby providing useful markers for sexing wild samples. Validated markers were evaluated in other populations and species of the genus Characidium, this exploration suggesting a rapid turnover of W-specific repetitive elements. Together, our analyses point to a complex origin for the sex chromosome of C. gomesi and highlight the utility of RAD-seq for studying the composition and evolution of sex chromosomes systems in wild populations.

Large-scale SNP discovery and construction of a high-density genetic map of Colossoma macropomum through genotyping-by-sequencing

Colossoma macropomum, or tambaqui, is the largest native Characiform species found in the Amazon and Orinoco river basins, yet few resources for genetic studies and the genetic improvement of tambaqui exist. In this study, we identified a large number of single-nucleotide polymorphisms (SNPs) for tambaqui and constructed a high-resolution genetic linkage map from a full-sib family of 124 individuals and their parents using the genotyping by sequencing method. In all, 68,584 SNPs were initially identified using minimum minor allele frequency (MAF) of 5%. Filtering parameters were used to select high-quality markers for linkage analysis. We selected 7,734 SNPs for linkage mapping, resulting in 27 linkage groups with a minimum logarithm of odds (LOD) of 8 and maximum recombination fraction of 0.35. The final genetic map contains 7,192 successfully mapped markers that span a total of 2,811 cM, with an average marker interval of 0.39 cM. Comparative genomic analysis between tambaqui and zebrafish revealed variable levels of genomic conservation across the 27 linkage groups which allowed for functional SNP annotations. The large-scale SNP discovery obtained here, allowed us to build a high-density linkage map in tambaqui, which will be useful to enhance genetic studies that can be applied in breeding programs.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Manila clam Ruditapes philippinarum

Genetic linkage maps are indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as mollusks. In this study, high-density linkage map was constructed for an ecologically and commercially important clam species, Ruditapes philippinarum. For the consensus linkage map, a total of 9658 markers spanning 1926.98 cM were mapped to 18 sex-averaged linkage groups, with an average marker distance of 0.42 cM. Based on the high-density linkage map, ten QTLs for growth-related traits and shell color were detected. The coverage and density of the current map are sufficient for us to effectively detect QTL for segregating traits, and two QTL positions were all coincident with the closest markers. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of R. philippinarum and other bivalve molluscan species.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Evolutionary Genetics of the Cavefish Astyanax mexicanus

Blind and depigmented fish belonging to the species Astyanax mexicanus are outstanding models for evolutionary genetics. During their evolution in the darkness of caves, they have undergone a number of changes at the morphological, physiological, and behavioral levels, but they can still breed with their river-dwelling conspecifics. The fertile hybrids between these two morphotypes allow forward genetic approaches, from the search of quantitative trait loci to the identification of the mutations underlying the evolution of troglomorphism. We review here the past 30years of evolutionary genetics on Astyanax: from the first crosses and the discovery of convergent evolution of different Astyanax cavefish populations to the most recent evolutionary transcriptomics and genomics studies that have provided researchers with potential candidate genes to be tested using functional genetic approaches. Although significant progress has been made and some genes have been identified, cavefish have not yet fully revealed the secret of their adaptation to the absence of light. In particular, the genetic determinism of their loss of eyes seems complex and still puzzles researchers. We also discuss future research directions, including searches for the origin of cave alleles and searches for selection genome-wide, as well as the necessary but missing information on the timing of cave colonization by surface fish.

Establishment and Optimization of Genomic Selection to Accelerate the Domestication and Improvement of Intermediate Wheatgrass

Intermediate wheatgrass (IWG) is a perennial species and has edible and nutritious grain and desirable agronomic traits, including large seed size, high grain yield, and biomass. It also has the potential to provide ecosystem services and an economic return to farmers. However, because of its allohexaploidy and self-incompatibility, developing molecular markers for genetic analysis and molecular breeding has been challenging. In the present study, using genotyping-by-sequencing (GBS) technology, 3436 genome-wide markers discovered in a biparental population with 178 genets, were mapped to 21 linkage groups (LG) corresponding to 21 chromosomes of IWG. Genomic prediction models were developed using 3883 markers discovered in a breeding population containing 1126 representative genets from 58 half-sib families. High predictive ability was observed for seven agronomic traits using cross-validation, ranging from 0.46 for biomass to 0.67 for seed weight. Optimization results indicated that 8 to 10 genets from each half-sib family can form a good training population to predict the breeding value of their siblings, and 1600 genome-wide markers are adequate to capture the genetic variation in the current breeding population for genomic selection. Thus, with the advances in sequencing-based marker technologies, it was practical to perform molecular genetic analysis and molecular breeding on a new and challenging species like IWG, and genomic selection could increase the efficiency of recurrent selection and accelerate the domestication and improvement of IWG.

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres - an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half-tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination - a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high-density maps in species with strong recombination interference and will enrich many existing and future mapping resources.

A Transcriptomic Analysis of Cave, Surface, and Hybrid Isopod Crustaceans of the Species Asellus aquaticus

Cave animals, compared to surface-dwelling relatives, tend to have reduced eyes and pigment, longer appendages, and enhanced mechanosensory structures. Pressing questions include how certain cave-related traits are gained and lost, and if they originate through the same or different genetic programs in independent lineages. An excellent system for exploring these questions is the isopod, Asellus aquaticus. This species includes multiple cave and surface populations that have numerous morphological differences between them. A key feature is that hybrids between cave and surface individuals are viable, which enables genetic crosses and linkage analyses. Here, we advance this system by analyzing single animal transcriptomes of Asellus aquaticus. We use high throughput sequencing of non-normalized cDNA derived from the head of a surface-dwelling male, the head of a cave-dwelling male, the head of a hybrid male (produced by crossing a surface individual with a cave individual), and a pooled sample of surface embryos and hatchlings. Assembling reads from surface and cave head RNA pools yielded an integrated transcriptome comprised of 23,984 contigs. Using this integrated assembly as a reference transcriptome, we aligned reads from surface-, cave- and hybrid- head tissue and pooled surface embryos and hatchlings. Our approach identified 742 SNPs and placed four new candidate genes to an existing linkage map for A. aquaticus. In addition, we examined SNPs for allele-specific expression differences in the hybrid individual. All of these resources will facilitate identification of genes and associated changes responsible for cave adaptation in A. aquaticus and, in concert with analyses of other species, will inform our understanding of the evolutionary processes accompanying adaptation to the subterranean environment.

The rise of Astyanax cavefish

Numerous animals have invaded subterranean caverns and evolved remarkably similar features. These features include loss of vision and pigmentation, and gains in nonvisual sensation. This broad convergence echoes smaller-scale convergence, in which members of the same species repeatedly evolve the same cave-associated phenotypes. The blind Mexican tetra of the Sierra de El Abra region of northeastern Mexico has a complex origin, having recurrently colonized subterranean environments through numerous invasions of surface-dwelling fish. These colonizations likely occurred ∼1-5 MYa. Despite evidence of historical and contemporary gene flow between cave and surface forms, the cave-associated phenotype appears to remain quite stable in nature. This model system has provided insight to the mechanisms of phenotypic regression, the genetic basis for constructive trait evolution, and the origin of behavioral novelties. Here, we document the rise of this model system from its discovery by a Mexican surveyor in 1936, to a powerful system for cave biology and contemporary genetic research. The recently sequenced genome provides exciting opportunities for future research, and will help resolve several long-standing biological problems. Developmental Dynamics 244:1031-1038, 2015. © 2015 Wiley Periodicals, Inc.

A nonclassical MHC class I U lineage locus in zebrafish with a null haplotypic variant

Three sequence lineages of MHC class I genes have been described in zebrafish (Danio rerio): U, Z, and L. The U lineage genes encoded on zebrafish chromosome 19 are predicted to provide the classical function of antigen presentation. This MHC class I locus displays significant haplotypic variation and is the only MHC class I locus in zebrafish that shares conserved synteny with the core mammalian MHC. Here, we describe two MHC class I U lineage genes, mhc1ula and mhc1uma, that map to chromosome 22. Unlike the U lineage proteins encoded on chromosome 19, Ula and Uma likely play a nonclassical role as they lack conservation of key peptide binding residues, display limited polymorphic variation, and exhibit tissue-specific expression. We also describe a null haplotype at this chromosome 22 locus in which the mhc1ula and mhc1uma genes are absent due to a ~30 kb deletion with no other MHC class I sequences present. Functional and non-functional transcripts of mhc1ula and mhc1uma were identified; however, mhc1uma transcripts were often not amplified or amplified at low levels from individuals possessing an apparently bona fide gene. These distinct U lineage genes may be restricted to the superorder Ostariophysi as similar sequences only could be identified from the blind cavefish (Astyanax mexicanus), fathead minnow (Pimephales promelas), goldfish (Carassius auratus), and grass carp (Ctenopharyngodon idella).

Genome Assembly Improvement and Mapping Convergently Evolved Skeletal Traits in Sticklebacks with Genotyping-by-Sequencing

Marine populations of the threespine stickleback (Gasterosteus aculeatus) have repeatedly colonized and rapidly adapted to freshwater habitats, providing a powerful system to map the genetic architecture of evolved traits. Here, we developed and applied a binned genotyping-by-sequencing (GBS) method to build dense genome-wide linkage maps of sticklebacks using two large marine by freshwater F2 crosses of more than 350 fish each. The resulting linkage maps significantly improve the genome assembly by anchoring 78 new scaffolds to chromosomes, reorienting 40 scaffolds, and rearranging scaffolds in 4 locations. In the revised genome assembly, 94.6% of the assembly was anchored to a chromosome. To assess linkage map quality, we mapped quantitative trait loci (QTL) controlling lateral plate number, which mapped as expected to a 200-kb genomic region containing Ectodysplasin, as well as a chromosome 7 QTL overlapping a previously identified modifier QTL. Finally, we mapped eight QTL controlling convergently evolved reductions in gill raker length in the two crosses, which revealed that this classic adaptive trait has a surprisingly modular and nonparallel genetic basis.