Allelic expression changes in Medaka (Oryzias latipes) hybrids between inbred strains derived from genetically distant populations

Variations in allele expressions between genetically distant populations are one of the most important factors which affects their morphological and physiological variations. These variations are caused by natural mutations accumulated in their habitats. It has been reported that allelic expression differences in the hybrids of genetically distant populations are different from parental strains. In that case, there is a possibility that allelic expression changes lead to novel phenotypes in hybrids. Based on genomic information of the genetically distant populations, quantification and comparison of allelic expression changes make importance of regulatory sequences (cis-acting factors) or upstream regulatory factors (trans-acting modulators) for these changes clearer. In this study, we focused on two Medaka inbred strains, Hd-rR and HNI, derived from genetically distant populations and their hybrids. They are highly polymorphic and we can utilize whole-genome information. To analyze allelic expression changes, we established a method to quantify and compare allele-specific expressions of 11 genes between the parental strains and their reciprocal hybrids. In intestines of reciprocal hybrids, allelic expression was either similar or different in comparison with the parental strains. Total expressions in Hd-rR and HNI were tissue-dependent in the case of HPRT1, with high up-regulation of Hd-rR allele expression in liver. The proportion of genes with differential allelic expression in Medaka hybrids seems to be the same as that in other animals, despite the high SNP rate in the genomes of the two inbred strains. It is suggested that each tissue of the strain difference in trans-acting modulators is more important than polymorphisms in cis-regulatory sequences in producing the allelic expression changes in reciprocal hybrids.     

Patterns of phenotypic covariation and correlation in modern humans as viewed from morphological integration

Proportionality of phenotypic and genetic distance is of crucial importance to adequately focus on population history and structure, and it depends on the proportionality of genetic and phenotypic covariance. Constancy of phenotypic covariances is unlikely without constancy of genetic covariation if the latter is a substantial component of the former. If phenotypic patterns are found to be relatively stable, the most probable explanation is that genetic covariance matrices are also stable. Factors like morphological integration account for such stability. Morphological integration can be studied by analyzing the relationships among morphological traits. We present here a comparison of phenotypic correlation and covariance structure among worldwide human populations. Correlation and covariance matrices between 47 cranial traits were obtained for 28 populations, and compared with design matrices representing functional and developmental constraints. Among-population differences in patterns of correlation and covariation were tested for association with matrices of genetic distances (obtained after an examination of 10 Alu-insertions) and with Mahalanobis distances (computed after craniometrical traits). All matrix correlations were estimated by means of Mantel tests. Results indicate that correlation and covariance structure in our species is stable, and that among-group correlation/covariance similarity is not related to genetic or phenotypic distance. Conversely, genetic and morphological distance matrices were highly correlated. Correlation and covariation patterns were largely associated with functional and developmental factors, which probably account for the stability of covariance patterns.     

Discriminant functions for linear comparisons of means when covariance matrices are unequal

Situations exist, as in the biological example of discriminant analysis for natural hybridization, cited in the text, where (a) not all populations have equal variances, and (b) comparisions based on single degrees of freedom must be planned. This paper presents a statistical methodology of estimating discriminant functions for linear comparisons among k(<2) multivariate normal populations, and of testing their significance, when these populations have unequal covariance matrices.     

A universally applicable method of isolating somatic cell hybrids by two-colour flow sorting

A generally applicable technique is described that permits easy identification and isolation of heterokaryons a few hours after fusion. It is based on the labelling of living cells with different fluorochromes, which, at appropriate concentrations do not affect viability or gene expression. Both fluorochromes are relatively stable and do not cross-contaminate unlabelled cells. The technique has a powerful potential in studies on gene regulation in somatic cell hybrids since heterofluorescent hybrids between any type of cells can be isolated directly from large populations of monofluorescent parental cells by using a cell sorter equipped with a single laser. Thus the technique avoids the need for genetically marked parental cells for selection.     

Hybridization between tetraploid species of Alopecurus L. (Poaceae): morphological studies of natural and artificial hybrids

The morphology of individuals of mixed natural populations of three species of Alopecurus, A. pratensis, A. geniculars and A. arundinaceus , is described. It is suggested that plants morphologically intermediate between the parental species are hybrids or of hybrid origin. Hybridization occurs more readily between A. pratensis and A. arundinaceus than between A. pratensis and A. geniculate , since intermediate forms were recovered from all mixed populations of A. pratensis and A. arundinaceus. The intermediate forms occur only in a few populations of A. pratensis and A. geniculars where an appropriate habitat is available. Evidence of backcrossing and therefore introgression comes from the continuum of intermediate forms that exists between the species. Scatter diagrams for individual populations show distinct differences in the degree of introgression. Observations on artificially produced hybrids are similar to those on naturally occurring hybrids.     

