Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits.

Bumblebees are a diverse group of globally important pollinators in natural ecosystems and for agricultural food production. With both eusocial and solitary life-cycle phases, and some social parasite species, they are especially interesting models to understand social evolution, behavior, and ecology. Reports of many species in decline point to pathogen transmission, habitat loss, pesticide usage, and global climate change, as interconnected causes. These threats to bumblebee diversity make our reliance on a handful of well-studied species for agricultural pollination particularly precarious. To broadly sample bumblebee genomic and phenotypic diversity, we de novo sequenced and assembled the genomes of 17 species, representing all 15 subgenera, producing the first genus-wide quantification of genetic and genomic variation potentially underlying key ecological and behavioral traits. The species phylogeny resolves subgenera relationships, whereas incomplete lineage sorting likely drives high levels of gene tree discordance. Five chromosome-level assemblies show a stable 18-chromosome karyotype, with major rearrangements creating 25 chromosomes in social parasites. Differential transposable element activity drives changes in genome sizes, with putative domestications of repetitive sequences influencing gene coding and regulatory potential. Dynamically evolving gene families and signatures of positive selection point to genus-wide variation in processes linked to foraging, diet and metabolism, immunity and detoxification, as well as adaptations for life at high altitudes. Our study reveals how bumblebee genes and genomes have evolved across the Bombus phylogeny and identifies variations potentially linked to key ecological and behavioral traits of these important pollinators.

Improved Pathogenic Variant Localization via a Hierarchical Model of Sub-regional Intolerance.

Different parts of a gene can be of differential importance to development and health. This regional heterogeneity is also apparent in the distribution of disease-associated mutations, which often cluster in particular regions of disease-associated genes. The ability to precisely estimate functionally important sub-regions of genes will be key in correctly deciphering relationships between genetic variation and disease. Previous methods have had some success using standing human variation to characterize this variability in importance by measuring sub-regional intolerance, i.e., the depletion in functional variation from expectation within a given region of a gene. However, the ability to precisely estimate local intolerance was restricted by the fact that only information within a given sub-region is used, leading to instability in local estimates, especially for small regions. We show that borrowing information across regions using a Bayesian hierarchical model stabilizes estimates, leading to lower variability and improved predictive utility. Specifically, our approach more effectively identifies regions enriched for ClinVar pathogenic variants. We also identify significant correlations between sub-region intolerance and the distribution of pathogenic variation in disease-associated genes, with AUCs for classifying de novo missense variants in Online Mendelian Inheritance in Man (OMIM) genes of up to 0.86 using exonic sub-regions and 0.91 using sub-regions defined by protein domains. This result immediately suggests that considering the intolerance of regions in which variants are found may improve diagnostic interpretation. We also illustrate the utility of integrating regional intolerance into gene-level disease association tests with a study of known disease-associated genes for epileptic encephalopathy.

The Encyclopedia of DNA elements (ENCODE): data portal update.

The Encyclopedia of DNA Elements (ENCODE) Data Coordinating Center has developed the ENCODE Portal database and website as the source for the data and metadata generated by the ENCODE Consortium. Two principles have motivated the design. First, experimental protocols, analytical procedures and the data themselves should be made publicly accessible through a coherent, web-based search and download interface. Second, the same interface should serve carefully curated metadata that record the provenance of the data and justify its interpretation in biological terms. Since its initial release in 2013 and in response to recommendations from consortium members and the wider community of scientists who use the Portal to access ENCODE data, the Portal has been regularly updated to better reflect these design principles. Here we report on these updates, including results from new experiments, uniformly-processed data from other projects, new visualization tools and more comprehensive metadata to describe experiments and analyses. Additionally, the Portal is now home to meta(data) from related projects including Genomics of Gene Regulation, Roadmap Epigenome Project, Model organism ENCODE (modENCODE) and modERN. The Portal now makes available over 13000 datasets and their accompanying metadata and can be accessed at: https://www.encodeproject.org/.

