Examining Emergence of Functional Gene Clustering in a Simulated Evolution

Recent research suggests that rather than being random, gene order may be coupled with gene functionality. These findings may be explained by mechanisms that require physical proximity such as co-expression and co-regulation. Alternatively, they may be due to evolutionary-dynamics forces, as expressed in genetic drift or linkage disequilibrium. This paper proposes a biologically plausible model for evolutionary development. Using the model, which includes natural selection and the development of gene networks and cellular organisms, the co-evolution of recombination rate and gene functionality is examined. The results presented here are compatible with previous biological findings showing that functionally related genes are clustered. These results imply that evolutionary pressure in a complex environment is sufficient for the emergence of gene order that is coupled with functionality. They shed further light on the mechanisms that may cause such gene clusters.     

Missense Variants Reveal Functional Insights Into the Human ARID Family of Gene Regulators

Missense variants are alterations to protein coding sequences that result in amino acid substitutions. They can be deleterious if the amino acid is required for maintaining structure or/and function, but are likely to be tolerated at other sites. Consequently, missense variation within a healthy population can mirror the effects of negative selection on protein structure and function, such that functional sites on proteins are often depleted of missense variants. Advances in high-throughput sequencing have dramatically increased the sample size of available human variation data, allowing for population-wide analysis of selective pressures. In this study, we developed a convenient set of tools, called 1D-to-3D, for visualizing the positions of missense variants on protein sequences and structures. We used these tools to characterize human homologues of the ARID family of gene regulators. ARID family members are implicated in multiple cancer types, developmental disorders, and immunological diseases but current understanding of their mechanistic roles is incomplete. Combined with phylogenetic and structural analyses, our approach allowed us to characterise sites important for protein-protein interactions, histone modification recognition, and DNA binding by the ARID proteins. We find that comparing missense depletion patterns among paralogs can reveal sub-functionalization at the level of domains. We propose that visualizing missense variants and their depletion on structures can serve as a valuable tool for complementing evolutionary and experimental findings.     

Functional Coding Variants in SLC6A15, a Possible Risk Gene for Major Depression

SLC6A15 is a neuron-specific neutral amino acid transporter that belongs to the solute carrier 6 gene family. This gene family is responsible for presynaptic re-uptake of the majority of neurotransmitters. Convergent data from human studies, animal models and pharmacological investigations suggest a possible role of SLC6A15 in major depressive disorder. In this work, we explored potential functional variants in this gene that could influence the activity of the amino acid transporter and thus downstream neuronal function and possibly the risk for stress-related psychiatric disorders. DNA from 400 depressed patients and 400 controls was screened for genetic variants using a pooled targeted re-sequencing approach. Results were verified by individual re-genotyping and validated non-synonymous coding variants were tested in an independent sample (N = 1934). Nine variants altering the amino acid sequence were then assessed for their functional effects by measuring SLC6A15 transporter activity in a cellular uptake assay. In total, we identified 405 genetic variants, including twelve non-synonymous variants. While none of the non-synonymous coding variants showed significant differences in case-control associations, two rare non-synonymous variants were associated with a significantly increased maximal ³H proline uptake as compared to the wildtype sequence. Our data suggest that genetic variants in the SLC6A15 locus change the activity of the amino acid transporter and might thus influence its neuronal function and the risk for stress-related psychiatric disorders. As statistically significant association for rare variants might only be achieved in extremely large samples (N >70,000) functional exploration may shed light on putatively disease-relevant variants.     

Polymorphic Variants of Extracellular Superoxide Dismutase Gene in a Romanian Population with Atheroma

Extracellular superoxide dismutase (EC-SOD) is the main SOD isoform in the arterial wall contributing to cardiovascular defense against oxidative stress by removing the superoxide anion. In our study, the Thr40Ala and Arg213Gly polymorphic variants of the EC-SOD gene (SOD ( 3 )) were investigated for associations with atherosclerosis and other related factors in 144 subjects with significant atheroma (having one, two, or three major coronary arteries with >50% obstruction, and/or peripheral artery lesions, and/or carotid artery stenosis demonstrated by angiography and echography) and in 150 subjects with no significant atheroma. For the Arg213Gly polymorphism, only five heterozygous subjects were found. Although the difference in the genotype distribution for the Thr40Ala polymorphism was not statistically significant between patients with atheroma (AA 49.3%, AG 34.7%, GG 16.0%) and those without significant atheroma (AA 41.3%, AG 43.3%, GG 15.3%), there was an association of the Thr40 allele with diabetes (P = 0.03) and hypertension (P = 0.04).     

