Towards a Genetic Linkage Map of the California Condor, an Endangered New World Vulture Species.

The development of a linkage map is an important component for promoting genetic and genomic studies in California condors, an endangered New World vulture species. Using a set of designed anonymous microsatellite markers, we genotyped a reference condor population involving 121 individuals. After marker validation and genotype filtering, the genetic linkage analysis was performed using 123 microsatellite loci. This resulted in the identification of 15 linkage groups/subgroups that formed a first-generation condor genetic map, while no markers linked to a lethal chondrodystrophy mutation were found. A panel of polymorphic markers that is instrumental in molecular parentage diagnostics and other genetic studies in the California condor was selected. Further condor conservation genomics research will be focused on updating the linkage map and integrating it with cytogenetic and BAC-based physical maps and ultimately with the genome sequence assembly.

Facultative Parthenogenesis in California Condors.

Parthenogenesis is a relatively rare event in birds, documented in unfertilized eggs from columbid, galliform, and passerine females with no access to males. In the critically endangered California condor, parentage analysis conducted utilizing polymorphic microsatellite loci has identified two instances of parthenogenetic development from the eggs of two females in the captive breeding program, each continuously housed with a reproductively capable male with whom they had produced offspring. Paternal genetic contribution to the two chicks was excluded. Both parthenotes possessed the expected male ZZ sex chromosomes and were homozygous for all evaluated markers inherited from their dams. These findings represent the first molecular marker-based identification of facultative parthenogenesis in an avian species, notably of females in regular contact with fertile males, and add to the phylogenetic breadth of vertebrate taxa documented to have reproduced via asexual reproduction.

MicroRNA-related polymorphisms and non-Hodgkin lymphoma susceptibility in the Multicenter AIDS Cohort Study.

BACKGROUND: MicroRNAs, small non-coding RNAs involved in gene regulation, are implicated in lymphomagenesis. We evaluated whether genetic variations in microRNA coding regions, binding sites, or biogenesis genes (collectively referred to as miRNA-SNPs) were associated with risk of AIDS-associated non-Hodgkin lymphoma (AIDS-NHL), and serum levels of four lymphoma-related microRNAs. METHODS: Twenty-five miRNA-SNPs were genotyped in 180 AIDS-NHL cases and 529 HIV-infected matched controls from the Multicenter AIDS Cohort Study (MACS), and real-time polymerase chain reaction was used to quantify serum microRNA levels. Adjusted odds ratios (ORs) estimated using conditional logistic regression evaluated associations between miRNA-SNPs and AIDS-NHL risk. A semi-Bayes shrinkage approach was employed to reduce likelihood of false-positive associations. Adjusted mean ratios (MR) calculated using linear regression assessed associations between miRNA-SNPs and serum microRNA levels. RESULTS: DDX20 rs197412, a non-synonymous miRNA biogenesis gene SNP, was associated with AIDS-NHL risk (OR=1.34 per minor allele; 95% CI: 1.02-1.75), and higher miRNA-222 serum levels nearing statistical significance (MR=1.21 per minor allele; 95% CI: 0.98-1.49). MiRNA-196a2 rs11614913 was associated with decreased central nervous system (CNS) AIDS-NHL (CT vs. CC OR=0.52; 95% CI: 0.27-0.99). The minor allele of HIF1A rs2057482, which creates a miRNA-196a2 binding site, was associated with systemic AIDS-NHL risk (OR=1.73 per minor allele; 95% CI: 1.12-2.67), and decreased CNS AIDS-NHL risk (OR=0.49 per minor allele; 95% CI: 0.25-0.94). CONCLUSIONS: This study suggests that a few miRNA-SNPs are associated with AIDS-NHL risk and may modulate miRNA expression. These results support a role for miRNA in AIDS-NHL and may highlight pathways to be targeted for risk stratification or therapeutics.

Shared HLA Class I and II Alleles and Clonally Restricted Public and Private Brain-Infiltrating αß T Cells in a Cohort of Rasmussen Encephalitis Surgery Patients.

