Doubled haploid production in fruit crops

The interest of fruit breeders in haploids and doubled haploids (DH), lies in the possibility of shortening the time needed to produce homozygous lines compared to conventional breeding. Haplo-diploidization through gametic embryogenesis allows single-step development of complete homozygous lines from heterozygous parents. In a conventional breeding programme, a pure line is developed after several generations of selfing. With fruit crops, characterized by a long reproductive cycle, a high degree of heterozygosity, large size, and, sometimes, self-incompatibility, there is no way to obtain haploidization through conventional methods. This paper reviews the current status of research on doubled haploid production in the main fruit crops: Citrus, Malus domestica, Pyrus communis, Pyrus pyrifolia, Prunus persica, Prunus avium, Prunus domestica, Prunus armeniaca, Vitis vinifera, Actinidia deliciosa, Olea europaea, Morus alba, Actinidia deliziosa, [Musa balbisiana (BB)], Carica papaya, Annona squamosa, Feijoa sellowiana, Opuntia ficus-indica, Eriobotrya japonica.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

The Wright-Fisher model with temporally varying selection and population size

The Wright-Fisher model is considered in the case where the population size is random and the magnitude of the selective advantage of one of the alleles varies with time. The central question addressed is the possibility of ultimate genetic polymorphism. Partial results are obtained in the general case and complete results in the case where the population size and selective advantage are not density dependent. Bounds on the fixation probability are obtained when the selective advantage is constant.     

Estimations of linkage disequilibrium,effective population size and ROH‐based inbreeding coefficients in Spanish Churra sheep using imputed high‐density SNP genotypes

In this study, the availability of the Ovine HD SNP BeadChip (HD‐chip) and the development of an imputation strategy provided an opportunity to further investigate the extent of linkage disequilibrium (LD) at short distances in the genome of the Spanish Churra dairy sheep breed. A population of 1686 animals, including 16 rams and their half‐sib daughters, previously genotyped for the 50K‐chip, was imputed to the HD‐chip density based on a reference population of 335 individuals. After assessing the imputation accuracy for beagle v4.0 (0.922) and fimpute v2.2 (0.921) using a cross‐validation approach, the imputed HD‐chip genotypes obtained with beagle were used to update the estimates of LD and effective population size for the studied population. The imputed genotypes were also used to assess the degree of homozygosity by calculating runs of homozygosity and to obtain genomic‐based inbreeding coefficients. The updated LD estimations provided evidence that the extent of LD in Churra sheep is even shorter than that reported based on the 50K‐chip and is one of the shortest extents compared with other sheep breeds. Through different comparisons we have also assessed the impact of imputation on LD and effective population size estimates. The inbreeding coefficient, considering the total length of the run of homozygosity, showed an average estimate (0.0404) lower than the critical level. Overall, the improved accuracy of the updated LD estimates suggests that the HD‐chip, combined with an imputation strategy, offers a powerful tool that will increase the opportunities to identify genuine marker‐phenotype associations and to successfully implement genomic selection in Churra sheep.     

A case of Usher syndrome type IIA caused by a rare USH2A homozygous frameshift variant with maternal uniparental disomy (UPD) in a Chinese family

Usher syndrome encompasses a group of genetically and clinically heterogeneous autosomal recessive disorders with hearing deficiencies and retinitis pigmentosa. The mechanisms underlying the Usher syndrome are highly variable. In the present study, a Chinese family with Usher syndrome was recruited. Whole exome sequencing (WES), Sanger sequencing, homozygosity mapping, short tandem repeat (STR) analysis and segregation analysis were performed. Functional domains of the pathogenic variant for USH2A were analysed. We identified a homozygous frameshift variant c.99_100insT (p.Arg34Serfs*41) in the USH2A gene in the proband that showed discordant segregation in the father. Further homozygosity mapping and STR analysis identified an unusual homozygous variant of proband that originated from maternal uniparental disomy (UPD). The p.Arg34Serfs*41 variant produced a predicted truncated protein that removes all functional domains of USH2A. The variant was not included in the 1000 Human Genomes Project database, ExAC database, HGMD or gnomAD database, but was included in the ClinVar databases as pathogenic. Although USH2A is an autosomal recessive disease, the effects of UPD should be informed in genetic counselling since the recurrence risk of an affected child is greatly reduced when the disease is due to the UPD mechanism. To test potential patients, WES, combined with STR analysis and homozygosity mapping, provides an accurate and useful strategy for genetic diagnosis. In summary, our discoveries can help further the understanding of the molecular pathogenesis of Usher syndrome type IIA to advance the prevention, diagnosis and therapy for this disorder.     

