Amplicon-based semiconductor sequencing of human exomes: performance evaluation and optimization strategies.

The Ion Proton platform allows to perform whole exome sequencing (WES) at low cost, providing rapid turnaround time and great flexibility. Products for WES on Ion Proton system include the AmpliSeq Exome kit and the recently introduced HiQ sequencing chemistry. Here, we used gold standard variants from GIAB consortium to assess the performances in variants identification, characterize the erroneous calls and develop a filtering strategy to reduce false positives. The AmpliSeq Exome kit captures a large fraction of bases (>94 %) in human CDS, ClinVar genes and ACMG genes, but with 2,041 (7 %), 449 (13 %) and 11 (19 %) genes not fully represented, respectively. Overall, 515 protein coding genes contain hard-to-sequence regions, including 90 genes from ClinVar. Performance in variants detection was maximum at mean coverage >120×, while at 90× and 70× we measured a loss of variants of 3.2 and 4.5 %, respectively. WES using HiQ chemistry showed ~71/97.5 % sensitivity, ~37/2 % FDR and ~0.66/0.98 F1 score for indels and SNPs, respectively. The proposed low, medium or high-stringency filters reduced the amount of false positives by 10.2, 21.2 and 40.4 % for indels and 21.2, 41.9 and 68.2 % for SNP, respectively. Amplicon-based WES on Ion Proton platform using HiQ chemistry emerged as a competitive approach, with improved accuracy in variants identification. False-positive variants remain an issue for the Ion Torrent technology, but our filtering strategy can be applied to reduce erroneous variants.

Expression pattern of the Kallmann syndrome gene in the olfactory system suggests a role in neuronal targeting.

Kallmann syndrome is a genetic disorder characterized by a defect in olfactory system development, which appears to be due to an abnormality in the migration of olfactory axons and gonadotropin releasing hormone (Gn-RH) producing neurons. The X-linked Kallmann syndrome gene shares significant similarities with molecules involved in neural development. We have now isolated the evolutionarily conserved chicken homologue of the Kallmann gene. In the developing and adult chicken, high levels of expression were found in the mitral cells of the olfactory bulb (the target of olfactory axons) and in the Purkinje cells of the cerebellar cortex, both areas affected in patients with Kallmann syndrome. We propose a model in which the Kallmann syndrome gene product is a signal molecule required for neuronal targeting throughout life.

Different chromosomal localization of the Clcn4 gene in Mus spretus and C57BL/6J mice.

We report the unprecedented finding of a gene with a different map position in two mouse strains. The Clcn4 gene was found to map to the X chromosome in the wild Mediterrean mouse, Mus spretus but to chromosome 7 in the inbred strain of laboratory mouse C57BL/6J. These data indicate that a recent evolutionary rearrangement occurred on the mouse sex chromosomes, very close to the pseudoautosomal region. Our data provide molecular evidence for a major divergence near the pseudoautosomal region, consistent with the hypothesis that hybrid sterility in these species results from abnormal pairing of sex chromosomes during male meiosis.

Identification and mapping of human cDNAs homologous to Drosophila mutant genes through EST database searching.

Cross-species comparison is an effective tool used to identify genes and study their function in both normal and pathological conditions. We have applied the power of Drosophila genetics to the vast resource of human cDNAs represented in the expressed sequence tag (EST) database (dbEST) to identify novel human genes of high biological interest. Sixty-six human cDNAs showing significant homology to genes causing Drosophila mutant phenotypes were identified by screening dbEST using the "text string' option, and their map position was determined using both fluorescence in situ hybridization (FISH) and radiation hybrid mapping. Comparison between these genes and their putative partners in Drosophila may provide important insights into their function in mammals. Furthermore, integration of these genes into the transcription map of the human genome contributes to the positional candidate approach for disease gene identification.

Kallmann syndrome gene on the X and Y chromosomes: implications for evolutionary divergence of human sex chromosomes.

The recently identified gene for X-linked Kallmann syndrome (hypogonadotropic hypogonadism and anosmia) has a closely related homologue on the Y chromosome. The X and Y copies of this gene are located in a large region of X/Y homology, on Xp22.3 and Yq11.2, respectively. Comparison of the structure of the X-linked Kallmann syndrome gene and its Y homologue shed light on the evolutionary history of this region of the human sex chromosomes. Our data show that the Y homologue is not functional. Comparative analysis of X/Y sequence identity at several loci on Xp22.3 and Yq11.2 suggests that the homology between these two regions is the result of a complex series of events which occurred in the recent evolution of sex chromosomes.

SLC7A7, encoding a putative permease-related protein, is mutated in patients with lysinuric protein intolerance.

