F2 Intercross Research Articles

Analysis of inheritance in a novel seedling stage zebra leaf mutant, Pusa Zebra 18 in rice (Oryza sativa L.)

A novel zebra leaf mutant, Pusa Zebra 18 derived spontaneously from F2 intercross population of Pusa 44 near-isogenic lines (NILs) exhibits distinguishable yellowish cross bands across the green leaf blade, specifically expressed in the seedling stage from 6 DAT, which diffuses progressively, leading to the recovery of the normal green leaf by 30 DAT. We elucidated the genetics and nature of inheritance for zebra trait in F2 and BC1F2 populations by crossing Pusa 6B with Pusa Zebra 18. The leaves of the F1 were a normal green color, indicating that the seedling stage zebra leaf trait in Pusa Zebra 18 is recessive. Chi-square goodness of fit analysis of these segregating populations showed that the segregation ratio was fitting to the Mendelian segregation ratio of 3:1. The seedling stage zebra leaf mutation in Pusa Zebra 18 is monogenic recessive in nature. Further, mapping, fine mapping, cloning, and functional characterization need to be carried out, which will help gain insight into the mechanism involved in chlorophyll biosynthesis and regulation studies in rice.

Genetic Analysis of Obesity-Induced Diabetic Nephropathy in BTBR Mice.

Diabetic nephropathy (DN) is the leading cause of end-stage renal disease in the United States and has a significant impact on human suffering. Leptin-deficient BTBR mice (BTBRob/ob) develop hallmark features of obesity-induced diabetic nephropathy, whereas leptin-deficient C57BL/6J (B6ob/ob) mice do not. To identify genetic loci that underlie this strain difference, we constructed an F2 intercross between BTBRob/ob and B6ob/ob mice. We isolated kidneys from 460 F2 mice and histologically scored them for percent mesangial matrix and glomerular volume (∼50 glomeruli per mouse), yielding ∼45,000 distinct measures in total. The same histological measurements were made in kidneys from B6 and BTBR mice, either lean or obese (Lepob/ob) at 4 and 10 weeks of age, allowing us to assess the contribution of strain, age, and obesity on glomerular pathology. All F2 mice were genotyped for ∼5,000 single nucleotide polymorphisms (SNPs), ∼2,000 of which were polymorphic between B6 and BTBR, enabling us to identify a quantitative trait locus (QTL) on chromosome 7, with a peak at ∼30 Mbp, for percent mesangial matrix, glomerular volume, and mesangial volume. The podocyte-specific gene, Nephrin (Nphs1) is physically located at the QTL and contains high-impact SNPs in BTBR, including several missense variants within the extracellular Ig-like domains.

Dominance is common in mammals and is associated with trans-acting gene expression and alternative splicing

BackgroundDominance and other non-additive genetic effects arise from the interaction between alleles, and historically these phenomena play a major role in quantitative genetics. However, most genome-wide association studies (GWAS) assume alleles act additively.ResultsWe systematically investigate both dominance—here representing any non-additive within-locus interaction—and additivity across 574 physiological and gene expression traits in three mammalian stocks: F2 intercross pigs, rat heterogeneous stock, and mice heterogeneous stock. Dominance accounts for about one quarter of heritable variance across all physiological traits in all species. Hematological and immunological traits exhibit the highest dominance variance, possibly reflecting balancing selection in response to pathogens. Although most quantitative trait loci (QTLs) are detectable as additive QTLs, we identify 154, 64, and 62 novel dominance QTLs in pigs, rats, and mice respectively that are undetectable as additive QTLs. Similarly, even though most cis-acting expression QTLs are additive, gene expression exhibits a large fraction of dominance variance, and trans-acting eQTLs are enriched for dominance. Genes causal for dominance physiological QTLs are less likely to be physically linked to their QTLs but instead act via trans-acting dominance eQTLs. In addition, thousands of eQTLs are associated with alternatively spliced isoforms with complex additive and dominant architectures in heterogeneous stock rats, suggesting a possible mechanism for dominance.ConclusionsAlthough heritability is predominantly additive, many mammalian genetic effects are dominant and likely arise through distinct mechanisms. It is therefore advantageous to consider both additive and dominance effects in GWAS to improve power and uncover causality.

C3H/HeNSlc mouse with low phospholipid transfer protein expression showed dyslipidemia