Novel family of terpene synthases evolved from trans-isoprenyl diphosphate synthases in a flea beetle.

Sesquiterpenes play important roles in insect communication, for example as pheromones. However, no sesquiterpene synthases, the enzymes involved in construction of the basic carbon skeleton, have been identified in insects to date. We investigated the biosynthesis of the sesquiterpene (6R,7S)-himachala-9,11-diene in the crucifer flea beetle Phyllotreta striolata, a compound previously identified as a male-produced aggregation pheromone in several Phyllotreta species. A (6R,7S)-himachala-9,11-diene-producing sesquiterpene synthase activity was detected in crude beetle protein extracts, but only when (Z,E)-farnesyl diphosphate [(Z,E)-FPP] was offered as a substrate. No sequences resembling sesquiterpene synthases from plants, fungi, or bacteria were found in the P. striolata transcriptome, but we identified nine divergent putative trans-isoprenyl diphosphate synthase (trans-IDS) transcripts. Four of these putative trans-IDSs exhibited terpene synthase (TPS) activity when heterologously expressed. Recombinant PsTPS1 converted (Z,E)-FPP to (6R,7S)-himachala-9,11-diene and other sesquiterpenes observed in beetle extracts. RNAi-mediated knockdown of PsTPS1 mRNA in P. striolata males led to reduced emission of aggregation pheromone, confirming a significant role of PsTPS1 in pheromone biosynthesis. Two expressed enzymes showed genuine IDS activity, with PsIDS1 synthesizing (E,E)-FPP, whereas PsIDS3 produced neryl diphosphate, (Z,Z)-FPP, and (Z,E)-FPP. In a phylogenetic analysis, the PsTPS enzymes and PsIDS3 were clearly separated from a clade of known coleopteran trans-IDS enzymes including PsIDS1 and PsIDS2. However, the exon-intron structures of IDS and TPS genes in P. striolata are conserved, suggesting that this TPS gene family evolved from trans-IDS ancestors.

DVL1 frameshift mutations clustering in the penultimate exon cause autosomal-dominant Robinow syndrome.

Robinow syndrome is a genetically heterogeneous disorder characterized by mesomelic limb shortening, genital hypoplasia, and distinctive facial features and for which both autosomal-recessive and autosomal-dominant inheritance patterns have been described. Causative variants in the non-canonical signaling gene WNT5A underlie a subset of autosomal-dominant Robinow syndrome (DRS) cases, but most individuals with DRS remain without a molecular diagnosis. We performed whole-exome sequencing in four unrelated DRS-affected individuals without coding mutations in WNT5A and found heterozygous DVL1 exon 14 mutations in three of them. Targeted Sanger sequencing in additional subjects with DRS uncovered DVL1 exon 14 mutations in five individuals, including a pair of monozygotic twins. In total, six distinct frameshift mutations were found in eight subjects, and all were heterozygous truncating variants within the penultimate exon of DVL1. In five families in which samples from unaffected parents were available, the variants were demonstrated to represent de novo mutations. All variant alleles are predicted to result in a premature termination codon within the last exon, escape nonsense-mediated decay (NMD), and most likely generate a C-terminally truncated protein with a distinct -1 reading-frame terminus. Study of the transcripts extracted from affected subjects' leukocytes confirmed expression of both wild-type and variant alleles, supporting the hypothesis that mutant mRNA escapes NMD. Genomic variants identified in our study suggest that truncation of the C-terminal domain of DVL1, a protein hypothesized to have a downstream role in the Wnt-5a non-canonical pathway, is a common cause of DRS.

COQ4 mutations cause a broad spectrum of mitochondrial disorders associated with CoQ10 deficiency.

