Phenotypic Mutation 'corpulent' (pdf version)
Allelecorpulent
Mutation Type critical splice donor site
Chromosome9
Coordinate22,486,492 bp (GRCm39)
Base Change G ⇒ T (forward strand)
Gene Bbs9
Gene Name Bardet-Biedl syndrome 9
Synonym(s) E130103I17Rik, EST 3159894
Chromosomal Location 22,387,011-22,799,576 bp (+) (GRCm39)
MGI Phenotype FUNCTION: [Summary is not available for the mouse gene. This summary is for the human ortholog.] This gene is downregulated by parathyroid hormone in osteoblastic cells, and therefore is thought to be involved in parathyroid hormone action in bones. The exact function of this gene has not yet been determined. Alternatively spliced transcript variants encoding different isoforms have been identified. [provided by RefSeq, Jan 2017]
Accession Number

NCBI RefSeq: NM_178415, NM_181316; MGI:2442833

MappedYes 
Amino Acid Change
Institutional SourceBeutler Lab
Gene Model predicted gene model for protein(s): [ENSMUSP00000043042 ] [ENSMUSP00000120927 ] [ENSMUSP00000122058 ] [ENSMUSP00000116629 ]   † probably from a misspliced transcript
AlphaFold Q811G0
SMART Domains Protein: ENSMUSP00000043042
Gene: ENSMUSG00000035919

DomainStartEndE-ValueType
Pfam:PHTB1_N 1 421 8e-168 PFAM
Pfam:PHTB1_C 439 814 8.3e-163 PFAM
Predicted Effect probably null
SMART Domains Protein: ENSMUSP00000120927
Gene: ENSMUSG00000035919

DomainStartEndE-ValueType
Pfam:PHTB1_N 1 417 1.1e-166 PFAM
Pfam:PHTB1_C 440 818 7e-158 PFAM
Predicted Effect probably null
SMART Domains Protein: ENSMUSP00000120927
Gene: ENSMUSG00000035919

DomainStartEndE-ValueType
Pfam:PHTB1_N 1 417 1.1e-166 PFAM
Pfam:PHTB1_C 440 818 7e-158 PFAM
Predicted Effect probably null
SMART Domains Protein: ENSMUSP00000122058
Gene: ENSMUSG00000035919

DomainStartEndE-ValueType
Pfam:PHTB1_N 1 421 8e-168 PFAM
Pfam:PHTB1_C 439 814 8.3e-163 PFAM
Predicted Effect probably null
SMART Domains Protein: ENSMUSP00000116629
Gene: ENSMUSG00000035919

DomainStartEndE-ValueType
Pfam:PHTB1_N 1 421 8e-168 PFAM
Pfam:PHTB1_C 439 814 8.3e-163 PFAM
Predicted Effect probably null
Meta Mutation Damage Score 0.9496 question?
Is this an essential gene? Probably essential (E-score: 0.870) question?
Phenotypic Category Autosomal Recessive
Candidate Explorer Status loading ...
Single pedigree
Linkage Analysis Data
Penetrance  
Alleles Listed at MGI

All Mutations and Alleles(18) : Gene trapped(16) Targeted(2)