Hybridization in coral reef fishes: introgression and bi-directional gene exchange in Thalassoma (family Labridae)

Hybrids in coral reef fishes have traditionally been described based on external features using meristic characters and colouration to identify putative parental contributors. This study utilised molecular genetic techniques to verify hybrid status and identify putative parental species for five hybrid specimens (Labridae: Thalassoma) collected from Holmes Reef in the Coral Sea. Phylogenetic analyses support hybrid origins of the specimens. Mitochondrial COI gene, nuclear S7 (intron 1) and nuclear copy of mitochondrial (NUMT) D-loop region corroborate the identity of T. quinquevittatum as the maternal and T. jansenii as the paternal contributor. Backcrossing to parental species by hybrids and bi-directional gene exchange between the Holmes Reef populations of T. jansenii and T. quinquevittatum was detected, suggesting that hybrids are fertile and able to reproduce successfully. F(1) hybrids display a mixture of the colouration attributes of the two parental species, but subsequent backcrossed individuals were unrecognisable as hybrids and displayed colouration of either parental species. A large numerical imbalance exists between the putative parental species at Holmes Reef, with T. quinquevittatum outnumbering T. jansenii by approximately 25:1. In this case study, hybridization appears to be driven by ecological rather than evolutionary factors.     

The evolution of the phenotypic covariance matrix: evidence for selection and drift in Melanoplus

Phenotypic variation in trait means is a common observation for geographically separated populations. Such variation is typically retained under common garden conditions, indicating that there has been evolutionary change in the populations, as a result of selection and/or drift. Much less frequently studied is variation in the phenotypic covariance matrix (hereafter, P matrix), although this is an important component of evolutionary change. In this paper, we examine variation in the phenotypic means and P matrices in two species of grasshopper, Melanoplus sanguinipes and M. devastator. Using the P matrices estimated for 14 populations of M. sanguinipes and three populations of M. devastator we find that (1) significant differences between the sexes can be attributed to scaling effects; (2) there is no significant difference between the two species; (3) there are highly significant differences among populations that cannot be accounted for by scaling effects; (4) these differences are a consequence of statistically significant patterns of covariation with geographic and environmental factors, phenotypic variances and covariances increasing with increased temperature but decreasing with increased latitude and altitude. This covariation suggests that selection has been important in the evolution of the P matrix in these populations Finally, we find a significant positive correlation between the average difference between matrices and the genetic distance between the populations, indicating that drift has caused some of the variation in the P matrices.     

Contributions of heterosis and epistasis to hybrid fitness

Early-generation hybrid fitness is difficult to interpret because heterosis can obscure the effects of hybrid breakdown. We used controlled reciprocal crosses and common garden experiments to distinguish between effects of heterosis and nuclear and cytonuclear epistasis among morphotypes and advanced-generation hybrid derivative populations in the Piriqueta caroliniana (Turneraceae) plant complex. Seed germination, growth, and sexual reproduction of first-generation hybrids, inbred parental lines, and outbred parental lines were compared under field conditions. Average vegetative performance was greater for hybrids than for inbred lines, and first-season growth was similar for hybrids and outbred parental lines. Hybrid survival surpassed that of inbred lines and was equal to or greater than outbred lines' survival, and more F(1) than parental plants reproduced. Reductions in hybrid fitness due to Dobzhansky-Muller incompatibilities (epistasis among divergent genetic elements) were expressed as differences in vegetative growth, survival, and reproduction between plants from reciprocal crosses for both F(1) and backcross hybrid generations. Comparing performance of hybrids against parental genotypes from intra- and interpopulation crosses allowed a more robust prediction of F(1) hybrids' success and more accurate interpretations of the genetic architecture of F(1) hybrid vigor.     

Expression of anatomical leaf traits in homoploid hybrids between deciduous and evergreen species of Vaccinium

