Next-Generation Sequencing of Hereditary Hemochromatosis-Related Genes: Novel Likely Pathogenic Variants Found in the Portuguese Population

Hereditary hemochromatosis (HH) is an autosomal recessive disorder characterized by excessive iron absorption resulting in pathologically increased body iron stores. It is typically associated with common HFE gene mutation (p.Cys282Tyr and p.His63Asp). However, in Southern European populations up to one third of HH patients do not carry the risk genotypes. This study aimed to explore the use of next-generation sequencing (NGS) technology to analyse a panel of iron metabolism-related genes (HFE, TFR2, HJV, HAMP, SLC40A1, and FTL) in 87 non-classic HH Portuguese patients. A total of 1241 genetic alterations were detected corresponding to 53 different variants, 13 of which were not described in the available public databases. Among them, five were predicted to be potentially pathogenic: three novel mutations in TFR2 [two missense (p.Leu750Pro and p.Ala777Val) and one intronic splicing mutation (c.967-1G>C)], one missense mutation in HFE (p.Tyr230Cys), and one mutation in the 5'-UTR of HAMP gene (c.-25G>A). The results reported here illustrate the usefulness of NGS for targeted iron metabolism-related gene panels, as a likely cost-effective approach for molecular genetics diagnosis of non-classic HH patients. Simultaneously, it has contributed to the knowledge of the pathophysiology of those rare iron metabolism-related disorders.

Molecular markers for diabetic nephropathy

Type 2 diabetes is one of the most common metabolic disorders in the world. Globally, the prevalence of this disorder is predicted to increase, along with the risk of developing diabetic related complications. One of those complications is diabetic nephropathy, defined by a progressive increase in proteinuria and a gradual decline in renal function. Approximately 25% to 30% of type 2 diabetic individuals develop this complication. However, its underlying genetic mechanisms remain unclear. Thus, the aim of this study is to contribute to the discovery of the genetic mechanisms involved in the development and progression of diabetic nephropathy, through the identification of relevant genetic variants in Portuguese type 2 diabetic individuals. The exomes of 36 Portuguese type 2 diabetic individuals were sequenced on the Ion ProtonTM Sequencer. From those individuals, 19 did not present diabetic nephropathy, being included in the control group, while the 17 individuals that presented the diabetic complication formed the case group. A statistical analysis was then performed to identify candidate common genetic variants, as well as genes accumulating rare variants that could be associated with diabetic nephropathy. From the search for common variants in the study population, the statistically significant (p-value ≤ 0.05) variants rs1051303 and rs1131620 in the LTBP4 gene, rs660339 in UCP2, rs2589156 in RPTOR, rs2304483 in the SLC12A3 gene and rs10169718 present in ARPC2, were considered as the most biologically relevant to the pathogenesis of diabetic nephropathy. The variants rs1051303 and rs1131620, as well as the variants rs660339 and rs2589156 were associated with protective effects in the development of the complication, while rs2304483 and rs10169718 were considered risk variants, being present in individuals with diagnosed diabetic nephropathy. In the rare variants approach, the genes with statistical significance (p-value ≤ 0.05) found, the STAB1 gene, accumulating 9 rare variants, and the CUX1 gene, accumulating 2 rare variants, were identified as the most relevant. Both genes were considered protective, with the accumulated rare variants mainly present in the group without the renal complication. The present study provides an initial analysis of the genetic evidence associated with the development and progression of diabetic nephropathy, and the results obtained may contribute to a deeper understanding of the genetic mechanisms associated with this diabetic complication.

Molecular characterization of human T cell leukemia virus type 1 subtypes in a group of infected individuals diagnosed in Portugal and Spain