Rasmussen encephalitis (RE) is a rare pediatric neuroinflammatory disease characterized by intractable seizures and unilateral brain atrophy. T cell infiltrates in affected brain tissue and the presence of circulating autoantibodies in some RE patients have indicated that RE may be an autoimmune disease. The strongest genetic links to autoimmunity reside in the MHC locus, therefore, we determined the human leukocyte antigen (HLA) class I and class II alleles carried by a cohort of 24 RE surgery cases by targeted in-depth genomic sequencing. Compared with a reference population the allelic frequency of three alleles, DQA1*04:01:01, DQB1*04:02:01, and HLA-C*07:02:01:01 indicated that they might confer susceptibility to the disease. It has been reported that HLA-C*07:02 is a risk factor for Graves disease. Further, eight patients in the study cohort carried HLA-A*03:01:01:01, which has been linked to susceptibility to multiple sclerosis. Four patients carried a combination of three HLA class II alleles that has been linked to type 1 diabetes (DQA1*05:01:01:01~DQB1*02:01:01~DRB1*03:01:01:01), and five patients carried a combination of HLA class II alleles that has been linked to the risk of contracting multiple sclerosis (DQA1*01:02:01:01, DQB1*06:02:01, DRB1*15:01:01:01). We also analyzed the diversity of αß T cells in brain and blood specimens from 14 of these RE surgery cases by sequencing the third complementarity regions (CDR3s) of rearranged T cell receptor ß genes. A total of 31 unique CDR3 sequences accounted for the top 5% of all CDR3 sequences in the 14 brain specimens. Thirteen of these sequences were found in sequencing data from healthy blood donors; the remaining 18 sequences were patient specific. These observations provide evidence for the clonal expansion of public and private T cells in the brain, which might be influenced by the RE patient's HLA haplotype.

Association of COMT and TPH-2 genes with DSM-5 based PTSD symptoms.

BACKGROUND: Dopaminergic and serotonergic systems have been implicated in PTSD. The present study evaluated the association of four catechol-O-methyltransferase (COMT) gene loci, and the joint effect of COMT and tryptophan hydroxylase 2 (TPH2) genes on PTSD symptoms. METHODS: Subjects included 200 Caucasian Armenian adults exposed to the 1988 Spitak earthquake from 12 multigenerational (3-5 generations) families. Instruments used included the UCLA PTSD Reaction Index based on DSM-5 criteria, and the Beck Depression Inventory. RESULTS: The adjusted heritabilitiy of vulnerability to DSM-5 based PTSD symptoms was 0.60 (p<10(-4)). There was a significant association of the COMT allele rs4633C with total PTSD (p<0.03), and D category (p<0.04) (negative alterations in cognitions and mood) severity scores, but not with C category (avoidance) scores. There was no genetic correlation between C and D category severity scores. COMT allele rs4633C and the TPH-2 allele rs11178997T together accounted for 7% of the variance in PTSD severity scores (p<0.001). None of the COMT alleles were associated with depression. LIMITATIONS: The ratings of earthquake exposure and prior trauma may have been subject to recall bias. The findings may not be generalizable to other ethnic/racial populations. CONCLUSION: COMT allele rs4633C may be causally related and/or is in linkage disequilibrium with gene(s) that are causally related to PTSD symptoms. Carriers of these COMT and the TPH-2 alleles may be at increased risk for PTSD. The findings provide biological support for dividing DSM-IV category C symptoms into DSM-5 categories C and D.

Rare Genomic Variants Link Bipolar Disorder with Anxiety Disorders to CREB-Regulated Intracellular Signaling Pathways.

Bipolar disorder is a common, complex, and severe psychiatric disorder with cyclical disturbances of mood and a high suicide rate. Here, we describe a family with four siblings, three affected females and one unaffected male. The disease course was characterized by early-onset bipolar disorder and co-morbid anxiety spectrum disorders that followed the onset of bipolar disorder. Genetic risk factors were suggested by the early onset of the disease, the severe disease course, including multiple suicide attempts, and lack of adverse prenatal or early life events. In particular, drug and alcohol abuse did not contribute to the disease onset. Exome sequencing identified very rare, heterozygous, and likely protein-damaging variants in eight brain-expressed genes: IQUB, JMJD1C, GADD45A, GOLGB1, PLSCR5, VRK2, MESDC2, and FGGY. The variants were shared among all three affected family members but absent in the unaffected sibling and in more than 200 controls. The genes encode proteins with significant regulatory roles in the ERK/MAPK and CREB-regulated intracellular signaling pathways. These pathways are central to neuronal and synaptic plasticity, cognition, affect regulation and response to chronic stress. In addition, proteins in these pathways are the target of commonly used mood-stabilizing drugs, such as tricyclic antidepressants, lithium, and valproic acid. The combination of multiple rare, damaging mutations in these central pathways could lead to reduced resilience and increased vulnerability to stressful life events. Our results support a new model for psychiatric disorders, in which multiple rare, damaging mutations in genes functionally related to a common signaling pathway contribute to the manifestation of bipolar disorder.

Therapeutic potential of the poly(ADP-ribose) polymerase inhibitor rucaparib for the treatment of sporadic human ovarian cancer.