Lysinuric protein intolerance (LPI, MIM 222700) is an autosomal recessive multisystem disorder found mainly in Finland and Italy. On a normal diet, LPI patients present poor feeding, vomiting, diarrhoea, episodes of hyperammoniaemic coma and failure to thrive. Hepatosplenomegaly, osteoporosis and a life-threatening pulmonary involvement (alveolar proteinosis) are also seen. LPI is caused by defective cationic amino acid (CAA) transport at the basolateral membrane of epithelial cells in kidney and intestine. Metabolic derangement is characterized by increased renal excretion of CAA, reduced CAA absorption from intestine and orotic aciduria. The gene causing LPI was assigned using linkage analysis to chromosome 14q11.2 near the T-cell receptor alpha/delta chains locus, and a critical region has been defined. We have identified two new transcripts (SLC7A8 and SLC7A7) homologous to amino acid transporters, highly expressed in kidney and mapping in the LPI critical region. Mutational analysis of both transcripts revealed that SLC7A7 (for solute carrier family 7, member 7) is mutated in LPI. In five Italian patients, we found either an insertion or deletion in the coding sequence, which provides evidence of a causative role of SLC7A7 in LPI. Furthermore, we detected a splice acceptor change resulting in a frameshift and premature translation termination in four unrelated Finnish patients. This mutation may represent the founder LPI allele in Finland.

Non-type I cystinuria caused by mutations in SLC7A9, encoding a subunit (bo,+AT) of rBAT.

Cystinuria (MIM 220100) is a common recessive disorder of renal reabsorption of cystine and dibasic amino acids. Mutations in SLC3A1, encoding rBAT, cause cystinuria type I (ref. 1), but not other types of cystinuria (ref. 2). A gene whose mutation causes non-type I cystinuria has been mapped by linkage analysis to 19q12-13.1 (Refs 3,4). We have identified a new transcript, encoding a protein (bo, +AT, for bo,+ amino acid transporter) belonging to a family of light subunits of amino acid transporters, expressed in kidney, liver, small intestine and placenta, and localized its gene (SLC7A9) to the non-type I cystinuria 19q locus. Co-transfection of bo,+AT and rBAT brings the latter to the plasma membrane, and results in the uptake of L-arginine in COS cells. We have found SLC7A9 mutations in Libyan-Jews, North American, Italian and Spanish non-type I cystinuria patients. The Libyan Jewish patients are homozygous for a founder missense mutation (V170M) that abolishes b o,+AT amino-acid uptake activity when co-transfected with rBAT in COS cells. We identified four missense mutations (G105R, A182T, G195R and G295R) and two frameshift (520insT and 596delTG) mutations in other patients. Our data establish that mutations in SLC7A9 cause non-type I cystinuria, and suggest that bo,+AT is the light subunit of rBAT.

Overexpression of Human Mutant PANK2 Proteins Affects Development and Motor Behavior of Zebrafish Embryos.

Pantothenate Kinase-Associated Neurodegeneration (PKAN) is a genetic and early-onset neurodegenerative disorder characterized by iron accumulation in the basal ganglia. It is due to mutations in Pantothenate Kinase 2 (PANK2), an enzyme that catalyzes the phosphorylation of vitamin B5, first and essential step in coenzyme A (CoA) biosynthesis. Most likely, an unbalance of the neuronal levels of this important cofactor represents the initial trigger of the neurodegenerative process, yet a complete understanding of the connection between PANK2 malfunctioning and neuronal death is lacking. Most PKAN patients carry mutations in both alleles and a loss of function mechanism is proposed to explain the pathology. When PANK2 mutants were analyzed for stability, dimerization capacity, and enzymatic activity in vitro, many of them showed properties like the wild-type form. To further explore this aspect, we overexpressed the wild-type protein, two mutant forms with reduced kinase activity and two retaining the catalytic activity in zebrafish embryos and analyzed the morpho-functional consequences. While the wild-type protein had no effects, all mutant proteins generated phenotypes that partially resembled those observed in pank2 and coasy morphants and were rescued by CoA and vitamin B5 supplementation. The overexpression of PANK2 mutant forms appears to be associated with perturbation in CoA availability, irrespective of their catalytic activity.

The molecular bases of cystinuria and lysinuric protein intolerance.

Cystinuria and lysinuric protein intolerance are inherited aminoacidurias caused by defective amino-acid transport activities linked to a family of heteromeric amino-acid transporters (HATs). HATs comprise two subunits: co-expression of subunits 4F2hc and y(+)LAT-1 induces the efflux of dibasic amino acids from cells, whereas co-expression of subunits rBAT and b(o,+)AT induces the renal reabsorption and intestinal absorption of cystine and dibasic amino acids at the brush border of epithelial cells. Recently, the role of b(o,+)AT (SLC7A9) in cystinuria (non Type I) and the role of y(+)LAT-1 (SLC7A7) in lysinuric protein intolerance have been demonstrated.

