Lifeng Tian

Genetic variants near TIMP3 and high-density lipoprotein–associated loci influence susceptibility to age-related macular degeneration

We executed a genome-wide association scan for age-related macular degeneration (AMD) in 2,157 cases and 1,150 controls. Our results validate AMD susceptibility loci near CFH (P < 10−75), ARMS2 (P < 10−59), C2/CFB (P < 10−20), C3 (P < 10−9), and CFI (P < 10−6). We compared our top findings with the Tufts/Massachusetts General Hospital genome-wide association study of advanced AMD (821 cases, 1,709 controls) and genotyped 30 promising markers in additional individuals (up to 7,749 cases and 4,625 controls). With these data, we identified a susceptibility locus near TIMP3 (overall P = 1.1 × 10−11), a metalloproteinase involved in degradation of the extracellular matrix and previously implicated in early-onset maculopathy. In addition, our data revealed strong association signals with alleles at two loci (LIPC, P = 1.3 × 10−7; CETP, P = 7.4 × 10−7) that were previously associated with high-density lipoprotein cholesterol (HDL-c) levels in blood. Consistent with the hypothesis that HDL metabolism is associated with AMD pathogenesis, we also observed association with AMD of HDL-c—associated alleles near LPL (P = 3.0 × 10−3) and ABCA1 (P = 5.6 × 10−4). Multilocus analysis including all susceptibility loci showed that 329 of 331 individuals (99%) with the highest-risk genotypes were cases, and 85% of these had advanced AMD. Our studies extend the catalog of AMD associated loci, help identify individuals at high risk of disease, and provide clues about underlying cellular pathways that should eventually lead to new therapies.

Low concordance of multiple variant-calling pipelines: practical implications for exome and genome sequencing

BackgroundTo facilitate the clinical implementation of genomic medicine by next-generation sequencing, it will be critically important to obtain accurate and consistent variant calls on personal genomes. Multiple software tools for variant calling are available, but it is unclear how comparable these tools are or what their relative merits in real-world scenarios might be.MethodsWe sequenced 15 exomes from four families using commercial kits (Illumina HiSeq 2000 platform and Agilent SureSelect version 2 capture kit), with approximately 120X mean coverage. We analyzed the raw data using near-default parameters with five different alignment and variant-calling pipelines (SOAP, BWA-GATK, BWA-SNVer, GNUMAP, and BWA-SAMtools). We additionally sequenced a single whole genome using the sequencing and analysis pipeline from Complete Genomics (CG), with 95% of the exome region being covered by 20 or more reads per base. Finally, we validated 919 single-nucleotide variations (SNVs) and 841 insertions and deletions (indels), including similar fractions of GATK-only, SOAP-only, and shared calls, on the MiSeq platform by amplicon sequencing with approximately 5000X mean coverage.ResultsSNV concordance between five Illumina pipelines across all 15 exomes was 57.4%, while 0.5 to 5.1% of variants were called as unique to each pipeline. Indel concordance was only 26.8% between three indel-calling pipelines, even after left-normalizing and intervalizing genomic coordinates by 20 base pairs. There were 11% of CG variants falling within targeted regions in exome sequencing that were not called by any of the Illumina-based exome analysis pipelines. Based on targeted amplicon sequencing on the MiSeq platform, 97.1%, 60.2%, and 99.1% of the GATK-only, SOAP-only and shared SNVs could be validated, but only 54.0%, 44.6%, and 78.1% of the GATK-only, SOAP-only and shared indels could be validated. Additionally, our analysis of two families (one with four individuals and the other with seven), demonstrated additional accuracy gained in variant discovery by having access to genetic data from a multi-generational family.ConclusionsOur results suggest that more caution should be exercised in genomic medicine settings when analyzing individual genomes, including interpreting positive and negative findings with scrutiny, especially for indels. We advocate for renewed collection and sequencing of multi-generational families to increase the overall accuracy of whole genomes.

Analysis of 589,306 genomes identifies individuals resilient to severe Mendelian childhood diseases

Genetic studies of human disease have traditionally focused on the detection of disease-causing mutations in afflicted individuals. Here we describe a complementary approach that seeks to identify healthy individuals resilient to highly penetrant forms of genetic childhood disorders. A comprehensive screen of 874 genes in 589,306 genomes led to the identification of 13 adults harboring mutations for 8 severe Mendelian conditions, with no reported clinical manifestation of the indicated disease. Our findings demonstrate the promise of broadening genetic studies to systematically search for well individuals who are buffering the effects of rare, highly penetrant, deleterious mutations. They also indicate that incomplete penetrance for Mendelian diseases is likely more common than previously believed. The identification of resilient individuals may provide a first step toward uncovering protective genetic variants that could help elucidate the mechanisms of Mendelian diseases and new therapeutic strategies.

Mutations in PDGFRB cause autosomal-dominant infantile myofibromatosis.

