Centromere localization for Bighead Carp (Aristichthys nobilis) through half-tetrad analysis in diploid gynogenetic families

Comparative transcriptome analyses and identification of candidate genes involved in vertebral abnormality of bighead carp Hypophthalmichthys nobilis

Body deformity occurs both in wild and farmed fishes, which is one of the most challenging problems for aquaculture industry. In most cases, such body deformities are linked to skeletal deformities. Currently, very limited information is available on skeletal deformities of farmed fish species which may be caused by genetic factor. In this study, we performed muscle and vertebra transcriptome analyses in body deformity and normality of bighead carp Hypophthalmichthys nobilis (from one meiotic gynogenesis family) using RNA-Seq. A total of 43,923 and 44,416 unigenes were predicted in muscles and vertebrae, respectively. Based on these data, we further explored the gene expression profiles in gynogenetic normal and abnormal bighead carp. No differentially expressed gene (DEG) was found in transcriptome data of muscles. Totally, 20 key DEGs were identified in transcriptome data of vertebrae, such as low density lipoprotein-related protein 2 (lrp2), bone morphogenetic protein 2B (bmp2b) and collagen alpha-1(IV) (col4a1). 12 potential pathways were also identified in vertebra transcriptome data, which were mainly involved in development, growth, cytoskeleton and energy metabolism, such as MAPK signaling pathway, regulation of actin cytoskeleton and TGF-beta signaling pathway. Results of this study will be informative for the understanding of genetic mechanisms for body shape formation and also provide potential candidate genes for selection program involved in body shape and skeletal development in H. nobilis.

Construction of High-Resolution RAD-Seq Based Linkage Map, Anchoring Reference Genome, and QTL Mapping of the Sex Chromosome in the Marine Medaka Oryzias melastigma

Medaka (Oryzias sp.) is an important fish species in ecotoxicology and considered as a model species due to its biological features including small body size and short generation time. Since Japanese medaka Oryzias latipes is a freshwater species with access to an excellent genome resource, the marine medaka Oryzias melastigma is also applicable for the marine ecotoxicology. In genome era, a high-density genetic linkage map is a very useful resource in genomic research, providing a means for comparative genomic analysis and verification of de novo genome assembly. In this study, we developed a high-density genetic linkage map for O. melastigma using restriction-site associated DNA sequencing (RAD-seq). The genetic map consisted of 24 linkage groups with 2,481 single nucleotide polymorphism (SNP) markers. The total map length was 1,784 cM with an average marker space of 0.72 cM. The genetic map was integrated with the reference-assisted chromosome assembly (RACA) of O. melastigma, which anchored 90.7% of the assembled sequence onto the linkage map. The values of complete Benchmarking Universal Single-Copy Orthologs were similar to RACA assembly but N50 (23.74 Mb; total genome length 779.4 Mb; gap 5.29%) increased to 29.99 Mb (total genome length 778.7 Mb; gap 5.2%). Using MapQTL analysis with SNP markers, we identified a major quantitative trait locus for sex traits on the Om10. The integration of the genetic map with the reference genome of marine medaka will serve as a good resource for studies in molecular toxicology, genomics, CRISPR/Cas9, and epigenetics.

Centromere localization in medaka fish based on half-tetrad analysis

Gene-centromere (G-C) mapping provides insight into vertebrate genome composition, structure and evolution. Although medaka fish are important experimental animals, no genome-wide G-C map of medaka has been constructed. In this study, we used 112 interspecific triploid hybrids and 152 DNA markers to make G-C maps of all 24 linkage groups (LGs). Under the assumption of 50% interference, 24 centromeres were localized onto all corresponding medaka LGs. Comparison with 21 centromere positions deduced from putative centromeric repeats revealed that 19 were localized inside the centromeric regions of the G-C maps, whereas two were not. Based on the centromere positions indicated in the G-C maps and those of centromeric repeats on each LG, we classified chromosomes as either biarmed or monoarmed; n = 24 = 10 metacentrics/submetacentrics + 14 subtelocentrics/acrocentrics, which is consistent with the results of previous karyological reports. This study helps to elucidate genome evolution mechanisms, and integrates physical and genetic maps with karyological information of medaka.

