Transcriptome-wide selection and validation of a solid set of reference genes for gene expression studies in the cephalopod mollusk Octopus vulgaris.

Octopus vulgaris is a cephalopod mollusk and an active marine predator that has been at the center of a number of studies focused on the understanding of neural and biological plasticity. Studies on the machinery involved in e.g., learning and memory, regeneration, and neuromodulation are required to shed light on the conserved and/or unique mechanisms that these animals have evolved. Analysis of gene expression is one of the most essential means to expand our understanding of biological machinery, and the selection of an appropriate set of reference genes is the prerequisite for the quantitative real-time polymerase chain reaction (qRT-PCR). Here we selected 77 candidate reference genes (RGs) from a pool of stable and relatively high-expressed transcripts identified from the full-length transcriptome of O. vulgaris, and we evaluated their expression stabilities in different tissues through geNorm, NormFinder, Bestkeeper, Delta-CT method, and RefFinder. Although various algorithms provided different assemblages of the most stable reference genes for the different kinds of tissues tested here, a comprehensive ranking revealed RGs specific to the nervous system (Ov-RNF7 and Ov-RIOK2) and Ov-EIF2A and Ov-CUL1 across all considered tissues. Furthermore, we validated RGs by assessing the expression profiles of nine target genes (Ov-Naa15, Ov-Ltv1, Ov-CG9286, Ov-EIF3M, Ov-NOB1, Ov-CSDE1, Ov-Abi2, Ov-Homer2, and Ov-Snx20) in different areas of the octopus nervous system (gastric ganglion, as control). Our study allowed us to identify the most extensive set of stable reference genes currently available for the nervous system and appendages of adult O. vulgaris.

Genetic load: genomic estimates and applications in non-model animals.

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Ancient genomes reveal early Andean farmers selected common beans while preserving diversity.

All crops are the product of a domestication process that started less than 12,000 years ago from one or more wild populations1,2. Farmers selected desirable phenotypic traits (such as improved energy accumulation, palatability of seeds and reduced natural shattering3) while leading domesticated populations through several more or less gradual demographic contractions2,4. As a consequence, the erosion of wild genetic variation5 is typical of modern cultivars, making them highly susceptible to pathogens, pests and environmental change6,7. The loss of genetic diversity hampers further crop improvement programmes to increase food production in a changing world, posing serious threats to food security8,9. Using both ancient and modern seeds, we analysed the temporal dynamics of genetic variation and selection during the domestication process of the common bean (Phaseolus vulgaris) in the southern Andes. Here, we show that most domestic traits were selected for before 2,500 years ago, with no or only minor loss of whole-genome heterozygosity. In fact, most of the changes at coding genes and linked regions that differentiate wild and domestic genomes are already present in the ancient genomes analysed here, and all ancient domestic genomes dated between 600 and 2,500 years ago are highly variable (at least as variable as modern genomes from the wild). Single seeds from modern cultivars show reduced variation when compared with ancient seeds, indicating that intensive selection within cultivars in the past few centuries probably partitioned ancestral variation within different genetically homogenous cultivars. When cultivars from different Andean regions are pooled, the genomic variation of the pool is higher than that observed in the pool of ancient seeds from north and central western Argentina. Considering that most desirable phenotypic traits are probably controlled by multiple polymorphic genes10, a plausible explanation of this decoupling of selection and genetic erosion is that early farmers applied a relatively weak selection pressure2 by using many phenotypically similar but genetically diverse individuals as parents. Our results imply that selection strategies during the past few centuries, as compared with earlier times, more intensively reduced genetic variation within cultivars and produced further improvements by focusing on a few plants carrying the traits of interest, at the cost of marked genetic erosion within Andean landraces.

Diving into divergence: Differentiation in swimming performances, physiology and gene expression between locally-adapted sympatric cichlid fishes.

Sympatric speciation occurs without geographical barriers and is thought to often be driven by ecological specialization of individuals that eventually diverge genetically and phenotypically. Distinct morphologies between sympatric populations occupying different niches have been interpreted as such differentiating adaptive phenotypes, yet differences in performance and thus likely adaptiveness between them were rarely tested. Here, we investigated if divergent body shapes of two sympatric crater lake cichlid species from Nicaragua, one being a shore-associated (benthic) species while the other prefers the open water zones (limnetic), affect cruising (Ucrit ) and sprinting (Usprint ) swimming abilities - performances particularly relevant to their respective lifestyles. Furthermore, we investigated species differences in oxygen consumption (MO2 ) across different swimming speeds and compare gene expression in gills and white muscle at rest and during exercise. We found a superior cruising ability in the limnetic Amphilophus zaliosus compared to the benthic Amphilophus astorquii, while sprinting was not different, suggesting that their distinct morphologies affect swimming performance. Increased cruising swimming ability in A. zaliosus was linked to a higher oxygen demand during activity (but not rest), indicating different metabolic rates during exercise - a hypothesis supported by coinciding gene expression patterns of gill transcriptomes. We identified differentially expressed genes linked to swimming physiology, regulation of swimming behaviour and oxygen intake. A combination of physiological and morphological differences may thus underlie adaptations to these species' distinct niches. This complex ecological specialization probably resulted in morphological and physiological trade-offs that contributed to the rapid establishment and maintenance of divergence with gene flow.