Over the past decade, Portugal and Spain received large numbers of immigrants from HTLV-1 endemic areas. Our aim was to investigate the diversity of subtypes circulating in these two countries and the introduction of new variants. We performed a molecular analysis of HTLV-1 strains in patients diagnosed since 1998. LTR and env proviral sequences from 26 individuals were analyzed to generate phylogenetic trees along with reference HTLV-1 subtypes from several geographic origins. Epidemiological and clinical data were recorded. Most subjects were immigrants (57.7%) from South America and Africa. All isolates belonged to the cosmopolitan A subtype. Most carried the transcontinental subgroup A, but five subjects carried subgroup D and one carried subgroup C, previously unreported in Europe. HTLV strains showed separate clusters linked to the patients' geographic origin. Although subjects with HTLV-1 infection tend not to be engaged in high-risk practices, silent dissemination of a broad diversity of HTLV-1 viruses may still occur.

Next-generation sequencing of hereditary hemochromatosis-related genes: novel likely pathogenic variants found in the Portuguese population

Hereditary hemochromatosis (HH) is an autosomal recessive disorder characterized by excessive iron absorption resulting in pathologically increased body iron stores. It is typically associated with common HFE gene mutation (p.Cys282Tyr and p.His63Asp). However, in Southern European populations up to one third of HH patients do not carry the risk genotypes. This study aimed to explore the use of next-generation sequencing (NGS) technology to analyse a panel of iron metabolism-related genes (HFE, TFR2, HJV, HAMP, SLC40A1, and FTL) in 87 non-classic HH Portuguese patients. A total of 1241 genetic alterations were detected corresponding to 53 different variants, 13 of which were not described in the available public databases. Among them, five were predicted to be potentially pathogenic: three novel mutations in TFR2 [two missense (p.Leu750Pro and p.Ala777Val) and one intronic splicing mutation (c.967-1GNC)], one missense mutation in HFE (p.Tyr230Cys), and one mutation in the 5′-UTR of HAMP gene(c.-25GNA). The results reported here illustrate the usefulness of NGS for targeted iron metabolism-related gene panels, as a likely cost-effective approach for molecular genetics diagnosis of non-classic HH patients. Simultaneously, it has contributed to the knowledge of the pathophysiology of those rare iron metabolism-related disorders.

Molecular Profiling of Shiga Toxin-Producing Escherichia coli and Enteropathogenic E. coli Strains Isolated from French Coastal Environments

Shiga toxin-producing Escherichia coli (STEC) and enteropathogenic E. coli (EPEC) strains may be responsible for food-borne infections in humans. Twenty-eight STEC and 75 EPEC strains previously isolated from French shellfish-harvesting areas and their watersheds and belonging to 68 distinguishable serotypes were characterized in this study. High-throughput real-time PCR was used to search for the presence of 75 E. coli virulence-associated gene targets, and genes encoding Shiga toxin (stx) and intimin (eae) were subtyped using PCR tests and DNA sequencing, respectively. The results showed a high level of diversity between strains, with 17 unique virulence gene profiles for STEC and 56 for EPEC. Seven STEC and 15 EPEC strains were found to display a large number or a particular combination of genetic markers of virulence and the presence of stx and/or eae variants, suggesting their potential pathogenicity for humans. Among these, an O26:H11 stx1a eae-β1 strain was associated with a large number of virulence-associated genes (n = 47), including genes carried on the locus of enterocyte effacement (LEE) or other pathogenicity islands, such as OI-122, OI-71, OI-43/48, OI-50, OI-57, and the high-pathogenicity island (HPI). One O91:H21 STEC strain containing 4 stx variants (stx1a, stx2a, stx2c, and stx2d) was found to possess genes associated with pathogenicity islands OI-122, OI-43/48, and OI-15. Among EPEC strains harboring a large number of virulence genes (n, 34 to 50), eight belonged to serotype O26:H11, O103:H2, O103:H25, O145:H28, O157:H7, or O153:H2.

Shaping the Dicot Fruit: Molecular and Genomic Approaches to Fruit Development