High serum levels of triglycerides (TG) and low levels of high-density lipoprotein cholesterol (HDL-C) increase the risk of coronary heart disease in humans. Herein, we first reported that the C3H/HeNSlc (C3H-S) mouse, a C3H/HeN-derived substrain, is a novel model for dyslipidemia. C3H-S showed hypertriglyceridemia and low total cholesterol (TC), HDL-C, and phospholipid (PL) concentrations. To identify the gene locus causing dyslipidemia in C3H-S, we performed genetic analysis. In F2 intercrosses between C3H-S mice and strains with normal serum lipids, the locus associated with serum lipids was identified as 163–168 Mb on chromosome 2. The phospholipid transfer protein (Pltp) gene was a candidate gene within this locus. Pltp expression and serum PLTP activity were markedly lower in C3H-S mice. Pltp expression was negatively correlated with serum TG and positively correlated with serum TC and HDL-C in F2 mice. Genome sequencing analysis revealed that an endogenous retrovirus (ERV) sequence called intracisternal A particle was inserted into intron 12 of Pltp in C3H-S. These results suggest that ERV insertion within Pltp causes aberrant splicing, leading to reduced Pltp expression in C3H-S. This study demonstrated the contribution of C3H-S to our understanding of the relationship between TG, TC, and PL metabolism via PLTP.

Genetic fine-mapping reveals single nucleotide polymorphism mutations in the MC1R regulatory region associated with duck melanism.

Birds are among the most colourful terrestrial vertebrates, with various plumage colours and patterns. We conducted a genome-wide association study (GWAS) on an intercross F2 population of Pekin ducks and mallards (n=722) and identified a 1.57-Mb genetic region (Chr11: 20,176,480-21,750,101 bp) related to duck melanism. Fine mapping by linkage disequilibrium (LD) and FST analysis narrowed the final candidate region to a region of 22,500 bp (Chr11: 20,677,500-20,700,000 bp) including three coding genes, TCF25, MC1R and TUBB3. Combined with transcriptome and qRT-PCR analysis, MC1R was identified as the unique genetic locus responsible for black plumage in ducks, and it was significantly more highly expressed in the feather bulbs of black ducks. We also identified 52G > A (Chr11: 20,696,354G > A) and 376G > A (Chr11: 20,696,678G > A) mutations in the MC1R coding region that have been widely studied in ducks. In addition, structural variations (SVs) were screened by nanopore sequencing, and no significant SV was found to be associated with the duck black plumage trait. However, we identified four novel single nucleotide polymorphisms in the MC1R regulator region (Chr11: 20,678,412G > A, Chr11: 20,679,236G > A, Chr11: 20,692,496 A > G and Chr11: 20,692,791 A > G) that had a strong association with the black plumage phenotype of ducks and combined with potential changes in transcription binding affinities. The luciferase reporter gene assay demonstrated that Chr11: 20,678,412G > A and Chr11: 20,679,236G > A led to significant promoter activity changes. Our research emphasizes the importance of MC1R regulatory region mutation in determining the duck black plumage phenotype, and these results expand our understanding of the genetic mechanism underlying duck plumage colour.

Genome-wide association study revealed the genomic regions associated with skin pigmentation in an Ogye x White Leghorn F2 chicken population

Skin color in chickens is an economically important trait that determines the first impression of a consumer toward a broiler and can ultimately affect consumer choice in the market. Therefore, identification of genomic regions associated with skin color is crucial for increasing the sales value of chickens. Although previous studies have attempted to reveal the genetic markers associated with the skin coloration in chickens, most were limited to investigations of candidate genes, such as melanin-related genes, and focused on case/control studies based on a single or small population. In this study, we performed a genome-wide association study (GWAS) on 770 F2 intercrosses produced by an experimental population of 2 chicken breeds, namely Ogye and White Leghorns, with different skin colors. The GWAS demonstrated that the L* value among the 3 skin color traits is highly heritable, and the genomic regions located on 2 chromosomes (20 and Z) were detected to harbor SNPs significantly associated with the skin color trait, accounting for most of the total genetic variance. Particular genomic regions spanning a ∼2.94 Mb region on GGA Z and a ∼3.58 Mb region on GGA 20 were significantly associated with skin color traits, and in these regions, certain candidate genes, including MTAP, FEM1C, GNAS, and EDN3, were found. Our findings could help elucidate the genetic mechanisms underlying chicken skin pigmentation. Furthermore, the candidate genes can be used to provide a valuable breeding strategy for the selection of specific chicken breeds with ideal skin coloration.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population.

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback.

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

