Chromosome-level genome sequence of the Genetically Improved Farmed Tilapia (GIFT, Oreochromis niloticus) highlights regions of introgression with O. mossambicus

BACKGROUND

The Nile tilapia (Oreochromis niloticus) is the third most important freshwater fish for aquaculture. Its success is directly linked to continuous breeding efforts focusing on production traits such as growth rate and weight. Among those elite strains, the Genetically Improved Farmed Tilapia (GIFT) programme initiated by WorldFish is now distributed worldwide. To accelerate the development of the GIFT strain through genomic selection, a high-quality reference genome is necessary.

RESULTS

Using a combination of short (10X Genomics) and long read (PacBio HiFi, PacBio CLR) sequencing and a genetic map for the GIFT strain, we generated a chromosome level genome assembly for the GIFT. Using genomes of two closely related species (O. mossambicus, O. aureus), we characterised the extent of introgression between these species and O. niloticus that has occurred during the breeding process. Over 11 Mb of O. mossambicus genomic material could be identified within the GIFT genome, including genes associated with immunity but also with traits of interest such as growth rate.

CONCLUSION

Because of the breeding history of elite strains, current reference genomes might not be the most suitable to support further studies into the GIFT strain. We generated a chromosome level assembly of the GIFT strain, characterising its mixed origins, and the potential contributions of introgressed regions to selected traits.

Breeding for robustness: investigating the genotype-by-environment interaction and micro-environmental sensitivity of Genetically Improved Farmed Tilapia (Oreochromis niloticus)

Robustness has become a highly desirable breeding goal in the globalized agricultural market. Both genotype‐by‐environment interaction (G × E) and micro‐environmental sensitivity are important robustness components of aquaculture production, in which breeding stock is often disseminated to different environments. The objectives of this study were (i) to quantify the degree of G × E by assessing the growth performance of Genetically Improved Farmed Tilapia (GIFT) across three countries (Malaysia, India and China) and (ii) to quantify the genetic heterogeneity of environmental variance for body weight at harvest (BW) in GIFT as a measure of micro‐environmental sensitivity. Selection for BW was carried out for 13 generations in Malaysia. Subsets of 60 full‐sib families from Malaysia were sent to China and India after five and nine generations respectively. First, a multi‐trait animal model was used to analyse the BW in different countries as different traits. The results indicate a strong G × E. Second, a genetically structured environmental variance model, implemented using Bayesian inference, was used to analyse micro‐environmental sensitivity of BW in each country. The analysis revealed the presence of genetic heterogeneity of both BW and its environmental variance in all environments. The presence of genetic variation in residual variance of BW implies that the residual variance can be modified by selection. Incorporating both G × E and micro‐environmental sensitivity information may help in selecting robust genotypes with high performance across environments and resilience to environmental fluctuations.

Development of Diversity Arrays Technology markers as a tool for rapid genomic assessment in Nile tilapia, Oreochromis niloticus

The development of genomic markers is described for Nile tilapia, Oreochromis niloticus, using the Diversity Arrays Technology (DArT) genotype-by-sequencing platform. A total of 13 215 single nucleotide polymorphism (SNP) markers and 12 490 silicoDArT (dominant) markers were identified from broodstock of two selective breeding programs [Genetically Improved Farmed Tilapia (GIFT) strain from Malaysia and the Abbassa strain from Egypt]. Over 10 000 SNPs were polymorphic in either strain, and 2985 and 3087 showed strain-specific polymorphisms for the GIFT and Abbassa strains respectively. We demonstrate the potential utility of these markers for rapid genomic screening and use in breeding programs.

MHC-mediated spatial distribution in brown trout (Salmo trutta) fry

Major histocompatibility complex (MHC) class I-linked microsatellite data and parental assignment data for a group of wild brown trout (Salmo trutta L.) provide evidence of closer spatial aggregation among fry sharing greater numbers of MHC class I alleles under natural conditions. This result confirms predictions from laboratory experiments demonstrating a hierarchical preference for association of fry sharing MHC alleles. Full-siblings emerge from the same nest (redd), and a passive kin association pattern arising from limited dispersal from the nest (redd effect) would predict that all such pairs would have a similar distribution. However, this study demonstrates a strong, significant trend for reduced distance between pairs of full-sibling fry sharing more MHC class I alleles reflecting their closer aggregation (no alleles shared, 311.5 ± (s.e.)21.03 m; one allele shared, 222.2 ± 14.49 m; two alleles shared, 124.9 ± 23.88 m; P<0.0001). A significant trend for closer aggregation among fry sharing more MHC class I alleles was also observed in fry pairs, which were known to have different mothers and were otherwise unrelated (ML-r = 0) (no alleles: 457.6 ± 3.58 m; one allele (422.4 ± 3.86 m); two alleles (381.7 ± 10.72 m); P<0.0001). These pairs are expected to have emerged from different redds and a passive association would then be unlikely. These data suggest that sharing MHC class I alleles has a role in maintaining kin association among full-siblings after emergence. This study demonstrates a pattern consistent with MHC-mediated kin association in the wild for the first time.

Gene flow and population history in high dispersal marine invertebrates: mitochondrial DNA analysis of Holothuria nobilis (Echinodermata: Holothuroidea) populations from the Indo-Pacific