Here, we investigate the potential role of the PARP inhibitor rucaparib (CO-338, formerly known as AG014699 and PF-01367338) for the treatment of sporadic ovarian cancer. We studied the growth inhibitory effects of rucaparib in a panel of 39 ovarian cancer cell lines that were each characterized for mutation and methylation status of BRCA1/2, baseline gene expression signatures, copy number variations of selected genes, PTEN status, and sensitivity to platinum-based chemotherapy. To study interactions with chemotherapy, we used multiple drug effect analyses and assessed apoptosis, DNA fragmentation, and γH2AX formation. Concentration-dependent antiproliferative effects of rucaparib were seen in 26 of 39 (67%) cell lines and were not restricted to cell lines with BRCA1/2 mutations. Low expression of other genes involved in homologous repair (e.g., BCCIP, BRCC3, ATM, RAD51L1), amplification of AURKA or EMSY, and response to platinum-based chemotherapy was associated with sensitivity to rucaparib. Drug interactions with rucaparib were synergistic for topotecan, synergistic, or additive for carboplatin, doxorubicin or paclitaxel, and additive for gemcitabine. Synergy was most pronounced when rucaparib was combined with topotecan, which resulted in enhanced apoptosis, DNA fragmentation, and γH2AX formation. Importantly, rucaparib potentiated chemotherapy independent of its activity as a single agent. PARP inhibition may be a useful therapeutic strategy for a wider range of ovarian cancers bearing deficiencies in the homologous recombination pathway other than just BRCA1/2 mutations. These results support further clinical evaluation of rucaparib either as a single agent or as an adjunct to chemotherapy for the treatment of sporadic ovarian cancer.

Systematic identification of H274Y compensatory mutations in influenza A virus neuraminidase by high-throughput screening.

Compensatory mutations contribute to the appearance of the oseltamivir resistance substitution H274Y in the neuraminidase (NA) gene of H1N1 influenza viruses. Here, we describe a high-throughput screening method utilizing error-prone PCR and next-generation sequencing to comprehensively screen NA genes for H274Y compensatory mutations. We found four mutations that can either fully (R194G, E214D) or partially (L250P, F239Y) compensate for the fitness deficiency of the H274Y mutant. The compensatory effect of E214D is applicable in both seasonal influenza virus strain A/New Caledonia/20/1999 and 2009 pandemic swine influenza virus strain A/California/04/2009. The technique described here has the potential to profile a gene at the single-nucleotide level to comprehend the dynamics of mutation space and fitness and thus offers prediction power for emerging mutant species.

Linkage of posterior amorphous corneal dystrophy to chromosome 12q21.33 and exclusion of coding region mutations in KERA, LUM, DCN, and EPYC.

PURPOSE: To identify the genetic basis of posterior amorphous corneal dystrophy (PACD) segregating in a large pedigree. METHODS: The authors performed clinical evaluation of a previously unreported pedigree with PACD, light and electron microscopic examination of an excised corneal button, genomewide linkage analysis, fine mapping linkage and haplotype analysis, and screening of four candidate genes (KERA, LUM, DCN, and EPYC). RESULTS: Twenty-one participants were determined to be affected based on the presence of characteristic clinical features of PACD; 15 affected and 39 unaffected individuals from a single pedigree enrolled in the study and provided DNA for analysis. Histopathologic examination of an excised corneal specimen from an affected individual demonstrated disorganized stromal lamellae and stromal staining with colloidal iron. Genomewide analysis demonstrated significant evidence of linkage to chromosome region 12q21.33 and evidence suggestive of linkage to chromosome region 8q22.3. Fine mapping of the chromosome 12 locus confirmed significant linkage; the largest multipoint log odds ratio score was 5.6 at D12S351. The linkage support interval was approximately 3.5 Mb (3.5 cM) in length between flanking markers D12S1812 and D12S95, roughly the entire chromosome band 12q21.33. No coding region mutations were identified in four candidate genes-KERA, LUM, DCN, EPYC-located in the chromosome 12 linkage support interval. CONCLUSIONS: Linkage and haplotype analyses identified 12q21.33 as a locus for PACD. However, no mutations were identified in the candidate genes (KERA, LUM, DCN, EPYC) within this region.

The value of avian genomics to the conservation of wildlife.