Primary coenzyme Q10 (CoQ10) deficiencies are rare, clinically heterogeneous disorders caused by mutations in several genes encoding proteins involved in CoQ10 biosynthesis. CoQ10 is an essential component of the electron transport chain (ETC), where it shuttles electrons from complex I or II to complex III. By whole-exome sequencing, we identified five individuals carrying biallelic mutations in COQ4. The precise function of human COQ4 is not known, but it seems to play a structural role in stabilizing a multiheteromeric complex that contains most of the CoQ10 biosynthetic enzymes. The clinical phenotypes of the five subjects varied widely, but four had a prenatal or perinatal onset with early fatal outcome. Two unrelated individuals presented with severe hypotonia, bradycardia, respiratory insufficiency, and heart failure; two sisters showed antenatal cerebellar hypoplasia, neonatal respiratory-distress syndrome, and epileptic encephalopathy. The fifth subject had an early-onset but slowly progressive clinical course dominated by neurological deterioration with hardly any involvement of other organs. All available specimens from affected subjects showed reduced amounts of CoQ10 and often displayed a decrease in CoQ10-dependent ETC complex activities. The pathogenic role of all identified mutations was experimentally validated in a recombinant yeast model; oxidative growth, strongly impaired in strains lacking COQ4, was corrected by expression of human wild-type COQ4 cDNA but failed to be corrected by expression of COQ4 cDNAs with any of the mutations identified in affected subjects. COQ4 mutations are responsible for early-onset mitochondrial diseases with heterogeneous clinical presentations and associated with CoQ10 deficiency.

Mutations in the gene encoding the E2 conjugating enzyme UBE2T cause Fanconi anemia.

Fanconi anemia (FA) is a rare genetic disorder characterized by genome instability, increased cancer susceptibility, progressive bone marrow failure (BMF), and various developmental abnormalities resulting from the defective FA pathway. FA is caused by mutations in genes that mediate repair processes of interstrand crosslinks and/or DNA adducts generated by endogenous aldehydes. The UBE2T E2 ubiquitin conjugating enzyme acts in FANCD2/FANCI monoubiquitination, a critical event in the pathway. Here we identified two unrelated FA-affected individuals, each harboring biallelic mutations in UBE2T. They both produced a defective UBE2T protein with the same missense alteration (p.Gln2Glu) that abolished FANCD2 monoubiquitination and interaction with FANCL. We suggest this FA complementation group be named FA-T.

Mutations in TUBGCP4 alter microtubule organization via the γ-tubulin ring complex in autosomal-recessive microcephaly with chorioretinopathy.

We have identified TUBGCP4 variants in individuals with autosomal-recessive microcephaly and chorioretinopathy. Whole-exome sequencing performed on one family with two affected siblings and independently on another family with one affected child revealed compound-heterozygous mutations in TUBGCP4. Subsequent Sanger sequencing was performed on a panel of individuals from 12 French families affected by microcephaly and ophthalmic manifestations, and one other individual was identified with compound-heterozygous mutations in TUBGCP4. One synonymous variant was common to all three families and was shown to induce exon skipping; the other mutations were frameshift mutations and a deletion. TUBGCP4 encodes γ-tubulin complex protein 4, a component belonging to the γ-tubulin ring complex (γ-TuRC) and known to regulate the nucleation and organization of microtubules. Functional analysis of individual fibroblasts disclosed reduced levels of the γ-TuRC, altered nucleation and organization of microtubules, abnormal nuclear shape, and aneuploidy. Moreover, zebrafish treated with morpholinos against tubgcp4 were found to have reduced head volume and eye developmental anomalies with chorioretinal dysplasia. In summary, the identification of TUBGCP4 mutations in individuals with microcephaly and a spectrum of anomalies in eye development, particularly photoreceptor anomalies, provides evidence of an important role for the γ-TuRC in brain and eye development.

Submicroscopic deletions at 13q32.1 cause congenital microcoria.