The fruit is one of the most complex and important structures produced by flowering plants, and understanding the development and maturation process of fruits in different angiosperm species with diverse fruit structures is of immense interest. In the work presented here, molecular genetics and genomic analysis are used to explore the processes that form the fruit in two species: The model organism Arabidopsis and the diploid strawberry Fragaria vesca. One important basic question concerns the molecular genetic basis of fruit patterning. A long-standing model of Arabidopsis fruit (the gynoecium) patterning holds that auxin produced at the apex diffuses downward, forming a gradient that provides apical-basal positional information to specify different tissue types along the gynoecium’s length. The proposed gradient, however, has never been observed and the model appears inconsistent with a number of observations. I present a new, alternative model, wherein auxin acts to establish the adaxial-abaxial domains of the carpel primordia, which then ensures proper development of the final gynoecium. A second project utilizes genomics to identify genes that regulate fruit color by analyzing the genome sequences of Fragaria vesca, a species of wild strawberry. Shared and distinct SNPs among three F. vesca accessions were identified, providing a foundation for locating candidate mutations underlying phenotypic variations among different F. vesca accessions. Through systematic analysis of relevant SNP variants, a candidate SNP in FveMYB10 was identified that may underlie the fruit color in the yellow-fruited accessions, which was subsequently confirmed by functional assays. Our lab has previously generated extensive RNA-sequencing data that depict genome-scale gene expression profiles in F. vesca fruit and flower tissues at different developmental stages. To enhance the accessibility of this dataset, the web-based eFP software was adapted for this dataset, allowing visualization of gene expression in any tissues by user-initiated queries. Together, this thesis work proposes a well-supported new model of fruit patterning in Arabidopsis and provides further resources for F. vesca, including genome-wide variant lists and the ability to visualize gene expression. This work will facilitate future work linking traits of economic importance to specific genes and gaining novel insights into fruit patterning and development.

Caracterização morfológica e molecular de isolados virais obtidos de larvas do complexo Heliothinae.

2016

Novel ciliate genetic code variants including the reassignment of all three stop codons to sense codons in Condylostoma magnum

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum. We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.

Estudio molecular en dos hermanas afectadas por ictiosis congénita autosómica recesiva : descripción de una nueva mutación causal en TGM1

Se describe la variante homocigota c.320-2A>G de TGM1 en dos hermanas con ictiosis congénita autosómica recesiva. El clonaje de los transcritos generados por esta variante permitió identificar tres mecanismos moleculares de splicing alternativos.

No evidence for association of ataxia-telangiectasia mutated gene T2119C and C3161G amino acid substitution variants with risk of breast cancer

Background There is evidence that certain mutations in the double-strand break repair pathway ataxia-telangiectasia mutated gene act in a dominant-negative manner to increase the risk of breast cancer. There are also some reports to suggest that the amino acid substitution variants T2119C Ser707Pro and C3161G Pro1054Arg may be associated with breast cancer risk. We investigate the breast cancer risk associated with these two nonconservative amino acid substitution variants using a large Australian population-based case–control study. Methods The polymorphisms were genotyped in more than 1300 cases and 600 controls using 5' exonuclease assays. Case–control analyses and genotype distributions were compared by logistic regression. Results The 2119C variant was rare, occurring at frequencies of 1.4 and 1.3% in cases and controls, respectively (P = 0.8). There was no difference in genotype distribution between cases and controls (P = 0.8), and the TC genotype was not associated with increased risk of breast cancer (adjusted odds ratio = 1.08, 95% confidence interval = 0.59–1.97, P = 0.8). Similarly, the 3161G variant was no more common in cases than in controls (2.9% versus 2.2%, P = 0.2), there was no difference in genotype distribution between cases and controls (P = 0.1), and the CG genotype was not associated with an increased risk of breast cancer (adjusted odds ratio = 1.30, 95% confidence interval = 0.85–1.98, P = 0.2). This lack of evidence for an association persisted within groups defined by the family history of breast cancer or by age. Conclusion The 2119C and 3161G amino acid substitution variants are not associated with moderate or high risks of breast cancer in Australian women.