An Unbiased Method for Clustering Bacterial Effectors Using Host Cellular Phenotypes

We present a novel method implementing unbiased high-content morphometric cell analysis to classify bacterial effector phenotypes. This clustering methodology represents a significant advance over more qualitative visual approaches and can also be used to classify, and therefore predict the likely function of, unknown effector genes from any microbial genome. As a proof of concept, we use this approach to investigate 23 genetic regions predicted to encode antimacrophage effectors located across the genome of the insect and human pathogen Photorhabdus asymbiotica. Statistical cluster analysis using multiple cellular measures categorized treated macrophage phenotypes into three major groups relating to their putative functionality: (i) adhesins, (ii) cytolethal toxins, and (iii) cytomodulating toxins. Further investigation into their effects on phagocytosis revealed that several effectors also modulate this function and that the nature of this modulation (increased or decreased phagocytosis) is linked to the phenotype cluster group. Categorizing potential functionalities in this way allows rapid functional follow-up of key candidates for more-directed cell biological or biochemical investigation. Such an unbiased approach to the classification of candidate effectors will be useful for describing virulence-related regions in a wide range of genomes and will be useful in assigning putative functions to the growing number of microbial genes whose function remains unclear from homology searching.     

Anti-Candidal Activity of Genetically Engineered Histatin Variants with Multiple Functional Domains

The human bodily defense system includes a wide variety of innate antimicrobial proteins. Histatins are small molecular weight proteins produced by the human salivary glands that exhibit antifungal and antibacterial activities. While evolutionarily old salivary proteins such as mucins and proline-rich proteins contain large regions of tandem repeats, relatively young proteins like histatins do not contain such repeated domains. Anticipating that domain duplications have a functional advantage, we genetically engineered variants of histatin 3 with one, two, three, or four copies of the functional domain by PCR and splice overlap. The resulting proteins, designated reHst3 1-mer, reHist3 2-mer, reHis3 3-mer and reHist3 4-mer, exhibited molecular weights of 4,062, 5,919, 7,777, and 9,634 Da, respectively. The biological activities of these constructs were evaluated in fungicidal assays toward Candida albicans blastoconidia and germinated cells. The antifungal activities per mole of protein increased concomitantly with the number of functional domains present. This increase, however, was higher than could be anticipated from the molar concentration of functional domains present in the constructs. The demonstrated increase in antifungal activity may provide an evolutionary explanation why such domain multiplication is a frequent event in human salivary proteins.     

Genetic and Molecular Functional Characterization of Variants within TNFSF13B,a Positional Candidate Preeclampsia Susceptibility Gene on 13q

Background

Preeclampsia is a serious pregnancy complication, demonstrating a complex pattern of inheritance. The elucidation of genetic liability to preeclampsia remains a major challenge in obstetric medicine. We have adopted a positional cloning approach to identify maternal genetic components, with linkages previously demonstrated to chromosomes 2q, 5q and 13q in an Australian/New Zealand familial cohort. The current study aimed to identify potential functional and structural variants in the positional candidate gene TNFSF13B under the 13q linkage peak and assess their association status with maternal preeclampsia genetic susceptibility.

Methodology/Principal Findings

The proximal promoter and coding regions of the positional candidate gene TNFSF13B residing within the 13q linkage region was sequenced using 48 proband or founder individuals from Australian/New Zealand families. Ten sequence variants (nine SNPs and one single base insertion) were identified and seven SNPs were successfully genotyped in the total Australian/New Zealand family cohort (74 families/480 individuals). Borderline association to preeclampsia (p = 0.0153) was observed for three rare SNPs (rs16972194, rs16972197 and rs56124946) in strong linkage disequilibrium with each other. Functional evaluation by electrophoretic mobility shift assays showed differential nuclear factor binding to the minor allele of the rs16972194 SNP, residing upstream of the translation start site, making this a putative functional variant. The observed genetic associations were not replicated in a Norwegian case/control cohort (The Nord-Trøndelag Health Study (HUNT2), 851 preeclamptic and 1,440 non-preeclamptic women).