Lab Alleles
AlleleSourceChrCoordTypePredicted EffectPPH Score
IGL01349:Bbs9 APN 9 22798979 missense probably benign 0.00
IGL01586:Bbs9 APN 9 22557293 missense possibly damaging 0.46
IGL01646:Bbs9 APN 9 22582221 nonsense probably null
IGL01654:Bbs9 APN 9 22402238 critical splice donor site probably null
IGL02172:Bbs9 APN 9 22490772 missense possibly damaging 0.65
IGL02212:Bbs9 APN 9 22723808 missense probably benign 0.02
IGL02444:Bbs9 APN 9 22555083 missense probably damaging 0.96
IGL02829:Bbs9 APN 9 22490780 missense probably damaging 0.98
IGL03385:Bbs9 APN 9 22555044 missense probably benign 0.19
Crapulence UTSW 9 22479060 missense probably damaging 1.00
R0038:Bbs9 UTSW 9 22415390 missense probably benign 0.30
R0243:Bbs9 UTSW 9 22425297 missense probably damaging 1.00
R0595:Bbs9 UTSW 9 22408111 missense probably benign
R0688:Bbs9 UTSW 9 22479015 missense probably damaging 0.98
R0726:Bbs9 UTSW 9 22705119 missense probably damaging 0.99
R0749:Bbs9 UTSW 9 22486497 splice site probably null
R0783:Bbs9 UTSW 9 22479010 missense possibly damaging 0.69
R1148:Bbs9 UTSW 9 22486396 splice site probably benign
R1532:Bbs9 UTSW 9 22798945 missense probably benign 0.00
R1783:Bbs9 UTSW 9 22570415 missense possibly damaging 0.85
R2285:Bbs9 UTSW 9 22590230 missense probably damaging 1.00
R2402:Bbs9 UTSW 9 22557359 missense probably benign 0.23
R2655:Bbs9 UTSW 9 22415348 missense probably damaging 1.00
R3428:Bbs9 UTSW 9 22479183 splice site probably benign
R3798:Bbs9 UTSW 9 22550065 missense probably damaging 1.00
R3806:Bbs9 UTSW 9 22798926 missense probably damaging 0.98
R4660:Bbs9 UTSW 9 22490063 missense probably benign 0.16
R4873:Bbs9 UTSW 9 22490011 missense probably benign 0.06
R4875:Bbs9 UTSW 9 22490011 missense probably benign 0.06
R5291:Bbs9 UTSW 9 22540293 missense probably damaging 1.00
R5364:Bbs9 UTSW 9 22486492 critical splice donor site probably null
R5502:Bbs9 UTSW 9 22415370 missense probably damaging 1.00
R5646:Bbs9 UTSW 9 22490011 missense probably benign 0.06
R5932:Bbs9 UTSW 9 22723627 missense probably damaging 1.00
R6222:Bbs9 UTSW 9 22479147 missense possibly damaging 0.88
R6451:Bbs9 UTSW 9 22479060 missense probably damaging 1.00
R6547:Bbs9 UTSW 9 22425365 missense probably benign 0.01
R6726:Bbs9 UTSW 9 22557260 missense probably benign 0.00
R6745:Bbs9 UTSW 9 22582132 missense probably benign 0.00
R6908:Bbs9 UTSW 9 22479019 missense probably damaging 0.96
R6919:Bbs9 UTSW 9 22723840 critical splice donor site probably null
R7102:Bbs9 UTSW 9 22490849 missense probably damaging 1.00
R7536:Bbs9 UTSW 9 22582096 missense probably damaging 1.00
R7712:Bbs9 UTSW 9 22582109 missense probably benign 0.34
R8177:Bbs9 UTSW 9 22425359 missense probably benign 0.08
R8190:Bbs9 UTSW 9 22590284 missense probably damaging 1.00
R8241:Bbs9 UTSW 9 22590214 missense probably benign 0.00
R8440:Bbs9 UTSW 9 22479046 missense probably damaging 0.99
R8733:Bbs9 UTSW 9 22582128 missense probably benign 0.03
R8737:Bbs9 UTSW 9 22590244 missense probably benign 0.17
R8854:Bbs9 UTSW 9 22490060 missense probably damaging 1.00
R8885:Bbs9 UTSW 9 22590234 missense possibly damaging 0.67
R9135:Bbs9 UTSW 9 22490005 missense probably damaging 1.00
X0027:Bbs9 UTSW 9 22566626 missense probably damaging 1.00
Mode of Inheritance Autosomal Recessive
Local Stock
Repository
Last Updated 2019-09-04 9:39 PM by Anne Murray
Record Created 2017-09-05 2:23 PM
Record Posted 2018-08-01
Phenotypic Description

Figure 1. Corpulent mice exhibited increased body weights compared to wild-type littermates. Scaled weights are shown. Abbreviations: REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.

Figure 2. Corpulent mice exhibited increased fat mass compared to wild-type littermates. Normalized data are shown. Abbreviations: REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 3. Corpulent mice exhibited reduced lean fat ratios compared to wild-type littermates. Normalized data are shown. Abbreviations: REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 4. Corpulent mice exhibited reduced pelivis lengths compared to wild-type littermates. Normalized data are shown. Abbreviations: REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.
Figure 5. Corpulent mice exhibited reduced tibia lengths compared to wild-type littermates. Normalized data are shown. Abbreviations: REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.

The corpulent phenotype was identified among N-ethyl-N-nitrosourea (ENU)-mutagenized G3 mice of the pedigree R5364, some of which showed increased body weights compared to wild-type littermates (Figure 1). Some mice also showed increased fat mass (Figure 2), reduced lean fat ratios (Figure 3), reduced pelvis lengths (Figure 4), and reduced tibia lengths (Figure 5).

Nature of Mutation

Figure 6. Linkage mapping of the increased fat mass phenotype using a recessive model of inheritance. Manhattan plot shows -log10 P values (Y-axis) plotted against the chromosome positions of 98 mutations (X-axis) identified in the G1 male of pedigree R5364. Normalized data are shown for single locus linkage analysis without consideration of G2 dam identity. Horizontal pink and red lines represent thresholds of P = 0.05, and the threshold for P = 0.05 after applying Bonferroni correction, respectively.