Gene(s) and individual feeding behavior: Exploring eco-evolutionary dynamics underlying left-right asymmetry in the scale-eating cichlid fish Perissodus microlepis.

The scale-eating cichlid fish Perissodus microlepis is a textbook example of bilateral asymmetry due to its left or right-bending heads and of negative frequency-dependent selection, which is proposed to maintain this stable polymorphism. The mechanisms that underlie this asymmetry remain elusive. Several studies had initially postulated a simple genetic basis for this trait, but this explanation has been questioned, particularly by reports observing a unimodal distribution of mouth shapes. We hypothesize that this unimodal distribution might be due to a combination of genetic and phenotypically plastic components. Here, we expanded on previous work by investigating a formerly identified candidate SNP associated to mouth laterality, documenting inter-individual variation in feeding preference using stable isotope analyses, and testing their association with mouth asymmetry. Our results suggest that this polymorphism is influenced by both a polygenic basis and inter-individual non-genetic variation, possibly due to feeding experience, individual specialization, and intraspecific competition. We introduce a hypothesis potentially explaining the simultaneous maintenance of left, right, asymmetric and symmetric mouth phenotypes due to the interaction between diverse eco-evolutionary dynamics including niche construction and balancing selection. Future studies will have to further tease apart the relative contribution of genetic and environmental factors and their interactions in an integrated fashion.

Towards understanding the genetic basis of mouth asymmetry in the scale-eating cichlid Perissodus microlepis.

How polymorphisms consisting in left-right asymmetries are produced and maintained in natural populations is a tantalizing question, which remains largely unanswered. The scale-eating cichlid fish Perissodus microlepis is a remarkable example of extreme ecological specialization achieved by morphological and behavioural laterality. Its asymmetric mouth is accompanied by a pronounced lateralized foraging behaviour, where a left-bending morph preferentially feeds on the scales of the right side of its prey, while the opposite is true for the right morph. This striking asymmetry made this fish a textbook example of the astounding degree of ecological specialization and negative frequency-dependent selection. Yet, the genetic basis underlying this spectacular laterality remains unknown. We addressed this question through analyses of wild-caught fish using high-throughput DNA sequencing data. A novel array of SNP markers was developed by ddRAD sequencing (ddRADseq) and the use of pooled DNA samples (PoolSeq). We obtained more than 155 000 SNPs using ddRADseq and 3 900 000 SNPs with PoolSeq. Among these, we identified one (ddRAD) SNP, and 38 or 378 (PoolSeq) windows that are differentiated between the left and right morphs accounting for spurious associations due to geographic structuring. This allowed us to uncover candidate genomic regions that potentially contain genes for this trait. Then, this interesting trait has a genetic basis that is likely to be influenced by multiple loci. This result contributes to a greater understanding of the genetic bases of left-right asymmetry and, ultimately, the evolutionary processes governing the maintenance of this striking case of laterality.

Are sympatrically speciating Midas cichlid fish special? Patterns of morphological and genetic variation in the closely related species Archocentrus centrarchus.

Established empirical cases of sympatric speciation are scarce, although there is an increasing consensus that sympatric speciation might be more common than previously thought. Midas cichlid fish are one of the few substantiated cases of sympatric speciation, and they formed repeated radiations in crater lakes. In contrast, in the same environment, such radiation patterns have not been observed in other species of cichlids and other families of fish. We analyze morphological and genetic variation in a cichlid species (Archocentrus centrarchus) that co-inhabits several crater lakes with the Midas species complex. In particular, we analyze variation in body and pharyngeal jaw shape (two ecologically important traits in sympatrically divergent Midas cichlids) and relate that to genetic variation in mitochondrial control region and microsatellites. Using these four datasets, we analyze variation between and within two Nicaraguan lakes: a crater lake where multiple Midas cichlids have been described and a lake where the source population lives. We do not observe any within-lake clustering consistent across morphological traits and genetic markers, suggesting the absence of sympatric divergence in A. centrarchus. Genetic differentiation between lakes was low and morphological divergence absent. Such morphological similarity between lakes is found not only in average morphology, but also when analyzing covariation between traits and degree of morphospace occupation. A combined analysis of the mitochondrial control region in A. centrarchus and Midas cichlids suggests that a difference between lineages in the timing of crater lake colonization cannot be invoked as an explanation for the difference in their levels of diversification. In light of our results, A. centrarchus represents the ideal candidate to study the genomic differences between these two lineages that might explain why some lineages are more likely to speciate and diverge in sympatry than others.