Conclusion/Significance

TNFSF13B has previously been suggested to contribute to the normal immunological adaption crucial for a successful pregnancy. Our observations support TNFSF13B as a potential novel preeclampsia susceptibility gene. We discuss a possible role for TNFSF13B in preeclampsia pathogenesis, and propose the rs16972194 variant as a candidate for further functional evaluation.     

Interpretation of Genomic Variants Using a Unified Biological Network Approach

The decreasing cost of sequencing is leading to a growing repertoire of personal genomes. However, we are lagging behind in understanding the functional consequences of the millions of variants obtained from sequencing. Global system-wide effects of variants in coding genes are particularly poorly understood. It is known that while variants in some genes can lead to diseases, complete disruption of other genes, called ‘loss-of-function tolerant’, is possible with no obvious effect. Here, we build a systems-based classifier to quantitatively estimate the global perturbation caused by deleterious mutations in each gene. We first survey the degree to which gene centrality in various individual networks and a unified ‘Multinet’ correlates with the tolerance to loss-of-function mutations and evolutionary conservation. We find that functionally significant and highly conserved genes tend to be more central in physical protein-protein and regulatory networks. However, this is not the case for metabolic pathways, where the highly central genes have more duplicated copies and are more tolerant to loss-of-function mutations. Integration of three-dimensional protein structures reveals that the correlation with centrality in the protein-protein interaction network is also seen in terms of the number of interaction interfaces used. Finally, combining all the network and evolutionary properties allows us to build a classifier distinguishing functionally essential and loss-of-function tolerant genes with higher accuracy (AUC = 0.91) than any individual property. Application of the classifier to the whole genome shows its strong potential for interpretation of variants involved in Mendelian diseases and in complex disorders probed by genome-wide association studies.     

New Gene Variants Associated with the Risk of Chronic HBV Infection

