Integrating structural variant calling, annotation and prioritization into whole genome analysis workflows: a practical application in the molecular diagnosis of neurodevelopmental disorders

Background: WGS is increasingly used as a first-line diagnostic test for patients with rare genetic diseases such as neurodevelopmental disorders (NDD). Clinical applications require a robust infrastructure to support processing, storage and analysis of WGS data. The identification and interpretation of SVs from WGS data also needs to be improved. Finally, there is a need for a prioritization system that enables downstream clinical analysis and facilitates data interpretation. Here, we present the results of a clinical application of WGS in a cohort of patients with NDD. Methods: We developed highly portable workflows for processing WGS data, including alignment, quality control, and variant calling of SNVs and SVs. A benchmark analysis of state-of-the-art SV detection tools was performed to select the most accurate combination for SV calling. A gene-based prioritization system was also implemented to support variant interpretation. Results: Using a benchmark analysis, we selected the most accurate combination of tools to improve SV detection from WGS data and build a dedicated pipeline. Our workflows were used to process WGS data from 77 NDD patient-parent families. The prioritization system supported downstream analysis and enabled molecular diagnosis in 32% of patients, 25% of which were SVs and suggested a potential diagnosis in 20% of patients, requiring further investigation to achieve diagnostic certainty. Conclusion: Our data suggest that the integration of SNVs and SVs is a main factor that increases diagnostic yield by WGS and show that the adoption of a dedicated pipeline improves the process of variant detection and interpretation.

Discovery Of A Rare Pterosaur Bone Bed In A Cretaceous Desert With Insights On Ontogeny And Behavior Of Flying Reptiles.

A pterosaur bone bed with at least 47 individuals (wing spans: 0.65-2.35 m) of a new species is reported from southern Brazil from an interdunal lake deposit of a Cretaceous desert, shedding new light on several biological aspects of those flying reptiles. The material represents a new pterosaur, Caiuajara dobruskii gen. et sp. nov., that is the southermost occurrence of the edentulous clade Tapejaridae (Tapejarinae, Pterodactyloidea) recovered so far. Caiuajara dobruskii differs from all other members of this clade in several cranial features, including the presence of a ventral sagittal bony expansion projected inside the nasoantorbital fenestra, which is formed by the premaxillae; and features of the lower jaw, like a marked rounded depression in the occlusal concavity of the dentary. Ontogenetic variation of Caiuajara dobruskii is mainly reflected in the size and inclination of the premaxillary crest, changing from small and inclined (∼ 115°) in juveniles to large and steep (∼ 90°) in adults. No particular ontogenetic features are observed in postcranial elements. The available information suggests that this species was gregarious, living in colonies, and most likely precocial, being able to fly at a very young age, which might have been a general trend for at least derived pterosaurs.

Tubular Variant Of Basal Cell Adenoma Shares Immunophenotypical Features With Normal Intercalated Ducts And Is Closely Related To Intercalated Duct Lesions Of Salivary Gland.

The morphological criteria for identification of intercalated duct lesions (IDLs) of salivary glands have been defined recently. It has been hypothesised that IDL could be a precursor of basal cell adenoma (BCA). BCAs show a variety of histological patterns, and the tubular variant is the one that presents the strongest resemblance with IDLs. The aim of this study was to analyse the morphological and immunohistochemical profiles of IDLs and BCAs classified into tubular and non-tubular subtypes, to determine whether or not IDL and tubular BCA represent distinct entities. Eight IDLs, nine tubular BCAs and 19 non-tubular BCAs were studied. All tubular BCAs contained IDL-like areas, which represented 20-70% of the tumour. In non-tubular BCA, IDL-like areas were occasional and small (<5%). One patient presented IDLs, tubular BCAs and IDL/tubular BCA combined lesions. Luminal ductal cells of IDLs and tubular BCAs exhibited positivity for CK7, lysozyme, S100 and DOG1. In the non-tubular BCA group, few luminal cells exhibited such an immunoprofile; they were mainly CK14-positive. Basal/myoepithelial cells of IDLs, tubular BCAs and non-tubular BCAs were positive for CK14, calponin, α-SMA and p63; they were more numerous in BCA lesions. IDL, tubular BCA and non-tubular BCA form a continuum of lesions in which IDLs are related closely to tubular BCA. In both, the immunoprofile of luminal and myoepithelial cells recapitulates the normal intercalated duct. The difference between the adenoma-like subset of IDLs and tubular BCA rests mainly on the larger numbers of myoepithelial cells in the latter. Our findings indicate that at least some BCAs can arise via IDLs.

Frequency Of The Allelic Variant C.1150t > C In Exon 10 Of The Fibroblast Growth Factor Receptor 3 (fgfr3) Gene Is Not Increased In Patients With Pathogenic Mutations And Related Chondrodysplasia Phenotypes.