The sea cucumber, Holothuria nobilis, has a long-lived planktotrophic larvae, and previous allozyme surveys have suggested that high dispersal is realized. In contrast, recent ecological studies indicate that dispersal is low. To reconcile these data, and to investigate the evolution of this Indo-Pacific species, we screened geographical variation in 559 bp of a mitochondrial gene (COI) in 360 samples from the Australasian region and La Réunion. Sequences from La Réunion differed by > 7% from others and may constitute another species. Haplotype diversity in other samples was high (0.942, SD = 0.007), but haplotypes were closely related (mean nucleotide diversity: 0.0075, SD = 0.0041). AMOVA, pairwise FST values and exact tests did not detect significant population structure. Nested clade analysis showed that one of two main clades was over-represented in west Australia, whereas the other was more common in the northern Great Barrier Reef. Isolation-by-distance was identified as the main determinant of population structure at several clade levels. Contiguous range expansion was inferred for evolutionary older clade levels and this may correspond to a late Pleistocene (88 000-193 000 years ago) population expansion inferred from haplotype mismatch distributions. Thus, the population genetic structures detected are likely to be formed prior to the last ice age, with some indications for high dispersal on shorter time scales.

Mitochondrial DNA variation in Indo-Pacific populations of the giant tiger prawn, Penaeus monodon

Surveys of mitochondrial DNA (mtDNA) variation in the giant tiger prawn, Penaeus monodon, using restriction fragment length polymorphisms have provided the first clear evidence that the Indo-West Pacific region is a site of accumulation of genetic diversity rather than a site of origin of genetic diversity. No haplotyes were found in common between a group of five southeast African populations and a group of five Australian (including Western Australia) and three southeast Asian populations. The dominant haplotype was different in the Australian and southeast Asian population groups. Genetic diversity (pi) was greatest in Indonesia (pi averaged 0.05), less in the Philippines and Australia (pi averaged 0.01), and markedly less in the southeast African and the West Australian populations (pi averaged 0.003). The high diversity of the southeast Asian populations resulted from the occurrence in those populations of a set of haplotypes found only in southeast Asia but derived from the southeast African haplotypes. These genetic variants therefore evolved in the Indian Ocean and later migrated into the Indo-West Pacific region. Low genetic variation in the geographically marginal populations in southeast Africa and Western Australia is considered to be the result of bottlenecks, but mismatch distributions suggest that large population sizes have been maintained in Indonesian populations for long periods.

Small-Scale Dispersion of Eggs and Sperm of the Crown-of-Thorns Starfish (Acanthaster planci) in a Shallow Coral Reef Habitat

The dispersal of eggs and sperm of crown-of-thorns starfish, Acanthaster planci (L.), was measured in the field using an array of collectors up to 10 m downstream of a spawning starfish. Hydrodynamic measurements, gamete dispersal numerical models, and the gamete cloud dispersal measurements for the first time quantified the relationship between hydrodynamic conditions and the dispersion of eggs and sperm in the field. In general, gamete concentrations fell rapidly and logarithmically with distance from the spawning starfish; egg concentrations at 3 m were 1% of those near the starfish. Simplified dispersal models showed a good correspondence with these field data, and confirmed the observation that eggs rose higher in the water column and spread more laterally at low current speeds over the short spatial scales being considered. Fertilization rates, predicted from laboratory measurements of fertilization success and the gamete concentrations measured in the field, were estimated to be 90-100% within 1 m and 70-100% at 10 m. These results are explained by high success rates of fertilization (fertilizing capacity) at the measured dilutions, and were similar to fertilization rates previously measured by others for crown-of-thorns starfish in the field. Although the eggs were observed to spread upwards into the water column due to turbulence, laboratory measurements of sinking rates showed eggs to be very slightly negatively buoyant (median fall velocity of 0.072 mm·s-1), whereas sperm were neutrally buoyant. A significant fraction of eggs also entered the seabed near the starfish; the proportion decreased with increasing current strength. This process may provide a mechanism for enhanced fertilization of these gametes and/or a mechanism for self-recruitment to a given reef population.

The Effects of Sperm Concentration, Sperm:Egg Ratio, and Gamete Age on Fertilization Success in Crown-of-Thorns Starfish (Acanthaster planci) in the Laboratory

Laboratory experiments varying gamete concentrations and gamete age demonstrated significant reductions in fertilization success of the starfish Acanthaster planci (L.) with decreasing sperm concentration and increasing age of both eggs and sperm. The effect of aging in sperm was faster than that of eggs, and the speed of sperm aging increased with increasing dilution of sperm. Fertilization success was high over a wide range of sperm: egg ratios but declined rapidly at ratios less than 50, particularly at low sperm concentrations. A. planci gametes aged more slowly, and the loss of fertilizing capacity of sperm with dilution (the respiratory dilution effect) was far less, than in sea urchins. These characteristics provide a mechanism for enhanced fertilization success at given sperm concentrations and at greater distances and times from the point of gamete release, and may explain the higher fertilization rates achieved over longer distances in the wild by A. planci relative to sea urchins. Gametes would remain competent for longer periods at more dilute concentrations and so better achieve long-distance fertilization. Gametes obtained at the end of the breeding season were qualitatively different from those obtained early in the breeding season and showed reduced fertilization success for a given combination of variables, and different fertilization dynamics.