Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs

The recent availability of a genome‐wide SNP array for the goat genome dramatically increases the power to investigate aspects of genetic diversity and to conduct genome‐wide association studies in this important domestic species. We collected and analysed genotypes from 52 088 SNPs in Boer, Cashmere and Rangeland goats that had both polled and horned individuals. Principal components analysis revealed a clear genetic division between animals for each population, and model‐based clustering successfully detected evidence of admixture that matched aspects of their recorded history. For example, shared co‐ancestry was detected, suggesting Boer goats have been introgressed into the Rangeland population. Further, allele frequency data successfully tracked the altered genetic profile that has taken place after 40 years of breeding Australian Cashmere goats using the Rangeland animals as the founding population. Genome‐wide association mapping of the POLL locus revealed a strong signal on goat chromosome 1. The 769‐kb critical interval contained the polled intersex syndrome locus, confirming the genetic basis in non‐European animals is the same as identified previously in Saanen goats. Interestingly, analysis of the haplotypes carried by a small set of sex‐reversed animals, known to be associated with polledness, revealed some animals carried the wild‐type chromosome associated with the presence of horns. This suggests a more complex basis for the relationship between polledness and the intersex condition than initially thought while validating the application of the goat SNP50 BeadChip for fine‐mapping traits in goat.     

Genetic mapping of the autosomal region involved in XX sex-reversal and horn development in goats

Contrary to other genetic disorders, the genetic study of sex determination anomalies in humans stumbles over the difficulty in observing large pedigrees. In goats, abnormalities in sex determination are intimately linked to a dominant Mendelian gene coding for the ``polled' (hornless) character, which could render this species an interesting animal model for the rare human cases of SRY-negative XX males. In this report, we describe genetic linkage between the polled/intersex synchome (PIS) and four microsatellite markers of the distal region of goat Chromosome 1 (CHI1), quite distinct from the bovine ``polled' region. According to comparative mapping data, no sex-determining gene has been described so far in homologous regions in the human. This genetic localization constitutes a first step towards identifying a new autosomal sex-determining gene in mammals. Received: 7 July 1995 / Accepted: 6 September 1995     

Characterization of a Y‐specific duplication/insertion of the anti‐Mullerian hormone type II receptor gene based on a chromosome‐scale genome assembly of yellow perch,Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long‐reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi‐C, we generated a high‐continuity chromosome‐scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome‐size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male‐specific duplicate of the anti‐Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex‐specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by⁺) from XX genetic females (amhr2by^?). Our high‐quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex‐determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.     

A complex structural variant at the KIT locus in cattle with the Pinzgauer spotting pattern

A specific white spotting phenotype, termed finching or line‐backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux‐Zillertaler cattle, two Austrian breeds. The so‐called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome‐wide association study using 27 white spotted and 16 solid‐coloured Tux‐Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole‐genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4‐kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310‐kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KIT^PINZ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux‐Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KIT^PINZ variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT^PINZ variant confirms admixture and the reported historical relationship of these short‐headed breeds with Austrian Tux‐Zillertaler and suggests a mutation event, occurring before breed formation.     

Identification of small and large genomic candidate variants in bovine pulmonary hypoplasia and anasarca syndrome

The pulmonary hypoplasia and anasarca syndrome (PHA) is a congenital lethal disorder, which until now has been reported in cattle and sheep. PHA is characterized by extensive subcutaneous fetal edema combined with hypoplasia or aplasia of the lungs and dysplasia of the lymphatic system. PHA is assumed to be of genetic etiology. This study presents the occurrence of PHA in two different cattle breeds and their genetic causation. Two PHA cases from one sire were observed in Slovenian Cika cattle. Under the assumption of monogenic inheritance, genome-wide homozygosity mapping scaled down the critical regions to 3% of the bovine genome including a 43.6 Mb-sized segment on chromosome 6. Whole-genome sequencing of one case, variant filtering against controls and genotyping of a larger cohort of Cika cattle led to the detection of a likely pathogenic protein-changing variant perfectly associated with the disease: a missense variant on chromosome 6 in ADAMTS3 (NM_001192797.1: c.1222C>T), which affects an evolutionary conserved residue (NP_001179726.1: p.(His408Tyr)). A single PHA case was found in Danish Holstein cattle and was whole-genome sequenced along with its parents. However, as there was no plausible private protein-changing variant, mining for structural variation revealed a likely pathogenic trisomy of the entire chromosome 20. The identified ADAMTS3 associated missense variant and the trisomy 20 are two different genetic causes, which shows a compelling genetic heterogeneity for bovine PHA.     