We investigated the anatomical expression of leaf traits in hybrids between evergreen Vaccinium vitis‐idaea and deciduous V. myrtillus. We compared parents from four populations with their respective F1 hybrids and tested whether (i) transgression can be the source of novel anatomical traits in hybrids; (ii) expression of transgressive traits is more probable for traits with similar values in parents and intermediate for more distinct values, as predicted by theory; and (iii) independent origin of hybrids leads to identical trait expression profiles among populations. We found that anatomical leaf traits can be divided into four categories based on their similarity to parents: intermediate, parental‐like, transgressive and non‐significant. Contrary to the common view, parental‐like trait values were equally important in shaping the hybrid profile, as were intermediate traits. Transgression was revealed in 17/144 cases and concerned mainly cell and tissue sizes. As predicted by theory, we observed transgressive segregation more often when there was little phenotypic divergence, but intermediate values when parental traits were differentiated. It is likely that cell and tissue sizes are phylogenetically more conserved due to stabilising selection, whereas traits such as leaf thickness and volume fraction of the intercellular spaces, showing a consistent intermediate pattern across populations, are more susceptible to directional selection. Hybrid populations showed little similarity in expression profile, with only three traits identically expressed across all populations. Thus local adaptation of parental species and specific genetic background may be of importance.     

Comparison of Populations Based on Attribute Data,when Covariance Matrices are Unequal

Mahalanobis—like distance indices B² and G² of BALAKRISHNAN and SANGHVI (1968) are employed to compare multinomial populations of proportions in human data. Recently, BOWERING and MISRA (1982) employed these methods, with modifications, in a different field viz. to compare fish stocks, based on meristic characters. The method of BALAKRISHNAN and SANGHVI (1968), however, assumes equality of covariance matrices, which is doubtful with multinomial populations of proportions. In this paper a statistical method for comparing such populations (where covariance matrices are not equal) is presented. Also, a feature of the G² method of BALAKRISHNAN and SANGHVI (1968) is discussed, which makes its statistical status a little more clear.     

A note on simulating null distributions for G matrix comparisons

Genetic variances and covariances, summarized in G matrices, are key determinants of the course of adaptive evolution. Consequently, understanding how G matrices vary among populations is critical to answering a variety of questions in evolutionary biology. A method has recently been proposed for generating null distributions of statistics pertaining to differences in G matrices among populations. The general approach facilitated by this method is likely to prove to be very important in studies of the evolution of G . We have identified an issue in the method that will cause it to create null distributions of differences in G matrices that are likely to be far too narrow. The issue arises from the fact that the method as currently used generates null distributions of statistics pertaining to differences in G matrices across populations by simulating breeding value vectors based on G matrices estimated from data, randomizing these vectors across populations, and then calculating null values of statistics from G matrices that are calculated directly from the variances and covariances among randomized vectors. This calculation treats breeding values as quantities that are directly measurable, instead of predicted from G matrices that are themselves estimated from patterns of covariance among kin. The existing method thus neglects a major source of uncertainty in G matrices, which renders it anti‐conservative. We first suggest a correction to the method. We then apply the original and modified methods to a very simple instructive scenario. Finally, we demonstrate the use of both methods in the analysis of a real data set.     

The recombination landscape in Arabidopsis thaliana F2 populations

Recombination during meiosis shapes the complement of alleles segregating in the progeny of hybrids, and has important consequences for phenotypic variation. We examined allele frequencies, as well as crossover (XO) locations and frequencies in over 7000 plants from 17 F(2) populations derived from crosses between 18 Arabidopsis thaliana accessions. We observed segregation distortion between parental alleles in over half of our populations. The potential causes of distortion include variation in seed dormancy and lethal epistatic interactions. Such a high occurrence of distortion was only detected here because of the large sample size of each population and the number of populations characterized. Most plants carry only one or two XOs per chromosome pair, and therefore inherit very large, non-recombined genomic fragments from each parent. Recombination frequencies vary between populations but consistently increase adjacent to the centromeres. Importantly, recombination rates do not correlate with whole-genome sequence differences between parental accessions, suggesting that sequence diversity within A. thaliana does not normally reach levels that are high enough to exert a major influence on the formation of XOs. A global knowledge of the patterns of recombination in F(2) populations is crucial to better understand the segregation of phenotypic traits in hybrids, in the laboratory or in the wild.     

Daphnia Galeata × D. Longispina Hybrids in Western Norway

We describe the occurrence of D. galeata×longispina hybrids in two lakes of western Norway. Parental species and interspecific hybrids were characterised by both nuclear and mitochondrial molecular markers. In one of the populations, hybrids were shown to dominate the population over several years. A few individuals in both populations were probably not F1 hybrids, but possibly backcrosses or F2 hybrids. Most (possibly all) F1 hybrids were of D. galeata maternal origin. In addition, interspecific hybrids could be identified based on morphological characters, which were intermediate between the parental species. Interspecific hybridisation between these two species is remarkable, since they are distantly related.     