Structure of the mouse growth/differentiation factor 9 gene.

Using a mouse GDF-9 cDNA as a probe, over 20 kb of sequence encompassing the GDF-9 gene was isolated from a mouse 129SvEv genomic library. Sequence analysis of the exons, exon-intron boundaries, and 5'- and 3'-non-translated regions was used to establish the structure of the mouse GDF-9 gene. The GDF-9 gene is encoded by two exons separated by a 2.9 kb intron. Multiple putative transcription start sites are detected between 31 and 57 bp upstream of the start site of translation and a putative polyadenylation signal lies 342 bp 3' of the end of translation. Knowledge of the GDF-9 gene structure will enable us to further understand the role of GDF-9 in ovarian physiology and development.

The expression pattern of a mouse doublesex-related gene is consistent with a role in gonadal differentiation.

The signal for somatic sex determination in mammals, Caenorhabditis elegans and Drosophila melanogaster is chromosomal, but the overall mechanisms do not appear to be conserved between the phyla. However it has been found quite recently that the C. elegans sex-determining gene Mab-3 contains a domain highly homologous to the Drosophila sex-determining gene doublesex (dsx) and shares a similar role. These data suggest that at least some aspects of the regulation of sex determination might be conserved. In humans, a doublesex-related gene (DMRT1) was identified at less than 30 kb from the critical region for sex reversal on chromosome 9p24 (TD9). In order to get insights into the role of DMRT1 in sex determination/differentiation, we have isolated DMRT1 mouse homologue (Dmrt1) and analysed its expression pattern. The gene is expressed in the genital ridges of both sexes during the sex-determining switch and it shows male/female dimorphism at late stages of sex differentiation.

Overexpression of wild-type and mutant mucolipin proteins in mammalian cells: effects on the late endocytic compartment organization.

Mucolipin-1 is a 65-kDa membrane protein encoded by the MCOLN1 gene, which is mutated in patients with mucolipidosis type IV (MLIV), a rare neurodegenerative lysosomal storage disorder. We studied the subcellular localization of wild-type and three different mutant forms (T232P, F408del and F465L) of mucolipin by expressing Myc-tagged proteins in HeLa cells. The overexpressed wild-type mucolipin colocalizes to late endocytic structures and induces an aberrant distribution of these compartments. F408del and F465L MLIV mutant proteins show a distribution similar to the wild-type protein, whereas T232P is retained in the endoplasmic reticulum. Among the mutants, only F408del induces a redistribution of the late endocytic compartment. These findings suggest that the overexpression of the mucolipin cation channel influences the dynamic equilibrium of late endocytic compartments.

Erythrocyte adducin differential properties in the normotensive and hypertensive rats of the Milan strain. Characterization of spleen adducin m-RNA.

Previous studies have shown that erythrocytes from the Milan hypertensive strain of rats (MHS) differ from erythrocytes from the control normotensive strain (MNS). These differences are determined within the stem cells, are genetically associated with the development of hypertension, and are similar to those found between the tubular cells of the two strains. Moreover they seem to be dependent upon the presence of the membrane skeleton proteins. In this paper we describe our studies aimed at identifying some precise protein difference between the membrane skeletons of the two strains, which may cause the cellular differences described above. Milan hypertensive strain and MNS rats were immunized with ghost or membrane skeleton extracts prepared from the other or their own strains. Only MHS rats immunized with MNS ghost or membrane skeleton extracts produced an antibody against a 105 KD protein in about 95% of the animals. This protein has been identified with the recently described cytoskeletal protein adducin on the following bases: the protein binds calmodulin (CaM) and protein kinase C (PKc) in a Ca2+ dependent way. It also binds phosphatidylserine, is the substrate of exogenous PKc, and finally it is purified by high salt extraction of Triton-X100 insoluble erythrocyte cytoskeletons followed by affinity chromatography on CaM-sepharose. Using this antibody the isolation from a mouse spleen library, the characterization and sequencing of a partial cDNA clone coding for this protein has been carried out. In conclusion adducin may be considered a very useful tool to test the hypothesis that the cellular differences between MHS and MNS may be caused by a difference in a membrane skeleton protein.

Detection of p53 mutations by single-strand conformation polymorphisms (SSCP) gel electrophoresis. A comparative study of radioactive and nonradioactive silver-stained SSCP analysis.