Some patients with chronic hepatitis B virus(HBV) infection failed to clear HBV, even persistently continue to produce antibodies to HBV. Here we performed a two stage genome wide association study in a cohort of Chinese patients designed to discover single nucleotide variants that associate with HBV infection and clearance of HBV. The first stage involved genome wide exome sequencing of 101 cases(HBsAg plus anti-HBs positive) compared with 102 control patients(antiHBs positive, HBsAg negative). Over 80% of individual sequences displayed 209 sequence coverage. Adapters,uncertain bases [10% or low-quality base calls([50%) were filtered and compared to the human reference genome hg19. In the second stage, 579 chronic HBV infected cases and 439 HBV clearance controls were sequenced with selected genes from the first stage. Although there were no significant associated gene variants in the first stage, two significant gene associations were discovered when the two stages were assessed in a combined analysis. One association showed rs506121-‘‘T' allele [within the dedicator of cytokinesis 8(DOCK8) gene] was higher in chronic HBV infection group than that in clearance group(P = 0.002, OR = 0.77, 95% CI [0.65, 0.91]). The second association involved rs2071676—A allele within the Carbonic anhydrase(CA9) gene that was significantly elevated in chronic HBV infection group compared to the clearance group(P = 0.0003, OR = 1.35, 95% CI [1.15, 1.58]). Upon replication these gene associations would suggest the influence of DOCK8 and CA9 as potential risk genetic factors in the persistence of HBV infection.     

CFTR Gene Variants and Genotypes in Russian Patients with CBAVD Syndrome

Russian Journal of Genetics - CBAVD syndrome is one of the common genetic causes of male infertility, associated with obstructive azoospermia, commonly resulting from pathogenic CFTR gene variants....     

Type II Transmembrane Serine Protease Gene Variants Associate with Breast Cancer

Type II transmembrane serine proteases (TTSPs) are related to tumor growth, invasion, and metastasis in cancer. Genetic variants in these genes may alter their function, leading to cancer onset and progression, and affect patient outcome. Here, 464 breast cancer cases and 370 controls were genotyped for 82 single-nucleotide polymorphisms covering eight genes. Association of the genotypes was estimated against breast cancer risk, breast cancer–specific survival, and survival in different treatment groups, and clinicopathological variables. SNPs in TMPRSS3 (rs3814903 and rs11203200), TMPRSS7 (rs1844925), and HGF (rs5745752) associated significantly with breast cancer risk (P_trend = 0.008–0.042). SNPs in TMPRSS1 (rs12151195 and rs12461158), TMPRSS2 (rs2276205), TMPRSS3 (rs3814903), and TMPRSS7 (rs2399403) associated with prognosis (P = 0.004–0.046). When estimating the combined effect of the variants, the risk of breast cancer was higher with 4–5 alleles present compared to 0–2 alleles (P = 0.0001; OR, 2.34; 95% CI, 1.39–3.94). Women with 6–8 survival-associating alleles had a 3.3 times higher risk of dying of breast cancer compared to women with 1–3 alleles (P = 0.001; HR, 3.30; 95% CI, 1.58–6.88). The results demonstrate the combined effect of variants in TTSPs and their related genes in breast cancer risk and patient outcome. Functional analysis of these variants will lead to further understanding of this gene family, which may improve individualized risk estimation and development of new strategies for treatment of breast cancer.     

EGFR Gene Variants Are Associated with Specific Somatic Aberrations in Glioma

A number of gene variants have been associated with an increased risk of developing glioma. We hypothesized that the reported risk variants may be associated with tumor genomic instability. To explore potential correlations between germline risk variants and somatic genetic events, we analyzed matched tumor and blood samples from 95 glioma patients by means of SNP genotyping. The generated genotype data was used to calculate genome-wide allele-specific copy number profiles of the tumor samples. We compared the copy number profiles across samples and found two EGFR gene variants (rs17172430 and rs11979158) that were associated with homozygous deletion at the CDKN2A/B locus. One of the EGFR variants (rs17172430) was also associated with loss of heterozygosity at the EGFR locus. Our findings were confirmed in a separate dataset consisting of matched blood and tumor samples from 300 glioblastoma patients, compiled from publically available TCGA data. These results imply there is a functional effect of germline EGFR variants on tumor progression.     

Genetic Differentiation of Transposable Elements under Mutation and Unbiased Gene Conversion

A model is developed to predict the extent of genetic differentiation in a family of transposable elements under the combined effects of genetic drift, transposition, mutation and unbiased gene conversion. The model is based on simplifying assumptions that are valid when transposition is always to new sites and copy number per site is low. In the absence of gene conversion, the degree of differentiation as measured by the probability of identity of different elements is the same as at a single locus with the same mutation rate but in a population of effective size Nc/2, where N is the population size and c is the number of copies per individual. The inclusion of unbiased gene conversion does not significantly change this result. If, as seems to be the case, families of transposable elements are relatively homogeneous, then the model implies either that mutation rates for transposable elements are much lower than at comparable single-copy loci or that some other force, such as natural selection or biased gene conversion, is at work. Transposition is a very ineffective force for homogenizing a family of transposable elements.     

Functional Analysis of Prognostic Gene Expression Network Genes in Metastatic Breast Cancer Models

Identification of conserved co-expression networks is a useful tool for clustering groups of genes enriched for common molecular or cellular functions [1]. The relative importance of genes within networks can frequently be inferred by the degree of connectivity, with those displaying high connectivity being significantly more likely to be associated with specific molecular functions [2]. Previously we utilized cross-species network analysis to identify two network modules that were significantly associated with distant metastasis free survival in breast cancer. Here, we validate one of the highly connected genes as a metastasis associated gene. Tpx2, the most highly connected gene within a proliferation network specifically prognostic for estrogen receptor positive (ER+) breast cancers, enhances metastatic disease, but in a tumor autonomous, proliferation-independent manner. Histologic analysis suggests instead that variation of TPX2 levels within disseminated tumor cells may influence the transition between dormant to actively proliferating cells in the secondary site. These results support the co-expression network approach for identification of new metastasis-associated genes to provide new information regarding the etiology of breast cancer progression and metastatic disease.