Sry-negative XX true hermaphroditism in a roe deer

A two-year-old roe deer was brought down in the course of a hunt in the north of Spain (Asturias). On physical examination the individual presented well-developed bared antlers, but surprisingly a female external genitalia. Several anatomical, histological and genetic analyses were performed in order to explain the observed phenotype. Necropsy evidenced ovary-like structures with follicles on the surface; histological analyses of testes evidenced positive immunolabel against testosterone in Leydig cells; genetic analyses showed that the sex of the individual was consistent with a female individual. PCR analysis failed to detect SRY sequences; no PIS deletion, which is responsible for XX sex-reversal in goats, was detected. On the basis of its presumptive normal female sexual karyotype (XX) and the presence of two functional abdominal bilateral testes and ovaries, the roe deer was finally diagnosed as possessing an XX hermaphroditism syndrome. However, as in many other cases, the specific reason for the occurrence of this case of hermaphroditism could not be determined.     

Genome‐wide association study confirms that the chromosome Z harbours a region responsible for rumplessness in Hongshan chickens

Rumplessness in Hongshan chickens has been reported as a novel sex‐linked characteristic. Re‐sequencing data suggest that a pseudogene on the Z chromosome, LOC431648, is affiliated with this phenotype. In this study, we chose 23 rumpless and 25 normal Hongshan chickens to localize the potential variation by means of a genome‐wide association study using a high density microarray. A region on the Z chromosome was found to be closely associated with rumplessness in Hongshan chickens. The region, located in gene LINGO2, was approximately 3 Mb away from pseudogene LOC431648. The function of this gene has not yet been studied in birds. Differential expression of the candidate genes in the tail feather follicles was not detected by q‐PCR, which suggests that the rumplessness trait could be attributed to other genetic mechanisms.     

Whole genome sequencing identifies a missense polymorphism in PADI6 associated with testicular/ovotesticular XX disorder of sex development in dogs

Disorders of sex development (DSDs) are congenital malformations defined as discrepancies between sex chromosomes and phenotypical sex. Testicular or ovotesticular XX DSDs are frequently observed in female dogs, while monogenic XY DSDs are less frequent. Here, we applied whole genome sequencing (WGS) to search for causative mutations in XX DSD females in French Bulldogs (FB) and American Staffordshire Terries (AST) and in XY DSD Yorkshire Terries (YT). The WGS results were validated by Sanger sequencing and ddPCR. It was shown that a missense SNP of the PADI6 gene, is significantly associated with the XX DSD (SRY-negative) phenotype in AST (P = 0.0051) and FB (P = 0.0306). On the contrary, we did not find any associated variant with XY DSD in YTs. Our study suggests that the genetic background of the XX DSD may be more complex and breed-specific.     

De novo assembly of a chromosome‐level reference genome of red‐spotted grouper (Epinephelus akaara) using nanopore sequencing and Hi‐C

The red‐spotted grouper Epinephelus akaara (E. akaara) is one of the most economically important marine fish in China, Japan and South‐East Asia and is a threatened species. The species is also considered a good model for studies of sex inversion, development, genetic diversity and immunity. Despite its importance, molecular resources for E. akaara remain limited and no reference genome has been published to date. In this study, we constructed a chromosome‐level reference genome of E. akaara by taking advantage of long‐read single‐molecule sequencing and de novo assembly by Oxford Nanopore Technology (ONT) and Hi‐C. A red‐spotted grouper genome of 1.135 Gb was assembled from a total of 106.29 Gb polished Nanopore sequence (GridION, ONT), equivalent to 96‐fold genome coverage. The assembled genome represents 96.8% completeness (BUSCO) with a contig N50 length of 5.25 Mb and a longest contig of 25.75 Mb. The contigs were clustered and ordered onto 24 pseudochromosomes covering approximately 95.55% of the genome assembly with Hi‐C data, with a scaffold N50 length of 46.03 Mb. The genome contained 43.02% repeat sequences and 5,480 noncoding RNAs. Furthermore, combined with several RNA‐seq data sets, 23,808 (99.5%) genes were functionally annotated from a total of 23,923 predicted protein‐coding sequences. The high‐quality chromosome‐level reference genome of E. akaara was assembled for the first time and will be a valuable resource for molecular breeding and functional genomics studies of red‐spotted grouper in the future.     

