Micrognathia, glossoptosis, and cleft palate comprise one of the most common malformation sequences, Robin sequence. It is a component of the TARP syndrome, talipes equinovarus, atrial septal defect, Robin sequence, and persistent left superior vena cava. This disorder is X-linked and severe, with apparently 100% pre- or postnatal lethality in affected males. Here we characterize a second family with TARP syndrome, confirm linkage to Xp11.23-q13.3, perform massively parallel sequencing of X chromosome exons, filter the results via a number of criteria including the linkage region, use a unique algorithm to characterize sequence changes, and show that TARP syndrome is caused by mutations in the RBM10 gene, which encodes RNA binding motif 10.
We further show that this previously uncharacterized gene is expressed in midgestation mouse embryos in the branchial arches and limbs, consistent with the human phenotype. We conclude that massively parallel sequencing is useful to characterize large candidate linkage intervals and that it can be used successfully to allow identification of disease-causing gene mutations. Introduction Advances in genomic technologies can markedly speed the analysis of genetic contributions to disease and make once-intractable questions tractable.
The technique of massively parallel sequencing with exon capture allows rapid assessment of mutations that cause human diseases. The syndrome we have chosen to study with this new technology is an X-linked pleiotropic developmental anomaly syndrome (MIM ) comprising micrognathia, glossoptosis, and cleft palate (currently described as Robin sequence), persistent left superior vena cava, atrial septal defect, and talipes equinovarus, which was originally called “Robin's syndrome” by Gorlin et al.
Download Ebook The Best Goodbye (Abbi Glines) PDF, EPUB, MOBI.
Subsequently, this disorder was designated as TARP syndrome based on the acronym formed by talipes equinovarus, atrial septal defect, Robin sequence, and persistence of the left superior vena cava (to avoid confusion of Robin sequence and “Robin's syndrome”). Linkage analysis of the family originally evaluated by Gorlin et al. Confirmed the X-linked inheritance by mapping the locus to Xp11.23-q13.3. Although this mapping confirmed the inheritance pattern and narrowed the candidate gene list, the region was nearly 28 Mb (from 46.42 Mb to 74.04 Mb; genome build 36) and contained more than 200 genes or transcripts, including several complex gene families of G antigen transcripts. Thus, the task of interrogating all of these genes for sequence variants and validating the variants was daunting. Here we report the identification of the second family manifesting this syndrome, describe the phenotype in three affected males, show that the second family has linkage that is consistent with that of the initial family, describe our use of chromosome X exon target capture and massively parallel (so-called “next-gen”) sequencing to identify RBM10 (MIM ) mutations in the two families, and show that the murine Rbm10 gene is expressed in midgestation embryos in a pattern consistent with the human phenotype. Target Selection and Sequencing DNA isolation, genotyping, and haplotype analysis were performed as previously described.
To sequence genes, we used solution hybridization selection (SureSelect, Agilent) to generate a single-end sequencing library according to the manufacturer's directions (Illumina). For the purpose of this project, the exon targets were defined as all coding DNA sequence (UCSC) gene exons between positions 2,710,679 bp and 154,500,000 bp, which is from XG (MIM ) to TMLHE (MIM ), inclusive, human genome build 36 (hg18). The target sequence was 2,784,426 bp, and the oligonucleotide library was designed to target 2,264,175 bp of this (81.3%).
Briefly, 3 μg of genomic DNA was sheared via acoustic fragmentation (model S2, Covaris) with the following parameters: duty cycle, 10%; intensity, 5; cycle burst, 200; time, 180 s; set mode frequency, sweeping; temperature, 4°C. The DNA was subjected to end repair and adaptor ligation (NEBNext DNA Sample Prep Reagent Set 1, New England Biolabs). Size selection of 200–300 bp fragments was performed by using electrophoresis (4% NuSieve GTG gel, Lonza Group).
The prepped library was amplified with primers 1.1 and 2.1 (Illumina). The amplified library was hybridized to biotinylated RNA library baits (Agilent), and targeted sequences were purified with magnetic beads (Dynal M-280 streptavidin, Invitrogen). Purified target sequences were amplified with proprietary PCR primers (Agilent) and DNA polymerase (Herculase II fusion, Stratagene). Interpretation of Coding Variations To predict possible deleterious effects of coding changes on the native state of the protein, we developed novel software to predict the deviation from the “normal” function relative to the human reference sequence and associated UCSC known gene models (hg18). This software computes the consequence of the change to the sequence and predicts whether the nucleotide variation leads to a silent, missense, nonsense, canonical splice-site, or frameshift change. CDPred, a module within the coding analysis software, is an algorithm that scores amino acid variants on the basis of evolutionary conservation in conserved protein coding domains. CDPred assigns all nonsense, splice-site, and frameshift changes in a protein the most damaging score (−30), under the hypothesis that these variants lead to truncated protein products or nonsense-mediated mRNA decay, both of which can be equivalent to a null protein.
This algorithm also predicts whether a missense change would be deleterious to the function of the protein, based on a log-likelihood score computed from a position-specific scoring matrix of the best-aligning protein-conserved domain model to the query protein. For a given missense change (relative to the reference codon), the scores are in a range of +20 to −20. The CDPred algorithm is initialized with a reference protein sequence and the variant (missense or nonsense) relative to the reference. An alignment of the reference protein sequence to a database of protein domain models (position-specific scoring matrices) was computed with RPS-BLAST (parameters: –e 0.00001 and –F T).
Sequence alignments with greater than 80% overlap of the length of the model with an E value less than 1e-5 were considered to be significant. The delta log-likelihood scores were computed as the difference in position-specific scores of the variant amino acid and the reference amino acid. The lower or more negative the score, the more deleterious was the predicted change. Positive scores above 3 may also be potentially damaging, but a strongly positive score generally results from a situation in which the normal human amino acid at a position is different from many or all other aligned species over a domain and in which the variant allele is closer to the ancestral form. For example, if an individual had the ancestral form of the FOXP2 gene (MIM ), that individual would potentially exhibit loss of speech; however, the delta score for the amino acid positions that are unique to all other humans would be positive because that individual more closely resembles the ancestral form. The CDPred software is available for evaluation and implementation at the website listed in the section of the paper. Mutation Confirmation and Segregation Studies PCR amplification and Sanger sequencing of DNA samples to confirm the mutations and cosegregation were performed as previously described.
PCR primer sequences are available upon request. Mutation numbering was performed according to Human Gene Variation Society nomenclature with reference sequence. Human studies were performed according to an approved human subjects research protocol of the NHGRI institutional review board, and informed consent was obtained. In Situ Analysis C57BL/6J mouse embryos were fixed overnight in 4% paraformaldehyde in PBS. Reverse-transcribed digoxigenin-conjugated probes were made from PCR products with polymerase binding site linkers. The following DNA source was used for probe synthesis: Rbm10, reverse transcription PCR from C57BL/6J skin total RNA (Rbm10-F-T3- GCGCGAATTAACCCTCACTAAAGGGCGGGATGGATTAGGCAGTGAC; Rbm10-R-T7-GCGCGTAATACGACTCACTATAGGGACAAAAGGAACATGATTTGAG). In situ hybridizations were performed by using published protocols with the following modifications.
After probe hybridization, ribonuclease A digestion was omitted, and tris-buffered saline was used in place of PBS. BM-purple substrate (Roche, Molecular Biochemicals) was used in place of 5-bromo-4-chloro-3-indolyl phosphate and nitroblue tetrazolium. Multiple embryos were evaluated at each stage of development analyzed: E9.5 (n = 6), E10.5 (n = 5), and E11.5 (n = 4). Murine procedures were performed in accordance with NIH guidelines and under NHGRI mouse protocol G94-7. Clinical Reports Individual III-3 () was born at 32 weeks gestation to a 22-year-old gravida 4 para 2–3 abortus 1 mother.
The pregnancy was complicated by oligohydramnios and premature rupture of the membranes. Delivery was breech. Although a fetal heartbeat was noted on examination within 2 min prior to delivery, APGAR scores were 0 at 1, 5, and 10 min. The baby boy died at 5 min of age. Upon postmortem examination, his weight was 1335 g and his crown heel length was 41 cm (both appropriate for gestational age). External examination showed apparently low-set ears, micrognathia, and a large cleft of the hard and soft palate.
The hands were broad, with ulnar deviation. Bilateral talipes equinovarus and rocker bottom feet were present. The right testicle was undescended.
Upon internal examination, he was found to have a large atrial septal defect, and his lung showed an underdeveloped third lobe with overall marked underdevelopment of alveoli. A Meckel diverticulum was present. Radiography showed a severely underdeveloped mandible. Microscopic description showed extramedullary hematopoiesis of the liver. Standard resolution karyotype was reported as normal, 46,XY. Pedigree of TARP Syndrome Family 2 Individual III-9 was born at 37.5 weeks gestation to a 26-year-old gravida 1 para 1 mother. The prenatal course was notable for intrauterine growth retardation late in pregnancy.
At birth, the baby boy weighed 2090 g (∼tenth centile) and his length was 46 cm (appropriate for gestational age). He was apneic, hypotonic, and cyanotic post delivery. He had micrognathia, cleft palate, glossoptosis, and hyaline membrane disease. He developed seizures at 24 hr. Head ultrasound showed basal ganglia hemorrhages and subdural hematoma. He was intubated, but he died at 8 days of age with liver failure, kidney failure, hyaline membrane disease, athetoid movements, and seizures. Autopsy was declined.
Individual III-15 was born at 28 weeks gestation to a 23-year-old gravida 2 para 1–2 mother. Few clinical details are available about this affected baby boy. At birth, he was found to have multiple anomalies including cleft palate (Robin sequence), possible cardiac defects, and bilateral talipes equinovarus. He died of his multiple congenital anomalies.
Results Linkage analysis of the original family identified an 11 cM region in Xp11.23-q13.3 with a peak LOD score of 2.75 at marker DXS1039. The borders of the region were defined by recombinant markers at DXS1003 and DXS8092. Haplotype analysis was performed in family 2 with these three markers and additional markers both within and outside of the region ( DXS8054, DXS1208, DXS7132, and DXS6800).
Haplotypes were consistent with affection or carrier status in all individuals tested, and there were no recombinations in family 2 that narrowed the region defined by family 1 (data not shown). The linked short tandem repeat polymorphism haplotype in the two families was distinct, suggesting that the mutations in the two families were likely to be distinct. The target-selected DNA libraries from one female heterozygote from each of the two families (females were used for optimal DNA quantity and quality) were sequenced on one lane each of a sequencing instrument (Illumina GAII) in single-end 36 bp configuration, which yielded 20,262,045; 18,775,942 reads (family 1; family 2) or 729,433,620; 675,933,912 bp of total sequence. Of this sequence, 43.8%; 45.2% could be uniquely aligned to the targeted exons.
This aligned sequence yielded a gross overall coverage of 115×; 110× of the target. The capture efficiency varied across the targets with 2,239,228; 2,234,963 bp (80.5%) of the target with ≥1× coverage, 2,136,202; 2,128,057 bp (76.5%) with ≥10× coverage, and 2,071,297; 2,059,012 bp (74.1%) with ≥20× coverage. The most probable genotype variant-calling software (J.K.T., N.F.H., J.C.M., and L.G.B.
Et al., unpublished data) was able to make high-confidence genotype calls on 1,956,070; 1,941,688 bp of this sequence (70.0%). We filtered the results for variants in genes on the basis of several attributes that we reasoned were appropriate for this disorder (). Heterozygosity was used because the test subjects were both unaffected female carriers for an X-linked trait. We reasoned that the variant should be severe because it caused a lethal phenotype in hemizygous males, so we filtered for nonsynonymous and stop alleles. Similarly, we also scanned the sequences for insertion or deletion variants, reasoning that most of these would cause frameshifting, null alleles. The criterion of novelty was applied in two ways. First, we filtered for variants not present in dbSNP, reasoning that a variant causing a rare phenotype should not be common.
Second, we filtered for the absence of a sequence variant in three control DNA samples. The controls included two samples from males with syndromic microphthalmia and a single sample from a female parent of a child with a previously uncharacterized X-linked lethal disorder. We also applied a filter that bounded the variants genomically within the defined linkage interval from the original family, which was confirmed in family 2.
The Number of Genes with One or More Variants Following Each Filtering Criterion The initial analysis focused on single base pair substitutions, which, with the application of the heterozygosity + nonsynonymous + nonsense + novelty filters (), showed a single nonsense mutation, c.1235G>A in RBM10 (RNA binding motif 10), which predicts p.Trp412X in family 2 (, lower left panels; ). No other nonsense mutations were identified in either family with this alignment and these filters.
A second analysis for small deletions and insertions was performed with the following filters: nonsynonymous + novelty + frameshifting. This allowed identification of a single mutation in family 1, which was c.1893_1894insA in RBM10, which predicts p.Pro632ThrfsX41 (, lower right panels; ). Each RBM10 variant was confirmed by PCR amplification and Sanger sequencing (, bottom) in multiple individuals within the respective family, and mutation status segregated with known carrier status in all individuals (cosegregation is shown for family 2 in; cosegregation for family 1 is not shown). We concluded that these were causative variants on the basis of the co-occurrence of a nonsense and a frameshift mutation in the same gene in two families with a highly similar and extremely rare phenotype.
Cartoon of the RBM10 Gene Structure, Alternative Splicing, Conserved Domains, and Mutations Found in Two Families with TARP Syndrome Because the RBM10 gene product was poorly characterized and mutations had not been described with any abnormal phenotype in humans or in animal models, we set out to characterize the expression of this gene in midgestation mouse embryos. The orthologous murine Rbm10 isoform 1 gene product is 96% identical to the human protein. We isolated an Rbm10 antisense probe from the mouse sequence via RT-PCR of whole mouse RNA isolated from embryonic day E16.5 skin.
Whole-mount in situ expression analysis of the murine Rbm10 gene at E9.5 and E10.5 of development showed a similar pattern of expression, with the most robust staining observed in the first branchial arch (which gives rise to the mandible), second branchial arch, developing limb buds, and tailbud (). Robust expression remained for E11.5 embryos in both the limb and tail bud regions, whereas expression in branchial arches one and two decreased at this stage (data not shown). This pattern of expression correlated well with the human malformations observed in TARP syndrome, which include severe micrognathia and limb defects. Discussion Here we show that the use of targeted exon capture with filters that included a genomic delimitation defined by a linkage region, the zygosity state, novelty, a severe deleterious mutation prediction, and a common mutation among two affected families can identify a gene mutated in a rare disorder. This approach was efficient and less time consuming compared with what would be required for PCR amplification and Sanger sequencing of the approximately 200 genes within the ∼28 Mb candidate interval. These data show two rare variants in RBM10, one in each of two families with a highly similar, rare, pleiotropic multiple congenital anomaly syndrome. We suggest that the delineation of nonsense or frameshift mutations in each of two families with a rare disorder is unlikely by chance alone.
A simple estimate of the probability of finding these variants in the same gene in two patients in a candidate gene region of 200 genes is 1/200, or 0.005. Notably, nonsense or stop variants in RBM10 were not found in the recent report of 208 patients with mental retardation.
In addition, the mouse expression data showed a striking correlation with the phenotypic manifestations in the human disorder. Finally, the variants cosegregated with the carrier status of the females in both families. We conclude that these data show that nonsense or null mutations in RBM10 cause syndromic Robin sequence, or TARP syndrome. The only previous report of RBM10 and human phenotypes that we identified was a case report of a girl with an X;15 translocation with hypertelorism, a small face, high forehead, small, low-set ears, ulnar deviation of the hands, agenesis of the corpus callosum, and a hypoplastic fifth toe and metatarsal. This X;15 translocation breakpoint was near RBM10, although it was not determined whether the expression of RBM10 was altered by this translocation.
This girl and the boys from families 1 and 2 with TARP did not have extensive phenotypic similarity; thus, it is unlikely that this translocation mediated its phenotypic effects primarily through disruption or dysregulation of RBM10. The RBM10 gene and its 930 amino acid protein product is a member of the RNA binding motif (RBM) gene family. The RBM gene family is large, but mutations in only a few of the family members cause recognized human disorders, including dilated cardiomyopathy (MIM ), which is caused by mutations in RBM20 (MIM ), and alopecia, neurologic defects, and endocrinopathy syndrome (MIM ), which are caused by mutations in RBM28 (MIM ). A number of RBM genes have been shown to be important for RNA processing, RNA splicing, apoptosis, and other diverse biologic roles. RBM10 has been shown to undergo typical X chromosome inactivation.
RBM10 is predicted to include a zinc finger motif, a G patch, and two RNA recognition motif (RRM) domains (). This architecture is found in a variety of RNA binding proteins, including different heterogeneous nuclear ribonucleoproteins (RNPs), protein components of small nuclear RNPs, and those implicated in the regulation of alternative splicing. A typical RRM domain has two RNA or DNA binding sites and one RRM dimerization site.
The RBM10 gene has two known alternative splice forms. Variant 1 codes for a 930 amino acid protein (, 23 coding exons), and variant 2 codes for an 852 amino acid protein (, 22 coding exons). The nonsense change (p.Trp412X in the long form) would truncate the protein just after the end of the second RRM motif, which may destabilize the overall structure of the RBM10 protein. The frameshift (p.Pro632ThrfsX41) would also affect the last third of the protein. Alternatively, these truncating mutations may induce nonsense-mediated mRNA decay. In either case, the loss of function of RBM10 in TARP syndrome demonstrates that this gene is critical for normal mammalian development.
Further elucidation of the role of this gene in early embryogenesis will lead to a better understanding of congenital anomalies that affect the face, heart, and limbs. The filtering of the sequence data shows that the RBM10 gene could have been identified without linkage data, because the two RBM10 mutations were the only variants that met all of the filtering criteria for the entire targeted set of exons in these two families. We recognize that the determination of appropriate filtering criteria is subjective and that the outcome may be difficult to predict. For example, our criterion that the mutations should be frameshifting or nonsense would not hold for all X-linked male lethal disorders. Our overall approach was to begin with stringent or conservative filters and relax them in successive analyses in an attempt to identify the mutated gene.
We were fortunate that the first filtering screen was successful, but we recognize that this will not always be the case. As noted in the methods section, this sequence experiment included five samples. The others were three patients who, among them, had two disorders distinct from TARP. We believe that we have identified the causative mutations in one of these two disorders by use of the same filtering used to identify the TARP mutations, but the other has failed to identify a causative gene mutation (data not shown). Much more data on similar projects will need to be generated to develop optimal approaches to filtering. Another potential cause of a failure to identify causative mutations by exon capture is oligonucleotide design and sequencing coverage.
This implementation of solution hybridization exon targeting included oligonucleotide designs for just over 80% of the exonic base pairs. After generating >100× of aligned sequence coverage for each sample, we had adequate coverage (≥10 quality base score depth) for detection of heterozygous changes in 80% of the base pairs of RBM10. This would suggest that for any pair of samples from patients with TARP syndrome, where each sample has one causative mutation in each of two distinct targets, there is about a 60%–65% chance that the mutations in both will be detected. In addition, nonexonic variants, synonymous variants that may affect splice enhancers, and other types of causative mutations are likely to be challenging to detect by exon targeting. Therefore, early successes in gene identification by target selection will not be representative of all disorders, and it will be important to develop new approaches to address these challenges.
Narrated by the central character Prue Sarn, whose life is blighted by having a harelip. Only the weaver, Kester Woodseaves, perceives her inner beauty but Prue cannot believe herself worthy of him. Prue is wrongly accused of murder and only one man can save her and take her away to the happiness she believes she can never possess because of her harelip. A forgotten classic set in rural Shropshire at the turn of the 19th century blends a simple, rustic love story with a profound sense of nature’s mystic truth. Prue Sarn is an original and appealing heroine of English literature as she triumphs over a physical handicap to win her heart’s desire. Skillfully woven through this story is the aura of the English countryside, its flora and fauna anticipating every turn of the plot.
396 pages, with a reading time of ~6.25 hours (99,009 words), and first published in 1924. This DRM-Free edition published by epubBooks, 2011. Excerpt It was at a love-spinning that I saw Kester first.
And if, in these new-fangled days, when strange inventions crowd upon us, when I hear tell there is even a machine coming into use in some parts of the country for reaping and mowing, if those that may happen will read this don’t know what a love-spinning was, they shall hear in good time. But though it was Jancis Beguildy’s love-spinning, she being three-and-twenty at that time and I being two years less, yet that is not the beginning of the story I have set out to tell. Kester says that all tales, true tales or romancings, go farther back than the days of the child; aye, farther even than the little babe in its cot of rushes. Maybe you never slept in a cot of rushes; but all of us did at Sarn.
There is such a plenty of rushes at Sarn, and old Beguildy’s missus was a great one for plaiting them on rounded barrel-hoops. Then they’d be set on rockers, and a nice clean cradle they made, soft and green, so that the babe could feel as big-sorted as a little caterpillar (painted butterflies-as-is-to-be, Kester calls them) sleeping in its cocoon. Kester’s very set about such things. Never will he say caterpillars. He’ll say, ‘There’s a lot of butter-flies-as-is-to-be on our cabbages, Prue.’ He won’t say ‘It’s winter.’ He’ll say, ‘Summer’s sleeping.’ And there’s no bud little enough nor sad-coloured enough for Kester not to callen it the beginnings of the blow. But the time is not yet come for speaking of Kester.
It is the story of us all at Sarn, of Mother and Gideon and me, and Jancis (that was so beautiful), and Wizard Beguildy, and the two or three other folk that lived in those parts, that I did set out to tell. There were but a few, and maybe always will be, for there’s a discouragement about the place. It may be the water lapping, year in and year out - everywhere you look and listen, water; or the big trees waiting and considering on your right hand and on your left; or the unbreathing quiet of the place, as if it was created but an hour gone, and not created for us. Or it may be that the soil is very poor and marshy, with little nature or goodness in the grass, which is ever so where reeds and rushes grow in plenty, and the flower of the paigle. Happen you call it cowslip, but we always named it the paigle, or keys of heaven. It was a wonderful thing to see our meadows at Sarn when the cowslip was in blow.
Gold-over they were, so that you would think not even an angel’s feet were good enough to walk there. You could make a tossy-ball before a thrush had gone over his song twice, for you’d only got to sit down and gather with both hands.
Every way you looked there was naught but gold, saving towards Sarn, where the woods began, and the great stretch of grey water, gleaming and wincing in the sun. Neither woods nor water looked darksome in that fine spring weather, with the leaves coming new, and buds the colour of corn in the birch-tops.
Only in our oak wood there was always a look of the back-end of the year, their young leaves being so brown. So there was always a breath of October in our May.
But it was a pleasant thing to sit in the meadows and look away to the far hills. The larches spired up in their quick green, and the cowslip gold seemed to get into your heart, and even Sarn Mere was nothing but a blue mist in a yellow mist of birch-tops. And there was such a dream on the place that if a wild bee came by, let alone a bumble, it startled you like a shout.
If a bee comes in at the window now to my jar of gillyflowers, I can see it all in clear colours, with Plash lying under the sunset, beyond the woods, looking like a jagged piece of bottle glass. Plash Mere was bigger than Sarn, and there wasn’t a tree by it, so where there were no hills beyond it you could see the clouds rooted in it on the far side, and I used to think they looked like the white water-lilies that lay round the margins of Sarn half the summer through. There was nothing about plash that was different from any other lake or pool. There was no troubling of the waters, as at Sarn, nor any village sounding its bells beneath the furthest deeps. It was true, what folks said of Sarn, that there was summat to be felt there. It was at Plash that the Beguildys lived, and it was at their dwelling, that was part stone house and part cave, that I got my book learning. It may seem a strange thing to you that a woman of my humble station should be able to write and spell, and put all these things into a book.
And indeed when I was a young wench there were not many great ladies, even, that could do much more scribing than to write a love-letter, and some could but just write such things as ‘This be quince and apple’ on their jellies, and others had ado to put their names in the marriage register. Many have come to me, time and again, to write their love-letters for them, and a bitter old task it is, to write other women’s love-letters out of your own burning heart.