Hemiclone diversity in the hybridogenetic frog Rana esculenta outside the area of clone formation: the view from protein electrophoresis

European water frog hybrids Rana esculenta reproduce hemiclonally, by hybridogenesis: In the germ line they exclude the genome of the parental species Rana lessonae and produce haploid, unrecombined gametes with a genome of the parental species Rana ridibunda . These hybrids coexist with and depend as sexual parasites on the host parental species R. lessonae (the L-E population system); matings with R. lessonae restore somatic hybridity in each generation of R. esculenta . We investigated 15 L-E system populations in northern Switzerland, which is outside R. ridibunda 's native range. Frequency of hybrids in samples varied from 8% in marsh ponds to 100% in gravel pits and forest ponds. Clonal diversity (variation among R. ridibunda genomes of hybrids), detected by six protein electrophoretic marker loci, revealed a total of eight hemiclones and locally ranged from uniclonal populations in southern parts of the survey region to six coexisting hemiclones in the north. All alleles distinguishing hemiclones occur commonly in the nearest native R. ridibunda populations of east-central Europe; the most probable source of clonal diversity in our samples is multiple clone formation by primary hybridizations in the sympatry area of R. ridibunda and R. lessonae and subsequent dispersal of hemiclonal lineages. A positive correlation between amount of clonal diversity and hybrid frequency, predicted by the Frozen Niche Variation (FNV) model (each hemiclone is characterized by a relatively narrow niche, coexistence is possible through niche partitioning), was not found; this contrasts with hemiclonally reproducing fish hybrids ( Poeciliopsis ). Historical factors, such as availability of different colonizing hemiclones may be strong enough to override the signal from operation of the FNV.     

Impact of Hybridization on Shape,Variation and Covariation of the Mouse Molar

Hybridization can generate phenotypes departing from the parental range through many processes: increased variance due to increased heterozygosity, decreased phenotypic integration and altered covariation structure due to changes in the expression of (co)variation generating developmental processes, and transgression that may arise from epistasis and compensatory genes. Morphometric assessment of shape differentiation, variance and covariance may shed light on these processes, especially for complex characters for which the genetic basis has not yet been assessed. The shape of the first upper molar was thus quantified in a cross between inbred strains of the two European subspecies of the house mouse (Mus musculus domesticus and M. m. musculus). Hybrids appeared as moderately transgressive. Morphological variance was increased in F2s, back to levels observed in wild populations. The pattern of variance–covariance was different between the two parental strains, but restored to a wild-type structure in F2s. Finally, F2s displayed a degree of morphological integration comparable to wild populations but lower than observed in the parental strains. This is interpreted as the result of the reshuffling of the standing genetic variation in hybrids that should have restored the expression of (co)variance generating processes made ineffective in parents due to random allele fixation of some loci. Inter-parental differentiation was more important in regions of the tooth developing early during embryogenesis, whereas transgression was more pronounced in late developing regions. Mutations may more easily generate important geometric differences early during tooth development, occurring as a cascade of signalization from the first cusp to initiate onward. Epistasis and constraints of neighboring parts may be more important later, promoting transgression and impeding differentiation.     

Developmental instability in hybrids between Lychnis viscaria and Lychnis alpina (Caryophyllaceae)

Developmental instability shown by increased fluctuating asymmetry can be caused by either genetic or environmental stress. Genomic coadaptation and heterozygosity are the genetic factors that are commonly assumed to increase the level of developmental stability. Therefore, in hybrid populations the level of fluctuating asymmetry (FA) can be lower due to higher heterozygosity or higher due to disruption of coadapted gene complexes, depending on the degree of divergence between hybridizing taxa. Here I present data on FA in petals from hybrids between Lychnis viscaria (Caryophyllaceae) and Lychnis alpina and from parental species grown in a common garden environment. Petal asymmetry of hybrids was clearly higher than that of either parental species grown in common environment. Between the two parental species petal asymmetry did not differ. The mean size of the petals in hybrids was about the same as in L. viscaria, thus indicating no heterotic effect. Therefore, it seems that hybrids between L. viscaria and L. alpina do suffer from the disruption of coadapted gene complexes as indicated by higher developmental instability.     

Genomic prediction of hybrid performance in maize with models incorporating dominance and population specific marker effects