Congenital microcoria (MCOR) is a rare autosomal-dominant disorder characterized by inability of the iris to dilate owing to absence of dilator pupillae muscle. So far, a dozen MCOR-affected families have been reported worldwide. By using whole-genome oligonucleotide array CGH, we have identified deletions at 13q32.1 segregating with MCOR in six families originating from France, Japan, and Mexico. Breakpoint sequence analyses showed nonrecurrent deletions in 5/6 families. The deletions varied from 35 kbp to 80 kbp in size, but invariably encompassed or interrupted only two genes: TGDS encoding the TDP-glucose 4,6-dehydratase and GPR180 encoding the G protein-coupled receptor 180, also known as intimal thickness-related receptor (ITR). Unlike TGDS which has no known function in muscle cells, GPR180 is involved in the regulation of smooth muscle cell growth. The identification of a null GPR180 mutation segregating over two generations with iridocorneal angle dysgenesis, which can be regarded as a MCOR endophenotype, is consistent with the view that deletions of this gene, with or without the loss of elements regulating the expression of neighboring genes, are the cause of MCOR.

Heterozygous reelin mutations cause autosomal-dominant lateral temporal epilepsy.

Autosomal-dominant lateral temporal epilepsy (ADLTE) is a genetic epilepsy syndrome clinically characterized by focal seizures with prominent auditory symptoms. ADLTE is genetically heterogeneous, and mutations in LGI1 account for fewer than 50% of affected families. Here, we report the identification of causal mutations in reelin (RELN) in seven ADLTE-affected families without LGI1 mutations. We initially investigated 13 ADLTE-affected families by performing SNP-array linkage analysis and whole-exome sequencing and identified three heterozygous missense mutations co-segregating with the syndrome. Subsequent analysis of 15 small ADLTE-affected families revealed four additional missense mutations. 3D modeling predicted that all mutations have structural effects on protein-domain folding. Overall, RELN mutations occurred in 7/40 (17.5%) ADLTE-affected families. RELN encodes a secreted protein, Reelin, which has important functions in both the developing and adult brain and is also found in the blood serum. We show that ADLTE-related mutations significantly decrease serum levels of Reelin, suggesting an inhibitory effect of mutations on protein secretion. We also show that Reelin and LGI1 co-localize in a subset of rat brain neurons, supporting an involvement of both proteins in a common molecular pathway underlying ADLTE. Homozygous RELN mutations are known to cause lissencephaly with cerebellar hypoplasia. Our findings extend the spectrum of neurological disorders associated with RELN mutations and establish a link between RELN and LGI1, which play key regulatory roles in both the developing and adult brain.

RNAIII of the Staphylococcus aureus agr system activates global regulator MgrA by stabilizing mRNA.

RNAIII, the effector of the agr quorum-sensing system, plays a key role in virulence gene regulation in Staphylococcus aureus, but how RNAIII transcriptionally regulates its downstream genes is not completely understood. Here, we show that RNAIII stabilizes mgrA mRNA, thereby increasing the production of MgrA, a global transcriptional regulator that affects the expression of many genes. The mgrA gene is transcribed from two promoters, P1 and P2, to produce two mRNA transcripts with long 5' UTR. Two adjacent regions of the mgrA mRNA UTR transcribed from the upstream P2 promoter, but not the P1 promoter, form a stable complex with two regions of RNAIII near the 5' and 3' ends. We further demonstrate that the interaction has several biological effects. We propose that MgrA can serve as an intermediary regulator through which agr exerts its regulatory function.

Angelman syndrome imprinting center encodes a transcriptional promoter.