Whole exome HiSeq sequencing of the G1 grandsire identified 98 mutations. All of the above phenotypes were  linked to a mutation in Bbs9:  a G to T transversion at base pair 22,575,196 (v38) on chromosome 9, or base pair 99,660 in the GenBank genomic region NC_000075 within the splice donor site of intron 8. The strongest association was found with a recessive model of inheritance to the fat mass phenotype, wherein two variant homozygotes departed phenotypically from 20 homozygous reference mice and 23 heterozygous mice with a P value of 2.945 x 10-21 (Figure 6).  A substantial semidominant effect was also observed (P = 3.348 x 10-5). 

The effect of the mutation at the cDNA and protein levels has not been examined, but the mutation is predicted to result in the use of a cryptic splice site in exon 8. The resulting transcript would have a 78-base pair deletion in exon 8. The deletion would cause an in-frame deletion of 26 amino acids beginning after amino acid 209 of the protein, which is normally 880 amino acids long.

         <--exon 7 exon 8-->                   intron 8-->        exon 9-->
1015 ……GAAAGTTACAA GTACCAGGTACTT……AAGCGGCTAGTT gtaagtctatacggct…… GTGGACTGGACT…… 1122
203  ……-E--S--Y--K --Y--Q--V--L-……-K--R--L--V-                    -V--D--W--T-…… 213
              correct            deleted                             correct
 

The donor splice site of intron 8, which is destroyed by the corpulent mutation, is indicated in blue lettering and the mutated nucleotide is indicated in red.

Illustration of Mutations in
Gene & Protein
Protein Prediction
Figure 7. Domain structure of BBS9. BBS9 has an N-terminal β-propeller fold followed by an amphipathic helical linker, a γ-adaptin ear (GAE) β-sandwich domain, a mixed α/β platform domain, and an α-helical domain. The corpulent mutation occurs within the splice donor site of intron 8.
Figure 8. Crystal structure of the β-propeller region of human BBS9. The WD40 β-propeller has seven anti-parallel, four-stranded β-sheets arranged around a central pore. Image created in UCSF Chimera and based on PDBID 4YD8.

BBS9 (alternatively, PTHB1) has an N-terminal β-propeller fold, an amphipathic helical linker, a γ-adaptin ear (GAE) domain, a mixed α/β platform domain, and a C-terminal α-helical domain [Figure 7; (1;2)].

The BBS9 β-propeller fold is a WD40 β-propeller, which has seven anti-parallel, four-stranded β-sheets arranged around a central pore [Figure 8; PDB:4YD8; (2)]. The BBS9 β-propeller region is proposed to mediate protein-protein interactions with other BBS proteins.  

The GAE domain is approximately 120-amino acids, and functions in the recruitment of accessory proteins that modulate the functions of GAE domain-containing proteins in membrane trafficking. The GAE domain forms an immunoglobulin-like β-sandwich fold composed of eight β-strands and two short α-helices (3-5). The binding site for the accessory proteins is in a shallow hydrophobic trough surrounded by charged (mainly basic) residues (4;5). In clathrin adaptors, the GAE and platform domains are an appendage domain that recruits either regulators of coat assembly or factors that promote correct targeting of the coated vesicle (6).

The corpulent mutation is a splice site mutation that is predicted to result in an in-frame deletion of 26 amino acids after amino acid 209. The deletion is within the β-propeller region.

Expression/Localization

BBS9 is ubiquitously expressed. BBS9 is localized to nonmembranous centriolar satellites in the cytoplasm as well as to the membrane of the cilium (7).

Background

Figure 9. Model of the function of the BBSome. The BBSome organizes IFT-A, IFT-B and kinesin motors into a functional complex. Entry of protein cargo to the cilium is regulated by active forms of Rab8, a master modulator for the ciliary protein trafficking. Rab8 is recruited to the basal body of primary cilium. The activities of Rab8 are then regulated by Rabin8, and its activity and basal body localization is modulated by the BBSome and Rab11. IFT particles dissociate after reaching the ciliary tip, and then the BBSome and DYF-2 coordinate to reorganize the entire IFT complex to ready it for retrograde transport.

Almost every cell in the body carries a single primary cilium at a certain stage of its life cycle. Primary cilia are required for planar cell polarity, phototransduction, olfaction, adipocyte differentiation, and Hedgehog signaling.