Identifying high performing hybrids is an essential part of every maize breeding program. Genomic prediction of maize hybrid performance allows to identify promising hybrids, when they themselves or other hybrids produced from their parents were not tested in field trials. Using simulations, we investigated the effects of marker density (10, 1, 0.3 marker per mega base pair, Mbp(-1)), convergent or divergent parental populations, number of parents tested in other combinations (2, 1, 0), genetic model (including population-specific and/or dominance marker effects or not), and estimation method (GBLUP or BayesB) on the prediction accuracy. We based our simulations on marker genotypes of Central European flint and dent inbred lines from an ongoing maize breeding program. To simulate convergent or divergent parent populations, we generated phenotypes by assigning QTL to markers with similar or very different allele frequencies in both pools, respectively. Prediction accuracies increased with marker density and number of parents tested and were higher under divergent compared with convergent parental populations. Modeling marker effects as population-specific slightly improved prediction accuracy under lower marker densities (1 and 0.3?Mbp(-1)). This indicated that modeling marker effects as population-specific will be most beneficial under low linkage disequilibrium. Incorporating dominance effects improved prediction accuracies considerably for convergent parent populations, where dominance results in major contributions of SCA effects to the genetic variance among inter-population hybrids. While the general trends regarding the effects of the aforementioned influence factors on prediction accuracy were similar for GBLUP and BayesB, the latter method produced significantly higher accuracies for models incorporating dominance.     

Hybrid bridges to gene flow: a case study in milkweeds (Asclepias)

Natural hybridization occurs throughout areas of sympatry for the North American milkweeds Asclepias exaltata and A. syriaca (Asclepiadaceae), even though the formation of F1 hybrid seed is a rare event. For introgressive hybridization to proceed, F1 and advanced hybrids must be released from reproductive barriers and successfully mate with one or both parental species. I investigated the mating system of natural hybrids between A. exaltata and A. syriaca in three populations in Shenandoah National Park, Virginia. Allozyme data and a maximum-likelihood procedure were used to estimate the frequency of six genotypic classes (parentals, F1, F2, and backcrosses) of the hybridizing populations, the pollinia received by hybrid plants, and the paternal parents of seeds produced by hybrids. F1 hybrids, backcross A. syriaca, and parental A. syriaca individuals were common in three hybrid populations. Even though self-pollinations and interhybrid pollinations were common, F2 seed production and the occurrence of F2 individuals were rare in hybrid populations. Hybrid plants received more pollen from A. syriaca than A. exaltata, which resulted in the production of more backcross-A. syriaca seed than backcross-A. exaltata seed. Asclepias exaltata was rare in the hybrid populations, but A. exaltata pollinia were received by hybrids and this species sired between 15% and 36% of the seeds produced on hybrids. The potential for introgression with A. exaltata populations is lower because this species is unsuccessful as the maternal parent in interspecific and backcross hand-pollinations. The asymetry of hybridization with A. syriaca as the maternal parent is further supported by the incorporation of maternally inherited chloroplast DNA markers in hybrids. Hybrid milkweeds frequently backcross with both parental species and may be released from the reproductive barriers that limit the formation of F1 hybrids in natural populations. The direction of interspecific gene flow and introgression in milkweeds is influenced by the reproductive biology of hybrids, the constituency of the surrounding population, and failure of some crosses to produce seeds. Finally, introgressive hybridization remains an important evolutionary force even when the initial formation of F1 hybrids in natural populations is rare.     

A metabolome‐based core hybridisation strategy for the prediction of rice grain weight across environments

Marker‐based prediction holds great promise for improving current plant and animal breeding efficiencies. However, the predictabilities of complex traits are always severely affected by negative factors, including distant relatedness, environmental discrepancies, unknown population structures, and indeterminate numbers of predictive variables. In this study, we utilised two independent F₁ hybrid populations in the years 2012 and 2015 to predict rice thousand grain weight (TGW) using parental untargeted metabolite profiles with a partial least squares regression method. A stable predictive model for TGW was built based on hybrids from the population in 2012 (r = 0.75) but failed to properly predict TGW for hybrids from the population in 2015 (r = 0.27). After integrating hybrids from both populations into the training set, the TGW of hybrids could be predicted but was largely dependent on population structures. Then, core hybrids from each population were determined by principal component analysis and the TGW of hybrids in both environments were successfully predicted (r > 0.60). Moreover, adjusting the population structures and numbers of predictive analytes increased TGW predictability for hybrids in 2015 (r = 0.72). Our study demonstrates that the TGW of F₁ hybrids across environments can be accurately predicted based on parental untargeted metabolite profiles with a core hybridisation strategy in rice. Metabolic biomarkers identified from early developmental stage tissues, which are grown under experimental conditions, may represent a workable approach towards the robust prediction of major agronomic traits for climate‐adaptive varieties.