Clusters of imprinted genes are often controlled by an imprinting center that is necessary for allele-specific gene expression and to reprogram parent-of-origin information between generations. An imprinted domain at 15q11-q13 is responsible for both Angelman syndrome (AS) and Prader-Willi syndrome (PWS), two clinically distinct neurodevelopmental disorders. Angelman syndrome arises from the lack of maternal contribution from the locus, whereas Prader-Willi syndrome results from the absence of paternally expressed genes. In some rare cases of PWS and AS, small deletions may lead to incorrect parent-of-origin allele identity. DNA sequences common to these deletions define a bipartite imprinting center for the AS-PWS locus. The PWS-smallest region of deletion overlap (SRO) element of the imprinting center activates expression of genes from the paternal allele. The AS-SRO element generates maternal allele identity by epigenetically inactivating the PWS-SRO in oocytes so that paternal genes are silenced on the future maternal allele. Here we have investigated functional activities of the AS-SRO, the element necessary for maternal allele identity. We find that, in humans, the AS-SRO is an oocyte-specific promoter that generates transcripts that transit the PWS-SRO. Similar upstream promoters were detected in bovine oocytes. This result is consistent with a model in which imprinting centers become DNA methylated and acquire maternal allele identity in oocytes in response to transiting transcription.

Mutations in CKAP2L, the human homolog of the mouse Radmis gene, cause Filippi syndrome.

Filippi syndrome is a rare, presumably autosomal-recessive disorder characterized by microcephaly, pre- and postnatal growth failure, syndactyly, and distinctive facial features, including a broad nasal bridge and underdeveloped alae nasi. Some affected individuals have intellectual disability, seizures, undescended testicles in males, and teeth and hair abnormalities. We performed homozygosity mapping and whole-exome sequencing in a Sardinian family with two affected children and identified a homozygous frameshift mutation, c.571dupA (p.Ile191Asnfs(∗)6), in CKAP2L, encoding the protein cytoskeleton-associated protein 2-like (CKAP2L). The function of this protein was unknown until it was rediscovered in mice as Radmis (radial fiber and mitotic spindle) and shown to play a pivotal role in cell division of neural progenitors. Sanger sequencing of CKAP2L in a further eight unrelated individuals with clinical features consistent with Filippi syndrome revealed biallelic mutations in four subjects. In contrast to wild-type lymphoblastoid cell lines (LCLs), dividing LCLs established from the individuals homozygous for the c.571dupA mutation did not show CKAP2L at the spindle poles. Furthermore, in cells from the affected individuals, we observed an increase in the number of disorganized spindle microtubules owing to multipolar configurations and defects in chromosome segregation. The observed cellular phenotypes are in keeping with data from in vitro and in vivo knockdown studies performed in human cells and mice, respectively. Our findings show that loss-of-function mutations in CKAP2L are a major cause of Filippi syndrome.

Molecular cloning, mRNA expression and biological activity of the pheromone biosynthesis activating neuropeptide (PBAN) from the European corn borer, Ostrinia nubilalis.

Pheromone biosynthesis activating neuropeptide (PBAN) is a member of the pyrokinin (FXPRLamide) insect neuropeptides. Here, we report the cloning of the gene Ostnu-PBAN from the E and Z pheromone strains of the European corn borer (ECB), Ostrinia nubilalis (Lepidoptera: Crambidae), a major pest of maize. The Ostnu-PBAN genomic sequence is > 5 kb in length and consists of six exons. The deduced amino acid sequence revealed a 200-residue precursor protein including a signal peptide, a 24-amino acid (aa) diapause hormone, a 37-aa PBAN and three other FXPRLamide neuropeptides. Our in vivo assays suggest that the 37-aa synthetic Ostnu-PBAN is hormonally active in the pheromone gland. It restores sex pheromone production to normal levels in mated females and decapitated virgins of both E and Z cultures. The results of a real-time PCR analysis indicated that Ostnu-PBAN mRNA levels reached a plateau in the brain-suboesophageal ganglion complexes 1 day after eclosion, and mating did not affect the mRNA expression. Three size classes of Ostnu-PBAN mRNA (1.9, 2.0 and 2.1 kb) were obtained, differing only in the length of the 3' untranslated region. However, there was no correlation between sequence divergence and the pheromone composition, voltinism or geographical origin (Hungary, Slovenia, Sweden, Turkey) of ECB moths.