p53 mutations are the most common genetic abnormality in humans tumors, but their clinical significance remains to be precisely elucidated. Conventional single-strand conformation polymorphism (SSCP) analysis, a well-established technique for detecting p53 mutations, uses radioactively labeled polymerase chain reaction (PCR) products, which migrate abnormally in the presence of mutations. We performed radioactive PCR-SSCP analysis in a series of 30 formalin-fixed, paraffin-embedded ovarian carcinomas and two cell lines (SW480 and Caov4) harboring known homozygous p53 mutations and compared the results with nonradioactive silver-stained SSCP. The purpose was to assess whether nonradioactive SSCP is suitable for detecting p53 mutations in a rapid, sensitive, cost-effective fashion, without the need of radioactive isotopes. We accomplished PCR amplification of p53 exons 5 through 8 in 26 carcinomas, and radioactive SSCP detected p53 mutations in 13 tumors; three mutations were localized in exon 5, six in exon 6, two in exon 7, and two in exon 8. All mutations were correctly identified with nonradioactive SSCP, except for one exon 8 mutation. To establish the sensitivity of nonradioactive SSCP, DNA samples of SW480 and Caov4 were mixed with increasing amounts (0-90%) of normal DNA and subjected to PCR-SSCP analysis. Mutations were detected until the concentration of SW480 and Caov4 was 15% and 10%, respectively, of the total sample. The results of our investigation demonstrate that nonradioactive silver-stained SSCP is a sensitive, rapid, and simple technique to detect p53 mutations, even in formalin-fixed tissues, and could be easily used to investigate large series of patients to assess the clinical significance of p53 mutations in human tumors.

Detection of beta-nerve growth factor mRNA in the human fetal brain.

Nerve growth factor (NGF) is a trophic molecule recently demonstrated to interact with different structures in the central nervous system. The expression of the beta-NGF mRNA from human fetal cortices at the 15-16th week of gestational age has been demonstrated and quantitated by polymerase chain reaction amplification of the specific cDNA. beta-NGF mRNA expression in the human brain coincides with the period of active differentiation and synaptogenesis, suggesting that the trophic agent plays a role in the cerebral development.

Human fetal brain beta-nerve growth factor cDNA: molecular cloning of 5' and 3' untranslated regions.

The mRNA of beta-nerve growth factor (beta-NGF) has been demonstrated to be present in the human brain, both in adult and fetal stages of development. However, its complete nucleotide sequence is unknown since a full-length cDNA has yet to be isolated. We report here the cloning and complete analysis of the human fetal brain beta-NGF transcript. cDNA synthesis was performed from total brain RNA and overlapping regions of beta-NGF cDNA were enzymatically amplified and sequenced. The central portion of the cDNA was amplified using primers designed on the known genomic DNA sequence. The 5' and 3' unknown regions were amplified by the anchored polymerase chain reaction (PCR) and by complementary DNA Ends-PCR respectively. The latter method is an original variation of Inverted PCR. The sequenced transcript is very similar to the most common form of beta-NGF mRNA present in the mouse central nervous system. In addition, both the 5' and 3' untranslated regions show a high degree of homology to the corresponding murine sequences.

Relationship between apparent optical properties and photosynthetic pigments in the sub-alpine Lake Iseo.

The aim of the present study was the evaluation of methods for estimating the content of bio-physical parameters in lake water on the basis of spectral reflectance measured above water surface, in particular the estimation of chlorophyll-a (chl-a) concentrations. Data sets considered refer to some sampling point located in the sub-alpine meso-eutrophic Lake Iseo, surveyed six times over the period March-July and once in November 1998, as these months were very important for the characterization of the springtime algal bloom, which affect the lake waters yearly. At each point station, limnological observations (chlorophyll, total suspended solids, Secchi disk depth) were conducted simultaneously with optical measurements. The latter consisted of water leaving radiance measured by means of a spectroradiometer above the water surface; moreover, a standard reflector radiance was also measured to obtain the water reflectance. Reflectance spectra were transformed according to two well-documented models and correlated to water quality parameters, to investigate their performances as retrieval algorithms under different conditions and referring to different analytical methods. Results outline the sensitivity of the models to chl-a concentrations, different phytoplankton composition, and the sampling depth.

Sequencing analysis of forty-eight human image cDNA clones similar to Drosophila mutant protein.

We have sequenced 48 human IMAGE cDNA clones selected from the public EST database (dbEST) for their significant homology to Drosophila mutant genes. A dynamically updated analysis report was produced by BlastX and BlastN analysis searches in the latest databases available. This analysis led us to estimate the grade of similarity with homologous genes isolated in other species. Bottlenecks were detected in the sequencing process and here we have presented our problem-solving approach. We think the value of this full-length sequencing project is an enrichment of the sequence database information that is currently available to the human genome community.