Identifying QTLs for Grain Size in a Colossal Grain Rice (Oryza sativa L.) Line, and Analysis of Additive Effects of QTLs.

Grain size is an important component of quality and harvest traits in the field of rice breeding. Although numerous quantitative trait loci (QTLs) of grain size in rice have been reported, the molecular mechanisms of these QTLs remain poorly understood, and further research on QTL observation and candidate gene identification is warranted. In our research, we developed a suite of F2 intercross populations from a cross of 9311 and CG. These primary populations were used to map QTLs conferring grain size, evaluated across three environments, and then subjected to bulked-segregant analysis-seq (BSA-seq). In total, 4, 11, 12 and 14 QTLs for grain length (GL), grain width (GW), 1000-grain weight (TGW), and length/width ratio (LWR), respectively, were detected on the basis of a single-environment analysis. In particular, over 200 splicing-related sites were identified by whole-genome sequencing, including one splicing-site mutation with G>A at the beginning of intron 4 on Os03g0841800 (qGL3.3), producing a smaller open reading frame, without the third and fourth exons. A previous study revealed that the loss-of-function allele caused by this splicing site can negatively regulate rice grain length. Furthermore, qTGW2.1 and qGW2.3 were new QTLs for grain width. We used the near-isogenic lines (NILs) of these GW QTLs to study their genetic effects on individuals and pyramiding, and found that they have additive effects on GW. In summary, these discoveries provide a valuable genetic resource, which will facilitate further study of the genetic polymorphism of new rice varieties in rice breeding.

Strategies of preserving genetic diversity while maximizing genetic response from implementing genomic selection in pulse breeding programs

Key messageGenomic selection maximizes genetic gain by recycling parents to germplasm pool earlier and preserves genetic diversity by restricting the number of fixed alleles and the relationship in pulse breeding programs.Using a stochastic computer simulation, we investigated the benefit of optimization strategies in the context of genomic selection (GS) for pulse breeding programs. We simulated GS for moderately complex to highly complex traits such as disease resistance, grain weight and grain yield in multiple environments with a high level of genotype-by-environment interaction for grain yield. GS led to higher genetic gain per unit of time and higher genetic diversity loss than phenotypic selection by shortening the breeding cycle time. The genetic gain obtained from selecting the segregating parents early in the breeding cycle (at F1 or F2 stages) was substantially higher than selecting at later stages even though prediction accuracy was moderate. Increasing the number of F1 intercross (F1i) families and keeping the total number of progeny of F1i families constant, we observed a decrease in genetic gain and increase in genetic diversity, whereas increasing the number of progeny per F1i family while keeping a constant number of F1i families increased the rate of genetic gain and had higher genetic diversity loss per unit of time. Adding 50 F2 family phenotypes to the training population increased the accuracy of genomic breeding values (GEBVs) and genetic gain per year and decreased the rate of genetic diversity loss. Genetic diversity could be preserved by applying a strategy that restricted both the percentage of alleles fixed and the average relationship of the group of selected parents to preserve long-term genetic improvement in the pulse breeding program.

Triatoma sordida (Hemiptera, Triatominae) from La Paz, Bolivia: an incipient species or an intraspecific chromosomal polymorphism?