Locus‐specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post‐admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3–47.3 Mb) and 18 (18.7–25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre‐ and post‐admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6–24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (F_st), an alternative indicator of pre‐admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and F_st was observed. We found two regions with significant signatures of post‐admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

Genome‐wide assessment of diversity and differentiation between original and modern Brown cattle populations

Identifying genomic regions involved in the differences between breeds can provide information on genes that are under the influence of both artificial and natural selection. The aim of this study was to assess the genetic diversity and differentiation among four different Brown cattle populations (two original vs. two modern populations) and to characterize the distribution of runs of homozygosity (ROH) islands using the Illumina Bovine SNP50 BeadChip genotyping data. After quality control, 34 735 SNPs and 106 animals were retained for the analyses. Larger heterogeneity was highlighted for the original populations. Patterns of genetic differentiation, multidimensional scaling, and the neighboring joining tree distinguished the modern from the original populations. The F_ST‐outlier identified several genes putatively involved in the genetic differentiation between the two groups, such as stature and growth, behavior, and adaptability to local environments. The ROH islands within both the original and the modern populations overlapped with QTL associated with relevant traits. In modern Brown (Brown Swiss and Italian Brown), ROH islands harbored candidate genes associated with milk production traits, in evident agreement with the artificial selection conducted to improve this trait in these populations. In original Brown (Original Braunvieh and Braunvieh), we identified candidate genes related with fat deposition, confirming that breeding strategies for the original Brown populations aimed to produce dual‐purpose animals. Our study highlighted the presence of several genomic regions that vary between Brown populations, in line with their different breeding histories.     

A 63-bp insertion in exon 2 of the porcine KIF21A gene is associated with arthrogryposis multiplex congenita

A recessive form of arthrogryposis multiplex congenita (AMC) was detected 20 years ago in the Swiss Large White (SLW) pig population. A diagnostic marker test enabled the identification of carrier animals, but the underlying causal mutation remains unknown. To identify the mutation underlying AMC, we collected SNP chip genotyping data for 11 affected piglets and 23 healthy pigs. Association testing using 47 829 SNPs confirmed that AMC maps to SSC5 (P = 9.4 × 10⁻¹³). Subsequent autozygosity mapping revealed a common 6.06 Mb region (from 66 757 970 to 72 815 151 bp) of extended homozygosity in 11 piglets affected by AMC. Using WGS data, we detected a 63-bp insertion compatible with the recessive inheritance of AMC in the second exon of KIF21A gene encoding Kinesin Family Member 21A. The 63-bp insertion is predicted to introduce a premature stop codon in KIF21A gene (p.Val41_Phe42insTer) that truncates 1614 amino acids (~97%) from the protein. We found that this deleterious allele still segregates at a frequency of 0.1% in the SLW pig population. Carrier animals can now be detected unambiguously and excluded from breeding.     

Brachygnathia,cardiomegaly and renal hypoplasia syndrome (BCRHS) in Merino sheep maps to a 1.1‐megabase region on ovine chromosome OAR2

A genome scan was conducted to map the autosomal recessive lethal disorder brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Poll Merino sheep. The scan involved 10 affected and 27 unaffected animals from a single Poll Merino/Merino sheep flock, which were genotyped with the Illumina Ovine SNP50 BeadChip. Association and homozygosity mapping analyses located the disorder in a region comprising 20 consecutive SNPs spanning 1.1 Mb towards the distal end of chromosome OAR2. All affected animals and none of the unaffected animals were homozygous for the associated haplotype in this region. These results provide the basis for identifying the causative mutation(s) and should enable the development of a DNA test to identify carriers in the Poll Merino sheep population. Understanding the molecular control of BCRHS may provide insight into the fundamental genetic control and regulation of the affected organ systems.