Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

TUMOR PROGRESSION, METASTATIC POTENTIAL AND GENOMIC INSTABILITY

To investigate the hypothesis that increased malignant potential correlates with increased levels of genetic instability, the following parameters of instability were measured: (1) spontaneous mutation rates for ouabain resistance in murine cell lines of different malignant potentials, (2) the background prevalence of 6-thioguanine (6-TG) resistance in clone 4 (highly metastatic) and clone 19 (poorly metastatic) of the K1735 murine melanoma, (3) the prevalence of ouabain resistant variants in three murine cell lines and their variants after exposure to the mutagen MNNG, (4) the rate of generation of major karyotypic abnormalities in B16 F1 (poorly metastatic) and B16 F10 (highly metastatic) murine melanoma, and (5) analysis of the G-banded karyotypes of cloned B16 F1 and B16 F10 melanoma.^ No correlation of increased spontaneous mutation rates with increased malignant potential was found in repeated experiments with three murine cell lines and their variants of different malignant potential. The background prevalence of g-TG resistance was not significantly different for the poorly and highly metastatic clones of K1735 melanoma. The studies with MNNG-induced mutation showed no increased sensitivity of the highly metastatic variants of the three murine cell lines to mutagenesis. Neither did the rate of generation of major karyotypic abnormalities correlate with malignant potential. However, certain karyotypic differences were demonstrated after G-banding of the B16 F1 and F10 melanomas.^ One hypothesis which is consistent with these results is that the rate of generation of genetic abnormalities need not be strongly related to the degree of malignant potential. An increased prevalence of genetic changes may merely reflect the accumulation of abnormalities while their rate of production remains constant. The presence of specific nonrandom changes likely is the main determinant of malignant potential rather than the rate of production of random changes. ^

The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Impact of range expansions on current human genomic diversity

The patterns of population genetic diversity depend to a large extent on past demographic history. Most human populations are known to have gone recently through a series of range expansions within and out of Africa, but these spatial expansions are rarely taken into account when interpreting observed genomic diversity, possibly because they are difficult to model. Here we review available evidence in favour of range expansions out of Africa, and we discuss several of their consequences on neutral and selected diversity, including some recent observations on an excess of rare neutral and selected variants in large samples. We further show that in spatially subdivided populations, the sampling strategy can severely impact the resulting genetic diversity and be confounded by past demography. We conclude that ignoring the spatial structure of human population can lead to some misinterpretations of extant genetic diversity.

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

Privacy-preserving genomic testing in the clinic: a model using HIV treatment

PURPOSE The implementation of genomic-based medicine is hindered by unresolved questions regarding data privacy and delivery of interpreted results to health-care practitioners. We used DNA-based prediction of HIV-related outcomes as a model to explore critical issues in clinical genomics. METHODS We genotyped 4,149 markers in HIV-positive individuals. Variants allowed for prediction of 17 traits relevant to HIV medical care, inference of patient ancestry, and imputation of human leukocyte antigen (HLA) types. Genetic data were processed under a privacy-preserving framework using homomorphic encryption, and clinical reports describing potentially actionable results were delivered to health-care providers. RESULTS A total of 230 patients were included in the study. We demonstrated the feasibility of encrypting a large number of genetic markers, inferring patient ancestry, computing monogenic and polygenic trait risks, and reporting results under privacy-preserving conditions. The average execution time of a multimarker test on encrypted data was 865 ms on a standard computer. The proportion of tests returning potentially actionable genetic results ranged from 0 to 54%. CONCLUSIONS The model of implementation presented herein informs on strategies to deliver genomic test results for clinical care. Data encryption to ensure privacy helps to build patient trust, a key requirement on the road to genomic-based medicine.Genet Med advance online publication 14 January 2016Genetics in Medicine (2016); doi:10.1038/gim.2015.167.

Identification of potentially pathogenic variants in candidate breast and ovarian cancer susceptibility genes in BRCAX patients