BackgroundTriatoma sordida is one of the main Chagas disease vectors in Brazil. In addition to Brazil, this species has already been reported in Bolivia, Argentina, Paraguay, and Uruguay. It is hypothesized that the insects currently identified as T. sordida are a species subcomplex formed by three cytotypes (T. sordida sensu stricto [s.s.], T. sordida La Paz, and T. sordida Argentina). With the recent description of T. rosai from the Argentinean specimens, it became necessary to assess the taxonomic status of T. sordida from La Paz, Bolivia, since it was suggested that it may represent a new species, which has taxonomic, evolutionary, and epidemiological implications. Based on the above, we carried out molecular and experimental crossover studies to assess the specific status of T. sordida La Paz.MethodsTo evaluate the pre- and postzygotic barriers between T. sordida La Paz and T. sordida s.s., experimental crosses and intercrosses between F1 hybrids and between F2 hybrids were conducted. In addition, cytogenetic analyses of the F1 and F2 hybrids were applied with an emphasis on the degree of pairing between the homeologous chromosomes, and morphological analyses of the male gonads were performed to evaluate the presence of gonadal dysgenesis. Lastly, the genetic distance between T. sordida La Paz and T. sordida s.s. was calculated for the CYTB, ND1, and ITS1 genes.ResultsRegardless of the gene used, T. sordida La Paz showed low genetic distance compared to T. sordida s.s. (below 2%). Experimental crosses resulted in offspring for both directions, demonstrating that there are no prezygotic barriers installed between these allopatric populations. Furthermore, postzygotic barriers were not observed either (since the F1 × F1 and F2 × F2 intercrosses resulted in viable offspring). Morphological and cytogenetic analyses of the male gonads of the F1 and F2 offspring demonstrated that the testes were not atrophied and did not show chromosome pairing errors.ConclusionBased on the low genetic distance (which configures intraspecific variation), associated with the absence of prezygotic and postzygotic reproductive barriers, we confirm that T. sordida La Paz represents only a chromosomal polymorphism of T. sordida s.s.Graphical abstract

A ROR2 coding variant is associated with craniofacial variation in domestic pigeons

Vertebrate craniofacial morphogenesis is a highly orchestrated process that is directed by evolutionarily conserved developmental pathways.1,2 Within species, canalized development typically produces modest morphological variation. However, as a result of millennia of artificial selection, the domestic pigeon displays radical craniofacial variation within a single species. One of the most striking cases of pigeon craniofacial variation is the short-beak phenotype, which has been selected in numerous breeds. Classical genetic experiments suggest that pigeon beak length is regulated by a small number of genetic factors, one of which is sex linked (Ku2 locus).3-5 However, the genetic underpinnings of pigeon craniofacial variation remain unknown. Using geometric morphometrics and quantitative trait locus (QTL) mapping on an F2 intercross between a short-beaked Old German Owl (OGO) and a medium-beaked Racing Homer (RH), we identified a single Z chromosome locus that explains a majority of the variation in beak morphology in the F2 population. Complementary comparative genomic analyses revealed that the same locus is strongly differentiated between breeds with short and medium beaks. Within the Ku2 locus, we identified an amino acid substitution in the non-canonical Wnt receptor ROR2 as a putative regulator of pigeon beak length. The non-canonical Wnt pathway serves critical roles in vertebrate neural crest cell migration and craniofacial morphogenesis.6,7 In humans, ROR2 mutations cause Robinow syndrome, a congenital disorder characterized by skeletal abnormalities, including a widened and shortened facial skeleton.8,9 Our results illustrate how the extraordinary craniofacial variation among pigeons can reveal genetic regulators of vertebrate craniofacial diversity.

GWAS and Post-GWAS High-Resolution Mapping Analyses Identify Strong Novel Candidate Genes Influencing the Fatty Acid Composition of the Longissimus dorsi Muscle in Pigs.

Fatty acid (FA) composition is one of the most important parameters for the assessment of meat quality in pigs. The FA composition in pork can also affect human health. Our aim was to identify quantitative trait loci (QTLs) and positional candidate genes affecting the FA profile of the longissimus dorsi muscle in a large F2 intercross between Landrace and Korean native pigs comprising 1105 F2 progeny by genome-wide association studies (GWAS) and post-GWAS high-resolution mapping analyses. We performed GWAS using the PorcineSNP60K BeadChip and a linear mixed model. Four genome-wide significant QTL regions in SSC8, SSC12, SSC14, and SSC16 were detected (p < 2.53 × 10−7). Several co-localizations of QTLs in SSC12 for oleic acid, linoleic acid, arachidonic acid, monounsaturated FAs, polyunsaturated FAs, and the polyunsaturated/saturated FA ratio were observed. To refine the QTL region in SSC12, a linkage and linkage disequilibrium analysis was applied and could narrow down the critical region to a 0.749 Mb region. Of the genes in this region, GAS7, MYH2, and MYH3 were identified as strong novel candidate genes based on further conditional association analyses. These findings provide a novel insight into the genetic basis of FA composition in pork and could contribute to the improvement of pork quality.