BACKGROUND: Genomic studies in non-domestic avian models, such as the California condor and white-throated sparrow, can lead to more comprehensive conservation plans and provide clues for understanding mechanisms affecting genetic variation, adaptation and evolution.Developing genomic tools and resources including genomic libraries and a genetic map of the California condor is a prerequisite for identification of candidate loci for a heritable embryonic lethal condition. The white-throated sparrow exhibits a stable genetic polymorphism (i.e. chromosomal rearrangements) associated with variation in morphology, physiology, and behavior (e.g., aggression, social behavior, sexual behavior, parental care).In this paper we outline the utility of these species as well as report on recent advances in the study of their genomes. RESULTS: Genotyping of the condor resource population at 17 microsatellite loci provided a better assessment of the current population's genetic variation. Specific New World vulture repeats were found in the condor genome. Using condor BAC library and clones, chicken-condor comparative maps were generated. A condor fibroblast cell line transcriptome was characterized using the 454 sequencing technology.Our karyotypic analyses of the sparrow in combination with other studies indicate that the rearrangements in both chromosomes 2m and 3a are complex and likely involve multiple inversions, interchromosomal linkage, and pleiotropy. At least a portion of the rearrangement in chromosome 2m existed in the common ancestor of the four North American species of Zonotrichia, but not in the one South American species, and that the 2m form, originally thought to be the derived condition, might actually be the ancestral one. CONCLUSION: Mining and characterization of candidate loci in the California condor using molecular genetic and genomic techniques as well as linkage and comparative genomic mapping will eventually enable the identification of carriers of the chondrodystrophy allele, resulting in improved genetic management of this disease.In the white-throated sparrow, genomic studies, combined with ecological data, will help elucidate the basis of genic selection in a natural population. Morphs of the sparrow provide us with a unique opportunity to study intraspecific genomic differences, which have resulted from two separate yet linked evolutionary trajectories. Such results can transform our understanding of evolutionary and conservation biology.

Associations between single nucleotide polymorphisms in double-stranded DNA repair pathway genes and familial breast cancer.

PURPOSE: DNA damage recognition and repair play a major role in risk for breast cancer. We investigated 104 single nucleotide polymorphisms (SNP) in 17 genes whose protein products are involved in double-stranded break repair (DSBR). EXPERIMENTAL DESIGN: We used a case-control design. Both the case individuals affected with breast cancer or with both breast and ovarian cancers and the controls had similar familial risk of breast cancer and were participants in a high-risk cancer registry. RESULTS: We found that 12 of the polymorphisms are associated with breast or breast and ovarian cancers, most notably rs16888927, rs16888997, and rs16889040, found in introns of RAD21, suggesting that SNPs in other genes in the DSBR pathway in addition to BRCA1 and BRCA2 may affect breast cancer risk. CONCLUSIONS: SNPs within or near several DSBR DNA repair pathway genes are associated with breast cancer in individuals from a high-risk population. In addition, our study reemphasizes the unique perspective that recruitment of cases and controls from family cancer registries has for gene discovery studies.

High-resolution functional profiling of a gammaherpesvirus RTA locus in the context of the viral genome.

Gammaherpesviruses Kaposi's sarcoma-associated herpesvirus and Epstein-Barr virus are associated with multiple human cancers. Our goal was to develop a quantitative, high-throughput functional profiling system to identify viral cis-elements and protein subdomains critical for virus replication in the context of the herpesvirus genome. In gamma-2 herpesviruses, the transactivating factor RTA is essential for initiation of lytic gene expression and viral reactivation. We used the RTA locus as a model to develop the functional profiling approach. The mutant murine gammaherpesvirus 68 viral library, containing 15-bp random insertions in the RTA locus, was passaged in murine fibroblast cells for multiple rounds of selection. The effect of each 15-bp insertion was characterized using fluorescent-PCR profiling. We identified 1,229 insertions in the 3,845-bp RTA locus, of which 393, 282, and 554 were critically impaired, attenuated, and tolerated, respectively, for viral growth. The functional profiling phenotypes were verified by examining several individual RTA mutant clones for transactivating function of the RTA promoter and transcomplementing function of the RTA-null virus. Thus, the profiling approach enabled us to identify several novel functional domains in the RTA locus in the context of the herpesvirus genome. Importantly, our study has demonstrated a novel system to conduct high-density functional genetic mapping. The genome-scale expansion of the genetic profiling approach will expedite the functional genomics research on herpesvirus.

CCL3 genotype and current depression increase risk of HIV-associated dementia.