Expanding the view of Clock and cycle gene evolution in Diptera.

We expanded the view of Clock (Clk) and cycle (cyc) gene evolution in Diptera by studying the fruit fly Anastrepha fraterculus (Afra), a Brachycera. Despite the high conservation of clock genes amongst insect groups, striking structural and functional differences of some clocks have appeared throughout evolution. Clk and cyc nucleotide sequences and corresponding proteins were characterized, along with their mRNA expression data, to provide an evolutionary overview in the two major groups of Diptera: Lower Diptera and Higher Brachycera. We found that AfraCYC lacks the BMAL (Brain and muscle ARNT-like) C-terminus region (BCTR) domain and is constitutively expressed, suggesting that AfraCLK has the main transactivation function, which is corroborated by the presence of poly-Q repeats and an oscillatory pattern. Our analysis suggests that the loss of BCTR in CYC is not exclusive of drosophilids, as it also occurs in other Acalyptratae flies such as tephritids and drosophilids, however, but it is also present in some Calyptratae, such as Muscidae, Calliphoridae and Sarcophagidae. This indicates that BCTR is missing from CYC of all higher-level Brachycera and that it was lost during the evolution of Lower Brachycera. Thus, we can infer that CLK protein may play the main role in the CLK\CYC transcription complex in these flies, like in its Drosophila orthologues.

Dissociation of the trimeric gp41 ectodomain at the lipid-water interface suggests an active role in HIV-1 Env-mediated membrane fusion.

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. The actual fusion process involves a switch from a homotrimeric prehairpin intermediate conformation, consisting of parallel coiled-coil helices, to a postfusion state where the ectodomains are arranged as a trimer of helical hairpins, adopting a six-helix bundle (6HB) state. Here, we show by solution NMR spectroscopy that a water-soluble 6HB gp41 ectodomain binds to zwitterionic detergents that contain phosphocholine or phosphatidylcholine head groups and phospholipid vesicles that mimic T-cell membrane composition. Binding results in the dissociation of the 6HB and the formation of a monomeric state, where its two α-helices, N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR), become embedded in the lipid-water interface of the virus and host cell. The atomic structure of the gp41 ectodomain monomer, based on NOE distance restraints and residual dipolar couplings, shows that the NHR and CHR helices remain mostly intact, but they completely lose interhelical contacts. The high affinity of the ectodomain helices for phospholipid surfaces suggests that unzippering of the prehairpin intermediate leads to a state where the NHR and CHR helices become embedded in the host cell and viral membranes, respectively, thereby providing a physical force for bringing these membranes into close juxtaposition before actual fusion.

Clinical utility of next-generation sequencing in the diagnosis of hereditary haemolytic anaemias.

Hereditary haemolytic anaemias are genetically and phenotypically heterogeneous disorders characterized by increased red cell destruction, with consequences ranging from innocuous to severe life-threatening anaemia. Diagnostic laboratories endeavour to assist clinicians reach the exact patient diagnosis by using tests principally based on morphological and biochemical techniques. However, these routine studies may be inconclusive, particularly in newborn infants and when transfusions have recently been administered. Large numbers and size of the potentially involved genes also impose a practical challenge for molecular diagnosis using routine sequencing approaches. To overcome these diagnostic shortcomings, we have utilized next-generation sequencing to provide a high-throughput, highly sensitive assay. We developed a panel interrogating 28 genes encoding cytoskeletal proteins and enzymes with sequencing coverage of the coding regions, splice site junctions, deep intronic and regulatory regions. We then evaluated 19 samples, including infants with unexplained extreme hyperbilirubinaemia and patients with transfusion-dependent haemolytic anaemia. Where possible, inheritance patterns of pathogenic mutations were determined by sequencing of immediate relatives. We conclude that this next-generation sequencing panel could be a cost-effective approach to molecular diagnosis of hereditary haemolytic anaemia, especially when the family history is uninformative or when routine laboratory testing fails to identify the causative haemolytic process.