Vibrational spectroscopy of phthalocyanine and naphthalocyanine in sandwich-type (na)phthalocyaninato and porphyrinato rare earth complexes : (Part 10) The infrared and Raman characteristics of phthalocyanine in heteroleptic bis(phthalocyaninato) rare earth complexes with decreased molecular symmetry

The infrared (IR) spectroscopic data for a series of eleven heteroleptic bis(phthalocyaninato) rare earth complexes MIII(Pc)[Pc(α-OC5H11)4] (M = Sm–Lu, Y) [H2Pc = unsubstituted phthalocyanine, H2Pc(α-OC5H11)4 = 1,8,15,22-tetrakis(3-pentyloxy)phthalocyanine] have been collected with 2 cm−1 resolution. Raman spectroscopic properties in the range of 500–1800 cm−1 for these double-decker molecules have also been comparatively studied using laser excitation sources emitting at 632.8 and 785 nm. Both the IR and Raman spectra for M(Pc)[Pc(α-OC5H11)4] are more complicated than those of homoleptic bis(phthalocyaninato) rare earth analogues due to the decreased molecular symmetry of these double-decker compounds, namely C4. For this series, the IR Pc√− marker band appears as an intense absorption at 1309–1317 cm−1, attributed to the pyrrole stretching. With laser excitation at 632.8 nm, Raman vibrations derived from isoindole ring and aza stretchings in the range of 1300–1600 cm−1 are selectively intensified. In contrast, when excited with laser radiation of 785 nm, the ring radial vibrations of isoindole moieties and dihedral plane deformations between 500 and 1000 cm−1 for M(Pc)[Pc(α-OC5H11)4] intensify to become the strongest scatterings. Both techniques reveal that the frequencies of pyrrole stretching, isoindole breathing, isoindole stretchings, aza stretchings and coupling of pyrrole and aza stretchings depend on the rare earth ionic size, shifting to higher energy along with the lanthanide contraction due to the increased ring-ring interaction across the series. The assignments of the vibrational bands for these compounds have been made and discussed in relation to other unsubstituted and substituted bis(phthalocyaninato) rare earth analogues, such as M(Pc)2 and M(OOPc)2 [H2OOPc = 2,3,9,10,16,17,23,24-octakis(octyloxy)phthalocyanine].

Time-Dependent Density Functional Molecular Orbital and Excited State Calculations on Bis(porphyrinyl)butadiynes in the Monocationic, Neutral, Monoanionic, and Dianionic Oxidation States

We report a theoretical study of the multiple oxidation states (1+, 0, 1−, and 2−) of a meso,meso-linked diporphyrin, namely bis[10,15,20-triphenylporphyrinatozinc(II)-5-yl]butadiyne (4), using Time-Dependent Density Functional Theory (TDDFT). The origin of electronic transitions of singlet excited states is discussed in comparison to experimental spectra for the corresponding oxidation states of the close analogue bis{10,15,20-tris[3‘,5‘-di-tert-butylphenyl]porphyrinatozinc(II)-5-yl}butadiyne (3). The latter were measured in previous work under in situ spectroelectrochemical conditions. Excitation energies and orbital compositions of the excited states were obtained for these large delocalized aromatic radicals, which are unique examples of organic mixed-valence systems. The radical cations and anions of butadiyne-bridged diporphyrins such as 3 display characteristic electronic absorption bands in the near-IR region, which have been successfully predicted with use of these computational methods. The radicals are clearly of the “fully delocalized” or Class III type. The key spectral features of the neutral and dianionic states were also reproduced, although due to the large size of these molecules, quantitative agreement of energies with observations is not as good in the blue end of the visible region. The TDDFT calculations are largely in accord with a previous empirical model for the spectra, which was based simplistically on one-electron transitions among the eight key frontier orbitals of the C4 (1,4-butadiyne) linked diporphyrins.