Genome‐wide association reveals the locus responsible for four‐horned ruminant

Phenotypic variability in horn characteristics, such as their size, number and shape, offers the opportunity to elucidate the molecular basis of horn development. The objective of this study was to map the genetic determinant controlling the production of four horns in two breeds, Jacob sheep and Navajo‐Churro, and examine whether an eyelid abnormality occurring in the same populations is related. Genome‐wide association mapping was performed using 125 animals from the two breeds that contain two‐ and four‐horned individuals. A case–control design analysis of 570 712 SNPs genotyped with the ovine HD SNP Beadchip revealed a strong association signal on sheep chromosome 2. The 10 most strongly associated SNPs were all located in a region spanning Mb positions 131.9–132.6, indicating the genetic architecture underpinning the production of four horns is likely to involve a single gene. The closest genes to the most strongly associated marker (OAR2_132568092) were MTX2 and the HOXD cluster, located approximately 93 Kb and 251 Kb upstream respectively. The occurrence of an eyelid malformation across both breeds was restricted to polled animals and those carrying more than two horns. This suggests the eyelid abnormality may be associated with departures from the normal developmental production of two‐horned animals and that the two conditions are developmentally linked. This study demonstrated the presence of separate loci responsible for the polled and four‐horned phenotypes in sheep and advanced our understanding of the complexity that underpins horn morphology in ruminants.     

Characterization of Pisrt1/Foxl2 in Ellobius lutescens and exclusion as sex-determining genes

The rodent Ellobius lutescens is an exceptional mammal which determines male sex constitutively without the SRY gene and, therefore, may serve as an animal model for human 46,XX female-to-male sex reversal. It was suggested that other factors of the network of sex-determining genes determine maleness in these animals. However, some sex-determining genes like SOX9 and SF1 have already been excluded by segregation analysis as primary sex-determining factors in E. lutescens. In this work, we have cloned and characterized two genes of the PIS (polled intersex syndrome) gene interval, which were reported as candidates in female-to-male sex reversal in hornless goats recently. The genes Foxl2 and Pisrt1 from that interval were identified in E. lutescens DNA and mapped to Chromosome 8. We have excluded linkage of Foxl2 and Pisrt1 loci with the sex of the animals. Hence, the involvement of this gene region in sex determination may be specific for goats and is not a general mechanism of XX sex reversal or XX male sex determination.The nucleotide sequence data reported in this article have been submitted to GenBank and have been assigned the accession number AY623815.     

Two loss-of-function alleles of the glutathione S-transferase (GST) gene cause anthocyanin deficiency in flower and fruit skin of peach (Prunus persica)

Flower and fruit colors are important agronomic traits. To date, there is no forward genetic evidence that the glutathione S-transferase (GST) gene is responsible for the white flower color in peach (Prunus persica). In this study, genetic analysis indicated that the white-flower trait is monogenetic, is recessive to the non-white allele, and shows pleiotropic effects with non-white-flowered types. The genetic locus underpinning this trait was mapped onto chromosome 3 between 0.421951 and 3.227115 Mb by using bulked segregant analysis in conjunction with whole-genome sequencing, and was further mapped between 0 and 1.178149 Mb by using the backcross 1 (BC₁) population. Finally, the locus was fine-mapped within 535.974- and 552.027-kb intervals by using 151 F₂ individuals and 75 individuals from a BC₁ self-pollinated (BC₁S₁) population, respectively. Pp3G013600, encoding a GST that is known to transport anthocyanin, was identified within the mapping interval. The analysis of genome sequence data showed Pp3G013600 in white flowers has a 2-bp insertion or a 5-bp deletion in the third exon. These variants likely render the GST non-functional because of early stop codons that reduce the protein length from 215 amino acids to 167 and 175 amino acids, respectively. Genetic markers based on these variants validated a complete correlation between the GST loss-of-function alleles and white flower in 128 peach accessions. This correlation was further confirmed by silencing of Pp3G013600 using virus-induced gene silencing technology, which reduced anthocyanin accumulation in peach fruit. The new knowledge from this study is useful for designing peach breeding programs to generate cultivars with white flower and fruit skin.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.     