Mutations in the FGFR3 gene cause the phenotypic spectrum of FGFR3 chondrodysplasias ranging from lethal forms to the milder phenotype seen in hypochondroplasia (Hch). The p.N540K mutation in the FGFR3 gene occurs in ∼70% of individuals with Hch, and nearly 30% of individuals with the Hch phenotype have no mutations in the FGFR3, which suggests genetic heterogeneity. The identification of a severe case of Hch associated with the typical mutation c.1620C > A and the occurrence of a c.1150T > C change that resulted in a p.F384L in exon 10, together with the suspicion that this second change could be a modulator of the phenotype, prompted us to investigate this hypothesis in a cohort of patients. An analysis of 48 patients with FGFR3 chondrodysplasia phenotypes and 330 healthy (control) individuals revealed no significant difference in the frequency of the C allele at the c.1150 position (p = 0.34). One patient carrying the combination `pathogenic mutation plus the allelic variant c.1150T > C' had a typical achondroplasia (Ach) phenotype. In addition, three other patients with atypical phenotypes showed no association with the allelic variant. Together, these results do not support the hypothesis of a modulatory role for the c.1150T > C change in the FGFR3 gene.

Variant Rh Alleles And Rh Immunisation In Patients With Sickle Cell Disease.

Alloimmunisation is a major complication in patients with sickle cell disease (SCD) receiving red blood cell (RBC) transfusions and despite provision of Rh phenotyped RBC units, Rh antibodies still occur. These antibodies in patients positive for the corresponding Rh antigen are considered autoantibodies in many cases but variant RH alleles found in SCD patients can also contribute to Rh alloimmunisation. In this study, we characterised variant RH alleles in 31 SCD patients who made antibodies to Rh antigens despite antigen-positive status and evaluated the clinical significance of the antibodies produced. RHD and RHCE BeadChip™ from BioArray Solutions and/or amplification and sequencing of exons were used to identify the RH variants. The serological features of all Rh antibodies in antigen-positive patients were analysed and the clinical significance of the antibodies was evaluated by retrospective analysis of the haemoglobin (Hb) levels before and after transfusion; the change from baseline pre-transfusion Hb and the percentage of HbS were also determined. We identified variant RH alleles in 31/48 (65%) of SCD patients with Rh antibodies. Molecular analyses revealed the presence of partial RHD alleles and variant RHCE alleles associated with altered C and e antigens. Five patients were compound heterozygotes for RHD and RHCE variants. Retrospective analysis showed that 42% of antibodies produced by the patients with RH variants were involved in delayed haemolytic transfusion reactions or decreased survival of transfused RBC. In this study, we found that Rh antibodies in SCD patients with RH variants can be clinically significant and, therefore, matching patients based on RH variants should be considered.

Rare earth element and yttrium geochemistry applied to the genetic study of cryolite ore at the Pitinga Mine (Amazon, Brazil)

This work aims at the geochemical study of Pitinga cryolite mineralization through REE and Y analyses in disseminated and massive cryolite ore deposits, as well as in fluorite occurrences. REE signatures in fluorite and cryolite are similar to those in the Madeira albite granite. The highest ΣREE values are found in magmatic cryolite (677 to 1345 ppm); ΣREE is lower in massive cryolite. Average values for the different cryolite types are 10.3 ppm, 6.66 ppm and 8.38 ppm (for nucleated, caramel and white types, respectively). Disseminated fluorite displays higher ΣREE values (1708 and 1526ppm) than fluorite in late veins(34.81ppm). Yttrium concentration is higher in disseminated fluorite and in magmatic cryolite. The evolution of several parameters (REEtotal, LREE/HREE, Y) was followed throughout successive stages of evolution in albite granites and associated mineralization. At the end of the process, late cryolite was formed with low REEtotal content. REE data indicate that the MCD was formed by, and the disseminated ore enriched by (additional formation of hydrothermal disseminated cryolite), hydrothermal fluids, residual from albite granite. The presence of tetrads is poorly defined, although nucleated, caramel and white cryolite types show evidence for tetrad effect.

Pharmacogenetics of OATP Transporters Reveals That SLCO1B1 c.388A > G Variant Is Determinant of Increased Atorvastatin Response

Aims: The relationship between variants in SLCO1B1 and SLCO2B1 genes and lipid-lowering response to atorvastatin was investigated. Material and Methods: One-hundred-thirty-six unrelated individuals with hypercholesterolemia were selected and treated with atorvastatin (10 mg/day/4 weeks). They were genotyped with a panel of ancestry informative markers for individual African component of ancestry (ACA) estimation by SNaPshot (R) and SLCO1B1 (c.388A>G, c.463C>A and c.521T>C) and SLCO2B1 (-71T>C) gene polymorphisms were identified by TaqMan (R) Real-time PCR. Results: Subjects carrying SLCO1B1 c.388GG genotype exhibited significantly high low-density lipoprotein (LDL) cholesterol reduction relative to c.388AA+c.388AG carriers (41 vs. 37%, p = 0.034). Haplotype analysis revealed that homozygous of SLCO1B1*15 (c.521C and c.388G) variant had similar response to statin relative to heterozygous and non-carriers. A multivariate logistic regression analysis confirmed that c.388GG genotype was associated with higher LDL cholesterol reduction in the study population (OR: 3.2, CI95%: 1.3-8.0, p < 0.05). Conclusion: SLCO1B1 c.388A>G polymorphism causes significant increase in atorvastatin response and may be an important marker for predicting efficacy of lipid-lowering therapy.

ABO genotyping in leukemia patients reveals new ABO variant alleles

The ABO blood group is the most important blood group system in transfusion medicine and organ transplantation. To date, more than 160 ABO alleles have been identified by molecular investigation. Almost all ABO genotyping studies have been performed in blood donors and families and for investigation of ABO subgroups detected serologically. The aim of the present study was to perform ABO genotyping in patients with leukemia. Blood samples were collected from 108 Brazilian patients with chronic myeloid leukemia (N = 69), chronic lymphoid leukemia (N = 13), acute myeloid leukemia (N = 15), and acute lymphoid leukemia (N = 11). ABO genotyping was carried out using allele specific primer polymerase chain reaction followed by DNA sequencing. ABO*001 was the most common allele found, followed by ABO*022 and by ABO*A103. We identified 22 new ABO*(variants) in the coding region of the ABO gene in 25 individuals with leukemia (23.2%). The majority of ABO variants was detected in O alleles (15/60.0%). In 5 of 51 samples typed as blood group O (9.8%), we found non-deletional ABO*O alleles. Elucidation of the diversity of this gene in leukemia and in other diseases is important for the determination of the effect of changes in an amino acid residue on the specificity and activity of ABO glycosyltransferases and their function. In conclusion, this is the first report of a large number of patients with leukemia genotyped for ABO. The findings of this study indicate that there is a high level of recombinant activity in the ABO gene in leukemia patients, revealing new ABO variants.

The International LAM Registry: A Component of an Innovative Web-Based Clinician, Researcher, and Patient-Driven Rare Disease Research Platform

Background: A relative friability to capture a sufficiently large patient population in any one geographic location has traditionally limited research into rare diseases. Methods and Results: Clinicians interested in the rare disease lymphangioleiomyomatosis (LAM) have worked with the LAM Treatment Alliance, the MIT Media Lab, and Clozure Associates to cooperate in the design of a state-of-the-art data coordination platform that can be used for clinical trials and other research focused on the global LAM patient population. This platform is a component of a set of web-based resources, including a patient self-report data portal, aimed at accelerating research in rare diseases in a rigorous fashion. Conclusions: Collaboration between clinicians, researchers, advocacy groups, and patients can create essential community resource infrastructure to accelerate rare disease research. The International LAM Registry is an example of such an effort.

TCF7L2 variant genotypes and type 2 diabetes risk in Brazil: significant association, but not a significant tool for risk stratification in the general population

Background: Genetic polymorphisms of the TCF7L2 gene are strongly associated with large increments in type 2 diabetes risk in different populations worldwide. In this study, we aimed to confirm the effect of the TCF7L2 polymorphism rs7903146 on diabetes risk in a Brazilian population and to assess the use of this genetic marker in improving diabetes risk prediction in the general population. Methods: We genotyped the single nucleotide polymorphisms (SNP) rs7903146 of the TCF7L2 gene in 560 patients with known coronary disease enrolled in the MASS II (Medicine, Angioplasty, or Surgery Study) Trial and in 1,449 residents of Vitoria, in Southeast Brazil. The associations of this gene variant to diabetes risk and metabolic characteristics in these two different populations were analyzed. To access the potential benefit of using this marker for diabetes risk prediction in the general population we analyzed the impact of this genetic variant on a validated diabetes risk prediction tool based on clinical characteristics developed for the Brazilian general population. Results: SNP rs7903146 of the TCF7L2 gene was significantly associated with type 2 diabetes in the MASS-II population (OR = 1.57 per T allele, p = 0.0032), confirming, in the Brazilian population, previous reports of the literature. Addition of this polymorphism to an established clinical risk prediction score did not increased model accuracy (both area under ROC curve equal to 0.776). Conclusion: TCF7L2 rs7903146 T allele is associated with a 1.57 increased risk for type 2 diabetes in a Brazilian cohort of patients with known coronary heart disease. However, the inclusion of this polymorphism in a risk prediction tool developed for the general population resulted in no improvement of performance. This is the first study, to our knowledge, that has confirmed this recent association in a South American population and adds to the great consistency of this finding in studies around the world. Finally, confirming the biological association of a genetic marker does not guarantee improvement on already established screening tools based solely on demographic variables.

EVOLUTION OF DISPARITY BETWEEN THE REGULAR AND VARIANT PHLOEM IN BIGNONIEAE (BIGNONIACEAE)

Premise of the study: The phloem is a plant tissue with a critical role in plant nutrition and signaling. However, little is still known about the evolution of this tissue. In lianas of the Bignoniaceae, two distinct types of phloem coexist: a regular and a variant phloem. The cells associated with these two phloem types are known to be anatomically different; however, it is still unclear what steps were involved in the evolution of such differences. Methods: Here we studied the anatomical development of the regular and variant phloem in representatives of all 21 genera of Bignonieae and used a phylogenetic framework to investigate the timing of changes associated with the evolution of each phloem type. Key results: We found that the variant phloem always appears in a determinate location, between the leaf orthostichies. Furthermore, the variant phloem was mostly occupied by very wide sieve tubes and generally included a higher concentration of fibers, indicating an increase in conduction and mechanical support. On the other hand, the regular phloem included much more parenchyma, more and wider rays, and tiny sieve tubes that resembled terminal sieve tubes from plants with seasonal formation of vascular tissues; these findings suggest reduced conduction and higher storage capacity in the regular phloem. Conclusions: Overall, differences between the regular and variant phloem increased over time, leading to further specialization in conduction in the variant phloem and an increase in storage specialization in the regular phloem.

First discovery of a rare polygyne colony in the stingless bee Melipona quadrifasciata (Apidae, Meliponini)

FAPESP

Search for the rare decay K(L)->pi(0)pi(0)gamma

The KTeV E799 experiment has conducted a search for the rare decay K(L)->pi(0)pi(0)gamma via the topology K(L)->pi(0)pi(0)(D)gamma (where pi(0)(D)->gamma e(+)e(-)). Because of Bose statistics of the pi(0) pair and the real nature of the photon, the K(L)->pi(0)pi(0)gamma decay is restricted to proceed at lowest order by the CP conserving direct emission (DE) of an E2 electric quadrupole photon. The rate of this decay is interesting theoretically since chiral perturbation theory predicts that this process vanishes at level O(p(4)). Therefore, this mode probes chiral perturbation theory at O(p(6)). In this paper we report a determination of an upper limit of 2.43x10(-7) (90% CL) for K(L)->pi(0)pi(0)gamma. This is approximately a factor of 20 lower than previous results.

Structural and optical properties of rare earth-doped (Ba(0.77)Ca(0.23))(1-x)(Sm, Nd, Pr, Yb)(x)TiO(3)

The structural, dielectric, and vibrational properties of pure and rare earth (RE)-doped Ba(0.77) Ca(0.23)TiO(3) (BCT23; RE = Nd, Sm, Pr, Yb) ceramics obtained via solid-state reaction were investigated. The pure and RE-doped BCT23 ceramics sintered at 1450 degrees C in air for 4 h showed a dense microstructure in all ceramics. The use of RE ions as dopants introduced lattice-parameter changes that manifested in the reduction of the volume of the unit cell. RE-doped BCT23 samples exhibit a more homogenous microstructure due to the absence of a Ti-rich phase in the grain boundaries as demonstrated by scanning electron microscopy imaging. The incorporation of REs led to perturbations of the local symmetry of TiO(6) octahedra and the creation of a new Raman mode. The results of Raman scattering measurements indicated that the Curie temperature of the ferroelectric phase transition depends on the RE ion and ion content, with the Curie temperature shifting toward lower values as the RE content increases, with the exception of Yb(3+) doping, which did not affect the ferroelectric phase transition temperature. The phase transition behavior is explained using the standard soft mode model. Electronic paramagnetic resonance measurements showed the existence of Ti vacancies in the structure of RE-doped BCT23. Defects are created via charge compensation mechanisms due to the incorporation of elements with a different valence state relative to the ions of the pure BCT23 host. It is concluded that the Ti vacancies are responsible for the activation of the Raman mode at 840 cm(-1), which is in agreement with lattice dynamics calculations. (c) 2011 American Institute of Physics. [doi:10.1063/1.3594710]

Hitting and returning to rare events for all alpha-mixing processes

We prove that for any a-mixing stationary process the hitting time of any n-string A(n) converges, when suitably normalized, to an exponential law. We identify the normalization constant lambda(A(n)). A similar statement holds also for the return time. To establish this result we prove two other results of independent interest. First, we show a relation between the rescaled hitting time and the rescaled return time, generalizing a theorem of Haydn, Lacroix and Vaienti. Second, we show that for positive entropy systems, the probability of observing any n-string in n consecutive observations goes to zero as n goes to infinity. (c) 2010 Elsevier B.V. All rights reserved.