Several events occur to facilitate transport of membrane proteins to the cilia: (i) sorting and packaging into carrier vesicles, (ii) docking and fusion of vesicles with the base of the cilium, and (iii) intraflagellar transport (IFT) from the cilia base to cilia tip. Rab8GTP promotes the extension of the ciliary membrane as well as docking and fusion of the carrier vesicles with the base of the cilium (7). The Rab8GTP localizes to the cilium whereby it promotes docking and fusion of vesicles to the base of the ciliary membrane. BBS9 is one of eight proteins (i.e., BBS1, BBS2, BBS4, BBS5, BBS7, BBS8, BBS9, and BBIP10) that form the core of the Bardet-Biedl syndrome protein complex (BBSome) required for ciliogenesis [Figure 9; (7;8)]. The BBSome transports and sorts membrane proteins into the primary cilium and maintains the unique protein composition of the ciliary membrane. The BBSome  associates with the ciliary membrane and binds to the Rab8 guanosyl exchange factor RAB3IP/Rabin8 (7). The BBsome is proposed to direct vesicular trafficking to the cilium by promoting Rabin8-associated activation of Rab8 (7). The exact function of BBS9 within the BBSome is unknown, but it may promote proper BBSome assembly and/or its ciliary localization (7).

Mutations in BBS9 are linked to Bardet-Bidel syndrome-9 [OMIM: #615986; (9-11)]. Patients with Bardet-Bidel syndrome-9 exhibit obesity, polydactyly, renal anomalies, retinopathy, and mental retardation (10). The phenotypes observed in Bardet-Bidel syndrome are proposed to be due to defects in vesicular transport in the cilium. Mutations in BBS9 have also been associated with premature ovarian failure, which causes female infertility (12).

Bbs9-deficient mice have not been generated/phenotypically characterized (MGI; accessed October 23, 2017). However, knockdown of BBS9 in zebrafish resulted in retina and brain developmental defects that were rescued by expression of human BBS9 mRNA (13). Mutations in the other members of the BBSome result in similar phenotypes, including increased body weights/obesity, defects in social function, retinal degeneration, renal cysts, olfaction defects, and male infertility due to cilium dysfunction [Bbs1 (14;15); Bbs2 (16); Bbs4 (15;16); Bbs5 (MGI; accessed October 27, 2017); Bbs7 (17); Bbs8 (18); and Bbip1 (8)].

Putative Mechanism

Cilia-associated obesity is proposed to result from an altered initiation of the ciliary transport at the basal area of the cilium, resulting in changes in adipogenesis, central signaling of food intake, and/or odor perception. During adipogenesis, ciliary proteins have putative roles in maintaining the function of pro-adipogenic factors (e.g., insulin-like growth factor 1 receptor (IGF1-R), CAATT/enhancer binding protein α-β (CEBP/A-B) and PPARγ), while also promoting anti-adipogenic signaling pathways (e.g., SHH, Wnt, and Notch) (19). The cilium may also function in neuronal food intake regulation. Cilia of hypothalamic neurons contain receptors involved in food intake regulation, including orexigenic neuropeptide Y (see the record oleo for information about neuropeptide Y receptor [NPY5R]) (20) and melanin-concentrating hormone receptor 1 (Mchr1) (21). A reduction in the number of cilia or a change in length may induce hyperphagia and obesity. Hyposmia and anosmia (reduction or ablation of odor perception, respectively) may lead to weight loss, weight gain or no change in weight. Food odor can stimulate appetite, food-seeking behavior and food ingestion, but strong or prolonged exposure to a food odor can have a satiating effect (22). The phenotype of the corpulent mice indicates loss of BBS9corpulent function.

Primers PCR Primer
corpulent_pcr_F: CACAGCATTTATGGACTGTGTAC
corpulent_pcr_R: TTCAATGAAAAGTGGTGCCAC

Sequencing Primer
corpulent_seq_F: CTTGCCTTTGCCACAGAT
corpulent_seq_R: GCCATTTACAGTTGCAGCACG
Genotyping

PCR program

1) 94°C 2:00
2) 94°C 0:30
3) 55°C 0:30
4) 72°C 1:00
5) repeat steps (2-4) 40x
6) 72°C 10:00
7) 4°C hold


The following sequence of 410 nucleotides is amplified (chromosome 9, + strand):


1   cacagcattt atggactgtg tacctttaaa gtttaaagtg atttttactt taaaacttat
61  gatgtcactt atcctttgcc atctcttcaa acaggtacca ggtacttgcc tttgccacag
121 atgcagacaa gaagcaggag atggaacagc aaaaacttgg ttctggaaag cggctagttg
181 taagtctata cggctcagga acgctactca tggtcacagt tgtagagaat tcctggaaga
241 aaattaaaaa tgtattttga caatttctct ttatccgctg tacaatgacc agggattgaa
301 ttgcctgttg ctgtgctgta gatcatgtct gtccacatgc cgaaatgcat cagaatacag
361 cgtggcgtgc tgcaactgta aatggctggg tggcaccact tttcattgaa 


Primer binding sites are underlined and the sequencing primers are highlighted; the mutated nucleotide is shown in red.

References
Science Writers Anne Murray
Illustrators Diantha La Vine
AuthorsEmre Turer, Zhao Zhang, Carol Wise, Jonathan Rios, and Bruce Beutler