奶山羊PIS遗传缺陷基因的检测及应用

文章旨在建立一种奶山羊无角间性综合征(Polled intersex syndrome,PIS)遗传缺陷基因检测方法。根据PIS基因序列(AF404302)分别设计PIS、PIS+、NEI 3对扩增引物,利用PCR技术鉴定奶山羊PIS遗传缺陷基因型。基因型为PIS PIS+与PIS PIS+的表型正常个体分别扩增出(141,300 bp)和(141,449,300 bp)的片段组合,隐性纯合间性山羊(PIS+PIS+)扩增出(449,300 bp)的片段组合。利用该方法检测一个224个体的奶山羊群体,结果显示:PIS PIS、PIS PIS+和PIS+PIS+个体分别为70、150和4个。该群体中PIS PIS+基因型频率高达66.9%,PIS+基因频率为35.3%,其子代群体出现间性山羊的可能性将超过12%。文章所开发的奶山羊PIS遗传缺陷基因检测方法可直接准确判别种公羊的基因型,从而避免缺陷基因携带者种公羊的使用。该方法易操作且准确性高,对奶山羊的标记辅助选择及奶山羊产业的健康发展具有重要意义。     

Whole‐genome sequencing reveals a large deletion in the MITF gene in horses with white spotted coat colour and increased risk of deafness

White spotting phenotypes in horses are highly valued in some breeds. They are quite variable and may range from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for white spotting phenotypes in horses. For the present study, we investigated an American Paint Horse family segregating a phenotype involving white spotting and blue eyes. Six of eight horses with the white‐spotting phenotype were deaf. We obtained whole‐genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~63‐kb deletion spanning exons 6–9 of the MITF gene (chr16:21 503 211–21 566 617). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. PCR‐based genotyping revealed that all eight available affected horses from the family carried the deletion. The finding of an MITF variant fits well with the syndromic phenotype involving both depigmentation and an increased risk for deafness and corresponds to human Waardenburg syndrome type 2A. Our findings will enable more precise genetic testing for depigmentation phenotypes in horses.     

Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.     

Deep sequencing of a candidate region harboring the SOX9 gene for the canine XX disorder of sex development

A disorder of sex development (DSD) in dogs with female sex chromosomes (78, XX), a lack of the SRY gene and the presence of testes or ovotestes is commonly diagnosed in numerous breeds. The molecular background of DSD is not fully recognized but has been linked to the copy number variation in the region harboring the SOX9 gene. We applied a genome‐wide association study and targeted next‐generation sequencing techniques to compare DSD and normal female dogs. The genome‐wide association study did not indicate a significant chromosome region. Targeted next‐generation sequencing of a 1.5‐Mb region on canine chromosome 9 harboring the SOX9 gene revealed two putatively DSD‐associated copy number variations 355 kb upstream and 691 kb downstream of SOX9, four blocks of low polymorphism and two blocks of an elevated heterozygosity. An initial next‐generation sequencing analysis showed an association with two SNPs, but validation in larger cohorts did not confirm this result. We identified a large homologous fragment (over 243.8 kb), named hfMAGI2, located upstream of SOX9, that overlaps a known copy number variation region. It shows a high sequence similarity with the 5′ flanking region of the MAGI2 gene located on canine chromosome 18 that encodes a protein involved in ovary formation during early embryonic development. Our study showed that the identified copy number variation region located upstream of the SOX9 gene contains potential regulatory sequences (long non‐coding RNA and hfMAGI2) and led to the assumption that a multiplication of this element may alter expression of the SOX9 gene, triggering the DSD phenotype.     

<Emphasis Type="Italic">R-spondin1</Emphasis> and <Emphasis Type="Italic">FOXL2</Emphasis>act into two distinct cellular types during goat ovarian differentiation

Background  

Up to now, two loci have been involved in XX sex-reversal in mammals following loss-of-function mutations, PIS (Polled Intersex Syndrome) in goats and R-spondin1 (RSPO1) in humans. Here, we analyze the possible interaction between these two factors during goat gonad development. Furthermore, since functional redundancy between different R-spondins may influence gonad development, we also studied the expression patterns of RSPO2, 3 and 4.