BACKGROUND: The prevalence of Human immunodeficiency virus (HIV)-associated dementia (HAD) has continued to rise even as incidence has fallen. Several host genetic variants have been identified that modify risk for HAD. However, the findings have not been replicated consistently and most studies did not consider the multitude of factors that might themselves confer risk for HAD. In the current study, we sought to replicate the findings of previous studies in a neurologically and behaviorally well-characterized cohort. METHODS: The sample consisted of 143 HIV+ individuals enrolled in the National NeuroAIDS Tissue Consortium (NNTC). Based on consensus diagnosis, 117 were considered neurologically normal upon study entry, and 26 had HAD. Seven single-nucleotide polymorphisms (SNPs) were genotyped within seven genes (CCL2, CCL3, CCL5, interleukin-1α [IL-1α], IL-10, stromal cell-derived factor 1, and tumor necrosis factor-α). Logistic regression analysis was used to predict group membership (normal vs HAD), with predictor variables including length of infection, age, current drug dependence, current depression, and genotype. RESULTS: The two groups were statistically similar with regards to demographic characteristics, current drug use, and disease factors. The HAD group had significantly greater number of individuals with current depression. Only one SNP, rs1130371 within the gene for CCL3, was entered into the analysis as the others showed symmetric distribution between groups. Logistic regression indicated that current depression and CCL3 genotype were significant predictors of HAD. Depression conferred a fivefold greater risk of HAD, while the TT genotype for CCL3 SNP (rs1130371) was associated with twofold risk for HAD. CONCLUSION: Depression and CCL3 genotype predicted HAD. The fact that SNPs previously found to be associated with HAD were not in our analysis, and that rs1130371 is in high linkage disequilibrium with neighboring genes indicates that more dense genotyping in significantly larger cohorts is required to further characterize the relationship between genotype and risk for HAD.

Quantitative association between HER-2/neu and steroid hormone receptors in hormone receptor-positive primary breast cancer.

BACKGROUND: HER-2/neu, which encodes a receptor tyrosine kinase, is amplified and overexpressed in 20%-25% of human breast cancers. Such tumors are often resistant to hormone therapy. Despite a general inverse association between HER-2/neu amplification/overexpression and estrogen receptor (ER) and/or progesterone receptor (PR) expression, a fraction of patients are both HER-2/neu- and hormone receptor (HR)-positive. The efficacy of hormone therapy in this group is currently a matter of debate. To better understand the relationship between HER-2/neu positivity and HR expression, we analyzed HER-2/neu, ER, and PR as continuous variables in breast cancer cell lines and two cohorts of primary breast cancer patients. METHODS: HER-2/neu and ER/PR expression was analyzed by enzyme-linked immunosorbent assay (ELISA) and enzyme immunoassay (EIA), respectively, in 14 human breast cancer cell lines, some of which had been transfected with the HER-2/neu gene. For the clinical study population, HER-2/neu protein levels were assessed by ELISA (cohort A, n = 665), and HER-2/neu gene copy number was determined using fluorescence in situ hybridization (cohort B, n = 894). ER/PR expression was analyzed by EIA (cohort A) or radioligand binding (cohort B). Associations between HER-2/neu and ER/PR expression were analyzed using Spearman's rho correlation and the chi-square test, and absolute levels were compared using the Mann-Whitney U test. All statistical tests were two-sided. RESULTS: HR-positive human breast cancer cell lines transfected with the HER-2/neu gene expressed statistically significantly lower levels of ER and PR than parental lines. In the clinical cohorts, levels of HER-2/neu overexpression and gene amplification were inversely correlated with ER/PR levels (Cohort A [n = 112]: for ER, r = -0.34, P<.001; for PR, r = -0.24, P =.010. Cohort B [n = 188]: for ER, r = -0.39, P<.001; for PR, r = -0.26, P<.001). Among patients with HR-positive tumors, HER-2/neu-positive tumors had statistically significantly lower ER/PR levels than HER-2/neu-negative ones (Cohort A: for ER, median = 25 fmol/mg [interquartile range [IQR] = 13-78] versus median = 38.5 fmol/mg [IQR = 17-99] and P =.031; for PR, median = 35 fmol/mg [IQR = 12-119] versus median = 88.5 fmol/mg [IQR = 22-236] and P<.001. Cohort B: for ER, median = 44 fmol/mg [IQR = 13-156] versus median = 92 fmol/mg [IQR = 35-235] and P<.001; for PR, median = 36 fmol/mg [IQR = 13-108] versus median = 84 fmol/mg [IQR = 24-250] and P<.001). Patients with higher levels of HER-2/neu overexpression or amplification had statistically significantly lower levels of ER/PR than patients with lower levels of HER-2/neu overexpression or amplification. CONCLUSION: Because absolute HR levels are strongly related to response to hormone therapy in primary and advanced breast cancer, reduced ER/PR expression may be one mechanism to explain the relative resistance of HER-2/neu-positive:HR-positive tumors to hormone therapy.