Mapping QTLs for Breast Muscle Weight in an F2 Intercross between Native Japanese Nagoya and White Plymouth Rock Chicken Breeds.

Nagoya (NAG), a native Japanese chicken breed, has high quality meat but low meat yield, whereas White Plymouth Rock (WPR), a parental breed of commercial broilers, has rapid growth but high body fat. We previously reported three quantitative trait loci (QTLs) for early postnatal growth in 239 F2 chickens between NAG and WPR breeds. In this study, using the same F2 chickens at 4 weeks of age, we performed genome-wide QTL analysis for breast muscle weight, fat weight and serum and liver levels of biochemical parameters. Two significant QTLs for pectoralis minor and/or major weights were revealed on chromosome 2 between 108 Mb and 127 Mb and chromosome 4 between 10 Mb and 68 Mb. However, no QTL for the other traits was detected. The two QTLs explained 7.0–11.1% of the phenotypic variances, and their alleles derived from WPR increased muscle weights. The chromosome 2 QTL may be a novel locus, whereas the chromosome 4 QTL coincided with a known QTL for meat quality. The findings provide information that is beneficial for genetic improvement of meat yield for the lean NAG breed and, furthermore, provide a better understanding of the genetic basis of chicken muscle development.

Genetic controls of Tas1r3-independent sucrose consumption in mice.

We have previously used crosses between C57BL/6ByJ (B6) and 129P3/J (129) inbred strains to map a quantitative trait locus (QTL) on mouse chromosome (Chr) 4 that affects behavioral and neural responses to sucrose. We have named it the sucrose consumption QTL 2 (Scon2), and shown that it corresponds to the Tas1r3 gene, which encodes a sweet taste receptor subunit TAS1R3. To discover other sucrose consumption QTLs, we have intercrossed B6 inbred and 129.B6-Tas1r3 congenic mice to produce F2 hybrids, in which Scon2 (Tas1r3) does not segregate, and hence does not contribute to phenotypical variation. Chromosome mapping using this F2 intercross identified two main-effect QTLs, Scon3 (Chr9) and Scon10 (Chr14), and an epistatically interacting QTL pair Scon3 (Chr9)-Scon4 (Chr1). Using serial backcrosses, congenic and consomic strains, we conducted high-resolution mapping of Scon3 and Scon4 and analyzed their epistatic interactions. We used mice with different Scon3 or Scon4 genotypes to understand whether these two QTLs influence sucrose intake via gustatory or postoral mechanisms. These studies found no evidence for involvement of the taste mechanisms, but suggested involvement of energy metabolism. Mice with the B6 Scon4 genotype drank less sucrose in two-bottle tests, and also had ahigher respiratory exchange ratio and lower energy expenditure under basal conditions (when they had only chow and water available). Our results provide evidence that Scon3 and Scon4 influence mouse-to-mouse variation in sucrose intake and that both likely act through a common postoral mechanism.

Fitness maps to a large-effect locus in introduced stickleback populations

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin (Eda), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution.

Independent evolution toward larger body size in the distinctive Faroe Island mice.

Most phenotypic traits in nature involve the collective action of many genes. Traits that evolve repeatedly are particularly useful for understanding how selection may act on changing trait values. In mice, large body size has evolved repeatedly on islands and under artificial selection in the laboratory. Identifying the loci and genes involved in this process may shed light on the evolution of complex, polygenic traits. Here, we have mapped the genetic basis of body size variation by making a genetic cross between mice from the Faroe Islands, which are among the largest and most distinctive natural populations of mice in the world, and a laboratory mouse strain selected for small body size, SM/J. Using this F2 intercross of 841 animals, we have identified 111 loci controlling various aspects of body size, weight and growth hormone levels. By comparing against other studies, including the use of a joint meta-analysis, we found that the loci involved in the evolution of large size in the Faroese mice were largely independent from those of a different island population or other laboratory strains. We hypothesize that colonization bottleneck, historical hybridization, or the redundancy between multiple loci have resulted in the Faroese mice achieving an outwardly similar phenotype through a distinct evolutionary path.

A &lt;i&gt;ROR2&lt;/i&gt; Coding Variant is Associated with Craniofacial Variation in Domestic Pigeons