Mutations within the BRCA1 and BRCA2 genes account for approximately 20% of hereditary breast cancers, with a further 10%–15% being attributable to rare mutations in moderate-risk genes and common variants in low-risk genes. The genes harbouring mutations in the remaining ∼65% of hereditary breast cancers are unknown. The identification of mutation carriers in hereditary breast and ovarian cancer (hboc) families is critical for determining who is most at risk of developing the disease and therefore who should be offered risk-reducing procedures or more intensive screening, or both.

Many of the high- and moderate-risk genes for hereditary breast cancers encode proteins that work in concert to maintain genomic stability and in dna damage signalling and repair. A novel BRCA1 protein complex identified within the research group whose target genes are involved in dna repair provided novel candidates for hboc susceptibility genes. These 12 candidate genes were sequenced in a cohort of 675 affected individuals from the Kathleen Cunningham Foundation Consortium for Research into Familial Breast Cancer (kConFab) with hereditary breast or ovarian cancer, but with no mutations in known susceptibility genes (BRCAx patients). This analysis identified 20 individuals (each from a different BRCAx family) with different potentially pathogenic variants across 6 of the candidate hboc susceptibility genes. The family members of each BRCAx index case were tested for the presence of the specific mutation identified in the proband to examine segregation with disease. To further expand on the potential role of the novel candidate hboc susceptibility genes identified in this study, the genetic variation of a second cohort of 520 Northern Irish BRCAx patients is being characterized using a 61-gene panel.

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.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

Background: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Contribution of rare damaging variants in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a heterogeneous and highly heritable neurodevelopmental disorder with a complex genetic architecture, consisting of a combination of common low-risk and more penetrant rare variants. This PhD project aimed to explore the contribution of rare variants in ASD susceptibility through NGS approaches in a cohort of 106 ASD families including 125 ASD individuals. Firstly, I explored the contribution of inherited rare variants towards the ASD phenotype in a girl with a maternally inherited pathogenic NRXN1 deletion. Whole exome sequencing of the trio family identified an increased burden of deleterious variants in the proband that could modulate the CNV penetrance and determine the disease development. In the second part of the project, I investigated the role of rare variants emerging from whole genome sequencing in ASD aetiology. To properly manage and analyse sequencing data, a robust and efficient variant filtering and prioritization pipeline was developed, and by its application a stringent set of rare recessive-acting and ultra-rare variants was obtained. As a first follow-up, I performed a preliminary analysis on de novo variants, identifying the most likely deleterious variants and highlighting candidate genes for further analyses. In the third part of the project, considering the well-established involvement of calcium signalling in the molecular bases of ASD, I investigated the role of rare variants in voltage-gated calcium channels genes, that mainly regulate intracellular calcium concentration, and whose alterations have been correlated with enhanced ASD risk. Specifically, I functionally tested the effect of rare damaging variants identified in CACNA1H, showing that CACNA1H variation may be involved in ASD development by additively combining with other high risk variants. This project highlights the challenges in the analysis and interpretation of variants from NGS analysis in ASD, and underlines the importance of a comprehensive assessment of the genomic landscape of ASD individuals.

Genomic And Chemical Insights Into Biosurfactant Production By The Mangrove-derived Strain Bacillus Safensis Ccma-560.

Many Bacillus species can produce biosurfactant, although most of the studies on lipopeptide production by this genus have been focused on Bacillus subtilis. Surfactants are broadly used in pharmaceutical, food and petroleum industry, and biological surfactant shows some advantages over the chemical surfactants, such as less toxicity, production from renewable, cheaper feedstocks and development of novel recombinant hyperproducer strains. This study is aimed to unveil the biosurfactant metabolic pathway and chemical composition in Bacillus safensis strain CCMA-560. The whole genome of the CCMA-560 strain was previously sequenced, and with the aid of bioinformatics tools, its biosurfactant metabolic pathway was compared to other pathways of closely related species. Fourier transform infrared (FTIR) and high-resolution TOF mass spectrometry (MS) were used to characterize the biosurfactant molecule. B. safensis CCMA-560 metabolic pathway is similar to other Bacillus species; however, some differences in amino acid incorporation were observed, and chemical analyses corroborated the genetic results. The strain CCMA-560 harbours two genes flanked by srfAC and srfAD not present in other Bacillus spp., which can be involved in the production of the analogue gramicidin. FTIR and MS showed that B. safensis CCMA-560 produces a mixture of at least four lipopeptides with seven amino acids incorporated and a fatty acid chain with 14 carbons, which makes this molecule similar to the biosurfactant of Bacillus pumilus, namely, pumilacidin. This is the first report on the biosurfactant production by B. safensis, encompassing the investigation of the metabolic pathway and chemical characterization of the biosurfactant molecule.