Sex-specific markers developed by next-generation sequencing confirmed an XX/XY sex determination system in bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix)

Sex-specific markers are powerful tools for identifying sex-determination system in various animals. Bighead carp (Hypophthalmichehys nobilis) and silver carp (Hypophthalmichthys molitrix) are two of the most important edible fish in Asia, which have a long juvenility period that can lasts for 4-5 years. In this study, we found one sex-specific marker by next-generation sequencing together with bioinformatics analysis in bighead carp. The male-specific markers were used to perform molecular sexing in the progenies of artificial gynogenetic diploids and found all progenies (n = 160) were females. Meanwhile, around 1 : 1 sex ratio was observed in a total of 579 juvenile offspring from three other families. To further extend the male-specific region, we performed genome walking and got a male-specific sequence of 8,661 bp. Five pairs of primers were designed and could be used to efficiently distinguish males from females in bighead carp and silver carp. The development of these male-specific markers and results of their molecular sexing in different populations provide strong evidence for a sex determination system of female homogametry or male heterogametry (XX/XY) in bighead carp and silver carp. To the best of our knowledge, this is the first report of effective sex-specific markers in these two large carp species.

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.

Genetic inviability is a major driver of type III survivorship in experimental families of a highly fecund marine bivalve

The offspring of most highly fecund marine fish and shellfish suffer substantial mortality early in the life cycle, complicating prediction of recruitment and fisheries management. Early mortality has long been attributed to environmental factors and almost never to genetic sources. Previous work on a variety of marine bivalve species uncovered substantial genetic inviability among the offspring of inbred crosses, suggesting a large load of early-acting deleterious recessive mutations. However, genetic inviability of randomly bred offspring has not been addressed. Here, genome-wide surveys reveal widespread, genotype-dependent mortality in randomly bred, full-sib progenies of wild-caught Pacific oysters (Crassostrea gigas). Using gene-mapping methods, we infer that 11-19 detrimental alleles per family render 97.9-99.8% of progeny inviable. The variable genomic positions of viability loci among families imply a surprisingly large load of partially dominant or additive detrimental mutations in wild adult oysters. Although caution is required in interpreting the relevance of experimental results for natural field environments, we argue that the observed genetic inviability corresponds with type III survivorship, which is characteristic of both hatchery and field environments and that our results, therefore, suggest the need for additional experiments under the near-natural conditions of mesocosms. We explore the population genetic implications of our results, calculating a detrimental mutation rate that is comparable to that estimated for conifers and other highly fecund perennial plants. Genetic inviability ought to be considered as a potential major source of low and variable recruitment in highly fecund marine animals.

Uncovering Cryptic Asexuality in Daphnia magna by RAD Sequencing

The breeding systems of many organisms are cryptic and difficult to investigate with observational data, yet they have profound effects on a species' ecology, evolution, and genome organization. Genomic approaches offer a novel, indirect way to investigate breeding systems, specifically by studying the transmission of genetic information from parents to offspring. Here we exemplify this method through an assessment of self-fertilization vs. automictic parthenogenesis in Daphnia magna. Self-fertilization reduces heterozygosity by 50% compared to the parents, but under automixis, whereby two haploid products from a single meiosis fuse, the expected heterozygosity reduction depends on whether the two meiotic products are separated during meiosis I or II (i.e., central vs. terminal fusion). Reviewing the existing literature and incorporating recombination interference, we derive an interchromosomal and an intrachromosomal prediction of how to distinguish various forms of automixis from self-fertilization using offspring heterozygosity data. We then test these predictions using RAD-sequencing data on presumed automictic diapause offspring of so-called nonmale producing strains and compare them with "self-fertilized" offspring produced by within-clone mating. The results unequivocally show that these offspring were produced by automixis, mostly, but not exclusively, through terminal fusion. However, the results also show that this conclusion was only possible owing to genome-wide heterozygosity data, with phenotypic data as well as data from microsatellite markers yielding inconclusive or even misleading results. Our study thus demonstrates how to use the power of genomic approaches for elucidating breeding systems, and it provides the first demonstration of automictic parthenogenesis in Daphnia.