Role of SRC-like adaptor protein (SLAP) in immune and malignant cell signaling.

SRC-like adaptor protein (SLAP) is an adaptor protein structurally similar to the SRC family protein kinases. Like SRC, SLAP contains an SH3 domain followed by an SH2 domain but the kinase domain has been replaced by a unique C-terminal region. SLAP is expressed in a variety of cell types. Current studies suggest that it regulates signaling of various cell surface receptors including the B cell receptor, the T cell receptor, cytokine receptors and receptor tyrosine kinases which are important regulator of immune and cancer cell signaling. SLAP targets receptors, or its associated components, by recruiting the ubiquitin machinery and thereby destabilizing signaling. SLAP directs receptors to ubiquitination-mediated degradation and controls receptors turnover as well as signaling. Thus, SLAP appears to be an important component in regulating signal transduction required for immune and malignant cells.

Differential increases of specific FMR1 mRNA isoforms in premutation carriers.

BACKGROUND: Over 40% of male and ∼16% of female carriers of a premutation FMR1 allele (55-200 CGG repeats) will develop fragile X-associated tremor/ataxia syndrome, an adult onset neurodegenerative disorder, while about 20% of female carriers will develop fragile X-associated primary ovarian insufficiency. Marked elevation in FMR1 mRNA transcript levels has been observed with premutation alleles, and RNA toxicity due to increased mRNA levels is the leading molecular mechanism proposed for these disorders. However, although the FMR1 gene undergoes alternative splicing, it is unknown whether all or only some of the isoforms are overexpressed in premutation carriers and which isoforms may contribute to the premutation pathology. METHODS: To address this question, we have applied a long-read sequencing approach using single-molecule real-time (SMRT) sequencing and qRT-PCR. RESULTS: Our SMRT sequencing analysis performed on peripheral blood mononuclear cells, fibroblasts and brain tissue samples derived from premutation carriers and controls revealed the existence of 16 isoforms of 24 predicted variants. Although the relative abundance of all mRNA isoforms was significantly increased in the premutation group, as expected based on the bulk increase in mRNA levels, there was a disproportionate (fourfold to sixfold) increase, relative to the overall increase in mRNA, in the abundance of isoforms spliced at both exons 12 and 14, specifically Iso10 and Iso10b, containing the complete exon 15 and differing only in splicing in exon 17. CONCLUSIONS: These findings suggest that RNA toxicity may arise from a relative increase of all FMR1 mRNA isoforms. Interestingly, the Iso10 and Iso10b mRNA isoforms, lacking the C-terminal functional sites for fragile X mental retardation protein function, are the most increased in premutation carriers relative to normal, suggesting a functional relevance in the pathology of FMR1-associated disorders.

Large-scale analysis of tandem repeat variability in the human genome.

Tandem repeats are short DNA sequences that are repeated head-to-tail with a propensity to be variable. They constitute a significant proportion of the human genome, also occurring within coding and regulatory regions. Variation in these repeats can alter the function and/or expression of genes allowing organisms to swiftly adapt to novel environments. Importantly, some repeat expansions have also been linked to certain neurodegenerative diseases. Therefore, accurate sequencing of tandem repeats could contribute to our understanding of common phenotypic variability and might uncover missing genetic factors in idiopathic clinical conditions. However, despite long-standing evidence for the functional role of repeats, they are largely ignored because of technical limitations in sequencing, mapping and typing. Here, we report on a novel capture technique and data filtering protocol that allowed simultaneous sequencing of thousands of tandem repeats in the human genomes of a three generation family using GS-FLX-plus Titanium technology. Our results demonstrated that up to 7.6% of tandem repeats in this family (4% in coding sequences) differ from the reference sequence, and identified a de novo variation in the family tree. The method opens new routes to look at this underappreciated type of genetic variability, including the identification of novel disease-related repeats.