Vertebrate craniofacial morphogenesis is a highly orchestrated process that is directed by evolutionarily conserved developmental pathways. Within species, canalized developmental programs typically produce only modest morphological variation. However, as a result of millennia of artificial selection, the domestic pigeon (Columba livia) displays radical variation in craniofacial morphology within a single species. One of the most striking cases of pigeon craniofacial variation is the short beak phenotype, which has been selected in numerous breeds. Classical genetic experiments suggest that pigeon beak length is regulated by a small number of genetic factors, one of which is sex-linked (Ku2 locus). However, the molecular genetic underpinnings of pigeon craniofacial variation remain unknown. To determine the genetic basis of the short beak phenotype, we used geometric morphometrics and quantitative trait loci (QTL) mapping on an F2 intercross between a short-beaked Old German Owl (OGO) and a medium-beaked Racing Homer (RH). We identified a single locus on the Z-chromosome that explains a majority of the variation in beak morphology in the RH x OGO F2 population. In complementary comparative genomic analyses, we found that the same locus is also strongly differentiated between breeds with short and medium beaks. Within the differentiated Ku2 locus, we identified an amino acid substitution in the non-canonical Wnt receptor ROR2 as a putative regulator of pigeon beak length. The non-canonical Wnt (planar cell polarity) pathway serves critical roles in vertebrate neural crest cell migration and craniofacial morphogenesis. In humans, the ROR2 locus has been identified in several large genome-wide association studies of three-dimensional (3D) normal-range facial variation and homozygous ROR2 mutations cause autosomal recessive Robinow syndrome, a rare congenital disorder characterized by skeletal abnormalities, including a widened and shortened facial skeleton. Our results illustrate how the extraordinary craniofacial variation among pigeons can reveal genetic regulators of vertebrate craniofacial diversity.

Epigenetics of Animal Personality: DNA Methylation Cannot Explain the Heritability of Exploratory Behavior in a Songbird.

SynopsisThe search for the hereditary mechanisms underlying quantitative traits traditionally focused on the identification of underlying genomic polymorphisms such as single-nucleotide polymorphisms. It has now become clear that epigenetic mechanisms, such as DNA methylation, can consistently alter gene expression over multiple generations. It is unclear, however, if and how DNA methylation can stably be transferred from one generation to the next and can thereby be a component of the heritable variation of a trait. In this study, we explore whether DNA methylation responds to phenotypic selection using whole-genome and genome-wide bisulfite approaches. We assessed differential erythrocyte DNA methylation patterns between extreme personality types in the Great Tit (Parus major). For this, we used individuals from a four-generation artificial bi-directional selection experiment and siblings from eight F2 inter-cross families. We find no differentially methylated sites when comparing the selected personality lines, providing no evidence for the so-called epialleles associated with exploratory behavior. Using a pair-wise sibling design in the F2 intercrosses, we show that the genome-wide DNA methylation profiles of individuals are mainly explained by family structure, indicating that the majority of variation in DNA methylation in CpG sites between individuals can be explained by genetic differences. Although we found some candidates explaining behavioral differences between F2 siblings, we could not confirm this with a whole-genome approach, thereby confirming the absence of epialleles in these F2 intercrosses. We conclude that while epigenetic variation may underlie phenotypic variation in behavioral traits, we were not able to find evidence that DNA methylation can explain heritable variation in personality traits in Great Tits.

Mapping the Effects of Genetic Variation on Chromatin State and Gene Expression Reveals Loci That Control Ground State Pluripotency.

Mouse embryonic stem cells (mESCs) cultured in the presence of LIF occupy a ground state with highly active pluripotency-associated transcriptional and epigenetic circuitry. However, ground state pluripotency in some inbred strain backgrounds is unstable in the absence of ERK1/2 and GSK3 inhibition. Using an unbiased genetic approach, we dissect the basis of this divergent response to extracellular cues by profiling gene expression and chromatin accessibility in 170 genetically heterogeneous mESCs. We map thousands of loci affecting chromatin accessibility and/or transcript abundance, including 10 QTL hotspots where genetic variation at a single locus coordinates the regulation of genes throughout the genome. For one hotspot, we identify a single enhancer variant ∼10 kb upstream of Lifr associated with chromatin accessibility and mediating a cascade of molecular events affecting pluripotency. We validate causation through reciprocal allele swaps, demonstrating the functional consequences of noncoding variation in gene regulatory networks that stabilize pluripotent states invitro.