Genomic Copy Number Alterations Of Primary And Secondary Metastasizing Pleomorphic Adenomas.

Metastasizing pleomorphic adenoma (MPA) is a rare tumour, and its mechanism of metastasis still is unknown. To date, there has been no study on MPA genomics. We analysed primary and secondary MPAs with array comparative genomic hybridization to identify somatic copy number alterations and affected genes. Tumour DNA samples from primary (parotid salivary gland) and secondary (scalp skin) MPAs were subjected to array comparative genomic hybridization investigation, and the data were analysed with NEXUS COPY NUMBER DISCOVERY. The primary MPA showed copy number losses affecting 3p22.2p14.3 and 19p13.3p123, and a complex pattern of four different deletions at chromosome 6. The 3p deletion encompassed several genes: CTNNB1, SETD2, BAP1, and PBRM1, among others. The secondary MPA showed a genomic profile similar to that of the primary MPA, with acquisition of additional copy number changes affecting 9p24.3p13.1 (loss), 19q11q13.43 (gain), and 22q11.1q13.33 (gain). Our findings indicated a clonal origin of the secondary MPA, as both tumours shared a common profile of genomic copy number alterations. Furthermore, we were able to detect in the primary tumour a specific pattern of copy number alterations that could explain the metastasizing characteristic, whereas the secondary MPA showed a more unbalanced genome.

Behaviour of albino and melanic variants of Biomphalaria glabrata Say, 1818 (Mollusca: Planorbidae) following infection by Schistosoma mansoni Sambon, 1907

The behaviour of the albino and melanic variants of Biomphalaria glabrata of Belo Horizonte (MG. Brazil) was studied comparatively, in terms of their respective susceptibilities to infection by Schistosoma mansoni of the same origin, through observation of the elimination of cercariae for a three-month period and the calculation of mortality and infection rates, in control and in infected snails. The number of amoebocytes, granulocytes and hyalinocytes in the circulating hemolymph during different periods of infection was analyzed. The evolution of the infection in the tissues was observed by means of histological cross-sections. The melanic variant showed greater susceptibility to infection and a higher mortality rate. The albino variant showed a higher number of circulating amoebocytes, both granulocytes and hyalinocytes. A higher number of degenerated sporocysts were seen in the histological cross-sections of the albino variant. The results suggest that the melanic variant of B. glabrata was more susceptible to infection by S. mansoni than was the albino variant.

Emergence of Klebsiella pneumoniae carrying the novel extended-spectrum 'beta'-lactamase gene variants 'bla IND. SHV-40', 'bla IND. TEM-116' and the class I integron-associated 'bla IND. GES-7' in Brazil

A clinical Klebsiella pneumoniae isolate carrying the extended-spectrum beta-lactamase gene variants bla(SHV-40), bla(TEM-116) and bla(GES-7) was recovered. Cefoxitin and ceftazidime activity was most affected by the presence of these genes and an additional resistance to trimethoprim-sulphamethoxazole was observed. The bla(GES-7) gene was found to be inserted into a class 1 integron. These results show the emergence of novel bla(TEM) and bla(SHV) genes in Brazil. Moreover, the presence of class 1 integrons suggests a great potential for dissemination of bla(GES) genes into diverse nosocomial pathogens. Indeed, the bla(GES-7) gene was originally discovered in Enterobacter cloacae in Greece and, to our knowledge, has not been reported elsewhere

Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility ( het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer ( HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated ""gene dumps'' and, perhaps, simultaneously, as "" gene factories''.