Genetic mapping of centromeres in the nine Citrus clementina chromosomes using half-tetrad analysis and recombination patterns in unreduced and haploid gametes

BACKGROUND

Mapping centromere locations in plant species provides essential information for the analysis of genetic structures and population dynamics. The centromere's position affects the distribution of crossovers along a chromosome and the parental heterozygosity restitution by 2n gametes is a direct function of the genetic distance to the centromere. Sexual polyploidisation is relatively frequent in Citrus species and is widely used to develop new seedless triploid cultivars. The study's objectives were to (i) map the positions of the centromeres of the nine Citrus clementina chromosomes; (ii) analyse the crossover interference in unreduced gametes; and (iii) establish the pattern of genetic recombination in haploid clementine gametes along each chromosome and its relationship with the centromere location and distribution of genic sequences.

RESULTS

Triploid progenies were derived from unreduced megagametophytes produced by second-division restitution. Centromere positions were mapped genetically for all linkage groups using half-tetrad analysis. Inference of the physical locations of centromeres revealed one acrocentric, four metacentric and four submetacentric chromosomes. Crossover interference was observed in unreduced gametes, with variation seen between chromosome arms. For haploid gametes, a strong decrease in the recombination rate occurred in centromeric and pericentromeric regions, which contained a low density of genic sequences. In chromosomes VIII and IX, these low recombination rates extended beyond the pericentromeric regions. The genomic region corresponding to a genetic distance < 5cM from a centromere represented 47% of the genome and 23% of the genic sequences.

CONCLUSIONS

The centromere positions of the nine citrus chromosomes were genetically mapped. Their physical locations, inferred from the genetic ones, were consistent with the sequence constitution and recombination pattern along each chromosome. However, regions with low recombination rates extended beyond the pericentromeric regions of some chromosomes into areas richer in genic sequences. The persistence of strong linkage disequilibrium between large numbers of genes promotes the stability of epistatic interactions and multilocus-controlled traits over successive generations but also maintains multi-trait associations. Identification of the centromere positions will allow the development of simple methods to analyse unreduced gamete formation mechanisms in a large range of genotypes and further modelling of genetic inheritance in sexual polyploidisation breeding schemes.

Microsatellite-centromere mapping in common carp through half-tetrad analysis in diploid meiogynogenetic families

Gene-centromere (G-C) mapping provides insights into the understanding of the composition, structure, and evolution of vertebrate genomes. Common carp (Cyprinus carpio) is an important aquaculture fish and has been proposed to undertake tetraploidization. In this study, we selected 214 informative microsatellite markers across 50 linkage groups of a common carp genetic map to perform gene-centromere mapping using half-tetrad analysis. A total of 199 microsatellites were segregated under the Mendelian expectations in at least one of the three gynogenetic families and were used for G-C distance estimation. The G-C recombination frequency (y) ranged from 0 to 0.99 (0.43 on average), corresponding to a fixation index (F) of 0.57 after one generation of gynogenesis. Large y values for some loci together with significant correlation between G-C distances and genetic linkage map distances suggested the presence of high interference in common carp. Under the assumption of complete interference, 50 centromeres were localized onto corresponding linkage groups (LGs) of common carp, with G-C distances of centromere-linked markers per LG ranging from 0 to 10.3 cM (2.9 cM on average). Based on the information for centromere positions, we proposed a chromosome formula of 2n = 100 = 58 m/sm + 42 t/st with 158 chromosome arms for common carp, which was similar to a study observed by cytogenetic method. The examination of crossover distributions along 10 LGs revealed that the proportion of crossover chromatids was overall higher than that of non-crossover chromatids in gynogenetic progenies, indicating high recombination levels across most LGs. Comparative genomics analyses suggested that the chromosomes of common carp have undergone extensive rearrangement after genome duplication. This study would be valuable to elucidate the mechanism of genome evolution and integrate physical and genetic maps in common carp.