Infantile myofibromatosis (IM) is a disorder of mesenchymal proliferation characterized by the development of nonmetastasizing tumors in the skin, muscle, bone, and viscera. Occurrence within families across multiple generations is suggestive of an autosomal-dominant (AD) inheritance pattern, but autosomal-recessive (AR) modes of inheritance have also been proposed. We performed whole-exome sequencing (WES) in members of nine unrelated families clinically diagnosed with AD IM to identify the genetic origin of the disorder. In eight of the families, we identified one of two disease-causing mutations, c.1978C>A (p.Pro660Thr) and c.1681C>T (p.Arg561Cys), in PDGFRB. Intriguingly, one family did not have either of these PDGFRB mutations but all affected individuals had a c.4556T>C (p.Leu1519Pro) mutation in NOTCH3. Our studies suggest that mutations in PDGFRB are a cause of IM and highlight NOTCH3 as a candidate gene. Further studies of the crosstalk between PDGFRB and NOTCH pathways may offer new opportunities to identify mutations in other genes that result in IM and is a necessary first step toward understanding the mechanisms of both tumor growth and regression and its targeted treatment.

The Cell Fate Determination Factor Dachshund Inhibits Androgen Receptor Signaling and Prostate Cancer Cellular Growth

Initially isolated as the dominant suppressor of the mutant epidermal growth factor receptor (ellipse), the Dachshund gene plays a key role in metazoan development regulating the Retinal Determination Gene Network. Herein, the DACH1 gene was expressed in normal prostate epithelial cells with reduced expression in human prostate cancer. DACH1 inhibited prostate cancer cellular DNA synthesis, growth in colony forming assays, and blocked contact-independent growth in soft agar assays. DACH1 inhibited androgen receptor (AR) activity, requiring a conserved DS Domain (Dachshund domain conserved with Ski/Sno) that bound NCoR/HDAC and was recruited to an androgen-responsive gene promoter. DACH1 inhibited ligand-dependent activity of AR mutations identified in patients with androgen-insensitive prostate cancer. The DS domain was sufficient for repression of the AR wild-type but failed to repress an AR acetylation site point mutant. These studies show a role for the Retinal Determination Gene Network in regulating cellular growth and signaling in prostate cancer.

Comprehensive analysis of gene expression in human retina and supporting tissues

Understanding the influence of gene expression on the molecular mechanisms underpinning human phenotypic diversity is fundamental to being able to predict health outcomes and treat disease. We have carried out whole transcriptome expression analysis on a series of eight normal human postmortem eyes by RNA sequencing. Here we present data showing that ∼80% of the transcriptome is expressed in the posterior layers of the eye and that there is significant differential expression not only between the layers of the posterior part of the eye but also between locations of a tissue layer. These differences in expression also extend to alternative splicing and splicing factors. Differentially expressed genes are enriched for genes associated with psychiatric, immune and cardiovascular disorders. Enrichment categories for gene ontology included ion transport, synaptic transmission and visual and sensory perception. Lastly, allele-specific expression was found to be significant forCFH,C3 andCFB, which are known risk genes for age-related macular degeneration. These expression differences should be useful in determining the underlying biology of associations with common diseases of the human retina, retinal pigment epithelium and choroid and in guiding the analysis of the genomic regions involved in the control of normal gene expression.

Analysis of Six Genetic Risk Factors Highly Associated with AMD in the Region Surrounding ARMS2 and HTRA1 on Chromosome 10, Region q26

Purpose. To determine the relationship of six genetic variants (rs10490924, rs3750848, del443ins54, rs3793917, rs11200638, and rs932275) localized to the ARMS2-HTRA1 region of chromosome 10, region q26, as risk factors for age-related macular degeneration (AMD), to define the haplotype structure of these six loci, and to confirm their genetic association with the disease. Methods. Caucasian patients (n = 482) were stratified into categories based on AREDS (Age-Related Eye Disease Study) grading criteria (groups 0 and 1 served as the control, groups 3 and 4 contained subjects with AMD, and group 2 was excluded from the analysis). The six genetic variants in the ARMS2-HTRA1 region were genotyped and analyzed both independently and as a joint haplotype for association in subjects with disease (n = 291) compared with the control (n = 191). Results. The six high-risk alleles all showed a statistically significant association with AMD (the most significant SNP was rs10490924 [P < or = 3.31 x 10(-5), OR = 1.86]; the least significant SNP was rs932275 [P < or = 9.15 x 10(-5), OR = 1.78]). Multimarker analysis revealed that all six markers were in strong linkage disequilibrium with each other, and the two major haplotypes that captured >98% of the genetic variation in the region were both significantly associated with the disease: One increased the risk of AMD and contained only risk alleles (P < or = 2.20 x 10(-5)), and the other haplotype decreased the risk of AMD and contained only wild-type alleles (P < or = 6.81 x 10(-5)). Furthermore, 36 individuals comprising both cases and controls were identified outside of these two major haplotypes, with at least one discordant marker. Conclusions. The results replicate the previously reported association between the high-risk alleles and AMD and independently confirm, for the first time, an association with AMD and the indel (del443ins54) polymorphism in a Caucasian population. Two major haplotypes that are associated with AMD and many minor novel haplotypes were identified. The novel haplotypes, identified from 36 cases and controls with discordant alleles spanning the ARMS2-HTRA1 region provide unique opportunities to gauge the relative phenotypic contributions of each of these genetic risk factors. With the identification of more discordant patients in the future, it may be possible to resolve the ongoing controversy as to which of the risk alleles and genes (ARMS2 vs. HTRA1) has the greatest impact on disease susceptibility. Future work should include the analysis of larger and more diverse populations, to further define the linkage structure of the region with a focus on phenotypic effects on AMD of the various haplotypes involving 10q26, as well as a functional analysis of the normal ARMS2 protein.

EYA1 Phosphatase Function Is Essential To Drive Breast Cancer Cell Proliferation Through Cyclin D1

The Drosophila Eyes Absent Homologue 1 (EYA1) is a component of the retinal determination gene network and serves as an H2AX phosphatase. The cyclin D1 gene encodes the regulatory subunits of a holoenzyme that phosphorylates and inactivates the pRb protein. Herein, comparison with normal breast showed that EYA1 is overexpressed with cyclin D1 in luminal B breast cancer subtype. EYA1 enhanced breast tumor growth in mice in vivo, requiring the phosphatase domain. EYA1 enhanced cellular proliferation, inhibited apoptosis, and induced contact-independent growth and cyclin D1 abundance. The induction of cellular proliferation and cyclin D1 abundance, but not apoptosis, was dependent upon the EYA1 phosphatase domain. The EYA1-mediated transcriptional induction of cyclin D1 occurred via the AP-1-binding site at -953 and required the EYA1 phosphatase function. The AP-1 mutation did not affect SIX1-dependent activation of cyclin D1. EYA1 was recruited in the context of local chromatin to the cyclin D1 AP-1 site. The EYA1 phosphatase function determined the recruitment of CBP, RNA polymerase II, and acetylation of H3K9 at the cyclin D1 gene AP-1 site regulatory region in the context of local chromatin. The EYA1 phosphatase regulates cell-cycle control via transcriptional complex formation at the cyclin D1 promoter.

Acetylation and nuclear receptor action

Acetylation is an essential post-translational modification featuring an acetyl group that is covalently conjugated to a protein substrate. Histone acetylation was first proposed nearly half a century ago by Dr. Vincent Allfrey. Subsequent studies have shown that the acetylated core histones are often associated with transcriptionally active chromatin. Acetylation at lysine residues of histone tails neutralizes the positive charge, which decreases their binding ability to DNA and increases the accessibility of transcription factors and coactivators to the chromatin template. In addition to histones, a number of non-histone substrates are acetylated. Acetylation of non-histone proteins governs biological processes, such as cellular proliferation and survival, transcriptional activity, and intracellular trafficking. We demonstrated that acetylation of transcription factors can regulate cellular growth. Furthermore, we showed that nuclear receptors (NRs) are acetylated at a phylogenetically conserved motif. Since our initial observations with the estrogen and androgen receptors, more than a dozen NRs have been shown to function as substrates for acetyltransferases with diverse functional consequences. This review focuses on the acetylation of NRs and the effect of acetylation on NR function. We discuss the potential role of acetylation in disease initiation and progression with an emphasis on tumorigenesis.

Activating Peroxisome Proliferator-Activated Receptor {gamma} Mutant Promotes Tumor Growth In vivo by Enhancing Angiogenesis

Peroxisome proliferator-activated receptor gamma (PPARgamma) is expressed in a variety of cancer cells. The addition of ligand activates the receptor by inducing a conformational change in the receptor, which can be recapitulated by mutation. To investigate the role of activated PPARgamma signaling in breast cancer, we compared the function of a constitutively active PPARgamma (PgammaCA) mutant with the wild-type PPARgamma in ErbB2-induced mammary tumorigenesis in vivo. Tumor cells transduced with either PPARgamma or PgammaCA were implanted into immunocompetent FVB mice. Enhanced tumor growth was observed in PgammaCA-transduced cells, which was associated with increased angiogenesis and endothelial stem cells as evidenced by increased number of cells stained with von Willebrand factor, c-Kit, CD133, and CD31. Genome-wide expression profiling identified a group of genes within the angiogenesis pathway, including Angptl4, as targets of activated PPARgamma; PgammaCA also induced Angptl4 protein secretion in ErbB2-transformed mammary epithelial cells. Angptl4 promoted vascular endothelial cell migration; conversely, immunodepletion of Angptl4 reduced PgammaCA-mediated cellular migration. Collectively, these studies suggest that activated PPARgamma induces Angptl4 to promote tumor growth through enhanced angiogenesis in vivo.