The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise.

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

High prevalence of Lynx rufus gammaherpesvirus 1 in wild Vermont bobcats.

Gammaherpesviruses (GHVs) are host specific DNA viruses that infect a large range of mammalian species. These viruses preferentially target host lymphocyte cell populations and infection may lead to morbidity or mortality in immunocompromised, co-infected, or non-adapted hosts. In this study, we tested for the presence of Lynx rufus gammaherpesvirus 1 (LruGHV1) in a northeastern United States population of wild bobcats (L. rufus). We estimated prevalence of infection and viral load in infected individuals using quantitative real-time PCR analysis of spleen DNA from 64 Vermont bobcats. We observed an overall prevalence of 64% using this methodology. Bobcat age was significantly positively associated with GHV infection status, and we noted a trend for higher viral loads in young animals, but prevalence and viral load were similar in male and female bobcats. A single LruGHV1 variant was identified from the sequencing of the viral glycoprotein B gene of Vermont bobcats. This gene sequence was 100% similar to that reported in Florida bobcats and slightly variant from other isolates identified in the Western USA. Our work suggests broad geographic distribution and high prevalence of LruGHV1 in bobcat populations across the United States with infection attributes that suggest horizontal transmission of the agent. Geographic differences in viral genotype may reflect historical migration and expansion events among bobcat populations.

Genome-wide association study of lung function phenotypes in a founder population.

BACKGROUND: Lung function is a long-term predictor of mortality and morbidity. OBJECTIVE: We sought to identify single nucleotide polymorphisms (SNPs) associated with lung function. METHODS: We performed a genome-wide association study (GWAS) of FEV1, forced vital capacity (FVC), and FEV1/FVC in 1144 Hutterites aged 6 to 89 years, who are members of a founder population of European descent. We performed least absolute shrinkage and selection operation regression to select the minimum set of SNPs that best predict FEV1/FVC in the Hutterites and used the GRAIL algorithm to mine the Gene Ontology database for evidence of functional connections between genes near the predictive SNPs. RESULTS: Our GWAS identified significant associations between FEV1/FVC and SNPs at the THSD4-UACA-TLE3 locus on chromosome 15q23 (P = 5.7 × 10(-8) to 3.4 × 10(-9)). Nine SNPs at or near 4 additional loci had P < 10(-5) with FEV1/FVC. Only 2 SNPs were found with P < 10(-5) for FEV1 or FVC. We found nominal levels of significance with SNPs at 9 of the 27 previously reported loci associated with lung function measures. Among a predictive set of 80 SNPs, 6 loci were identified that had a significant degree of functional connectivity (GRAIL P < .05), including 3 clusters of ß-defensin genes, 2 chemokine genes (CCL18 and CXCL12), and TNFRSF13B. CONCLUSION: This study identifies genome-wide significant associations and replicates results of previous GWASs. Multimarker modeling implicated for the first time common variation in genes involved in antimicrobial immunity in airway mucosa that influences lung function.

Human leukocyte antigen-G expression in differentiated human airway epithelial cells: lack of modulation by Th2-associated cytokines.

BACKGROUND: Human leukocyte antigen (HLA)-G is a nonclassical class I antigen with immunomodulatory roles including up-regulation of suppressor T regulatory lymphocytes. HLA-G was recently identified as an asthma susceptibility gene, and expression of a soluble isoform, HLA-G5, has been demonstrated in human airway epithelium. Increased presence of HLA-G5 has been demonstrated in bronchoalveolar lavage fluid recovered from patients with mild asthma; this suggests a role for this isoform in modulating airway inflammation though the mechanisms by which this occurs is unclear. Airway inflammation associated with Th2 cytokines such as IL-4 and IL-13 is a principal feature of asthma, but whether these cytokines elicit expression of HLA-G is not known. METHODS: We examined gene and protein expression of both soluble (G5) and membrane-bound (G1) HLA-G isoforms in primary differentiated human airway epithelial cells collected from normal lungs and grown in air-liquid interface culture. Cells were treated with up to 10 ng/ml of either IL-4, IL-5, or IL-13, or 100 ng/ml of the immunomodulatory cytokine IL-10, or 10,000 U/ml of the Th1-associated cytokine interferon-beta, for 24 hr, after which RNA was isolated for evaluation by quantitative PCR and protein was collected for Western blot analysis. RESULTS: HLA-G5 but not G1 was present in dAEC as demonstrated by quantitative PCR, western blot and confocal microscopy. Neither G5 nor G1 expression was increased by the Th2-associated cytokines IL-4, IL-5 or IL-13 over 24 hr, nor after treatment with IL-10, but was increased 4.5 ± 1.4 fold after treatment with 10,000 U/ml interferon-beta. CONCLUSIONS: These data demonstrate the constitutive expression of a T lymphocyte regulatory molecule in differentiated human airway epithelial cells that is not modulated by Th2-associated cytokines.

The maternal HLA-G 1597ΔC null mutation is associated with increased risk of pre-eclampsia and reduced HLA-G expression during pregnancy in African-American women.

The non-classical major histocompatibility complex molecule, human leukocyte antigen (HLA)-G, is thought to contribute to maternal immune tolerance and successful placentation during pregnancy. Genetic polymorphisms in HLA-G are known to influence expression levels as well as the relative expression of individual protein isoforms. As diminished or aberrant HLA-G expression patterns may contribute to the development of certain pregnancy complications, we sought to investigate the association between functional HLA-G polymorphisms and the risk of pre-eclampsia (PE) in African-American women. The association between maternal and fetal genotype at six HLA-G polymorphisms and risk of PE was assessed in 372 pregnancies (314 normotensive; 58 pre-eclamptic). We observed an elevated risk of PE (P = 0.00027) in pregnancies where the mother carried the 1597ΔC allele, a null allele that abolishes expression of full-length HLA-G isoforms. Furthermore, the frequency of the maternal 1597ΔC allele was highest in the subset of pre-eclamptic pregnancies that were delivered preterm, suggesting an association between the null allele and the severity of PE. We then replicated the association between higher maternal 1597ΔC allele frequency and increased severity of PE (P = 0.038) in an independent sample of 533 African-American women. Finally, to investigate the mechanistic basis of this association, we measured circulating soluble HLA-G (sHLA-G) concentrations in maternal serum collected during pregnancy in 51 healthy, normotensive African-American control women and found significantly lower levels in women carrying the 1597ΔC allele (P = 0.012). These results demonstrate that maternal HLA-G genotype is significantly associated with risk of PE in African-American women and is predictive of circulating sHLA-G levels during pregnancy.

HLA-G polymorphisms and soluble HLA-G protein levels in women with recurrent pregnancy loss from Basrah province in Iraq.

HLA-G is a nonclassical, class I major histocompatibility complex (MHC) gene that exhibits immunomodulatory properties and likely plays a role in the maintenance of successful pregnancy. In this study, we investigated the role of HLA-G polymorphisms on risk for recurrent pregnancy loss (RPL) and on circulating levels of soluble (s)HLA-G in Iraqi women. DNA and plasma were obtained from blood samples collected at 9-12 weeks gestation from 50 women with RPL and 50 healthy pregnant women in Basrah province, Iraq. As measured by ELISA, median sHLA-G levels were significantly lower in the RPL cases compared to healthy controls (21.4 vs. 38.8 U/ml, respectively; P=0.025), and decreased with increasing maternal age (P=0.0051). However, HLA-G allele and haplotype frequencies did not differ significantly between cases and controls (P values ≥0.12 for all tests). In contrast, homozygosity for the C allele (CC) at a tri-allelic promoter polymorphism, -725C/G/T, was associated with lower concentrations of sHLA-G compared to the CG or CT genotypes (median levels 21.1 vs. 40.1 vs. 42.6 U/ml, respectively; P=0.0089). These results demonstrate that HLA-G genotype influences circulating sHLA-G levels during pregnancy but is not significantly associated with risk of RPL.

IFNG genotype and sex interact to influence the risk of childhood asthma.

BACKGROUND: Asthma is a complex disease characterized by sex-specific differences in incidence, prevalence, and severity, but little is known about the molecular basis of these sex-based differences. OBJECTIVE: To investigate the genetic architecture of sex differences in asthma risk, we evaluated (1) associations between polymorphisms in the IFNG gene and childhood-onset asthma in combined and sex-specific samples and (2) interactions between polymorphisms and sex on asthma risk. METHODS: Main and sex-interaction effects of IFNG genetic diversity on asthma risk and IFN-γ levels were examined in a birth cohort of children at high risk for asthma and allergic diseases. Replication of the genetic association was assessed in an independent sample of asthma cases. RESULTS: Significant genotype-sex interactions on asthma were observed for 2 IFNG single nucleotide polymorphisms, rs2069727 and rs2430561, which were in strong linkage disequilibrium with each other. In contrast, none of the 10 IFNG single nucleotide polymorphisms showed significant main effects on asthma. The observed genotype-sex interaction on asthma was characterized by nonadditivity; that is, heterozygous boys had the highest risk for asthma, and heterozygous girls had the lowest risk. The interaction effect was robust to other asthma risk factors but was limited to children who experienced wheezing illnesses with viral infections during the first 3 years of life. Genotype-sex interactions were also observed in the IFN-γ response to LPS in the first year of life. Finally, the sex-interaction effect was replicated in an independent population of childhood asthma cases. CONCLUSIONS: These results provide insight into the genetic basis of sex differences in asthma and highlight the potential importance of interactions among sex, genotype, and environmental factors in asthma pathogenesis.

Sex-specific genetic architecture of human disease.

Sexual dimorphism in anatomical, physiological and behavioural traits are characteristics of many vertebrate species. In humans, sexual dimorphism is also observed in the prevalence, course and severity of many common diseases, including cardiovascular diseases, autoimmune diseases and asthma. Although sex differences in the endocrine and immune systems probably contribute to these observations, recent studies suggest that sex-specific genetic architecture also influences human phenotypes, including reproductive, physiological and disease traits. It is likely that an underlying mechanism is differential gene regulation in males and females, particularly in sex steroid-responsive genes. Genetic studies that ignore sex-specific effects in their design and interpretation could fail to identify a significant proportion of the genes that contribute to risk for complex diseases.

Ancient polymorphism and functional variation in the primate MHC-DQA1 5' cis-regulatory region.

Precise regulation of MHC gene expression is critical to vertebrate immune surveillance and response. Polymorphisms in the 5' proximal promoter region of the human class II gene HLA-DQA1 have been shown to influence its transcriptional regulation and may contribute to the pathogenesis of autoimmune diseases. We investigated the evolutionary history of this cis-regulatory region by sequencing the DQA1 5' proximal promoter region in eight nonhuman primate species. We observed unexpectedly high levels of sequence variation and multiple strong signatures of balancing selection in this region. Specifically, the considerable DQA1 promoter region diversity was characterized by abundant shared (or trans-species) polymorphism and a pronounced lack of fixed differences between species. The majority of transcription factor binding sites in the DQA1 promoter region were polymorphic within species, and these binding site polymorphisms were commonly shared among multiple species despite evidence for negative selection eliminating a significant fraction of binding site mutations. We assessed the functional consequences of intraspecific promoter region diversity using a cell line-based reporter assay and detected significant differences among baboon DQA1 promoter haplotypes in their ability to drive transcription in vitro. The functional differentiation of baboon promoter haplotypes, together with the significant deviations from neutral sequence evolution, suggests a role for balancing selection in the evolution of DQA1 transcriptional regulation in primates.

An imprinted locus epistatically influences Nstr1 and Nstr2 to control resistance to nerve sheath tumors in a neurofibromatosis type 1 mouse model.

Cancer is a complex disease in which cells acquire many genetic and epigenetic alterations. We have examined how three types of alterations, mutations in tumor suppressor genes, changes in an imprinted locus, and polymorphic loci, interact to affect tumor susceptibility in a mouse model of neurofibromatosis type 1 (NF1). Mutations in tumor suppressor genes such as TP53 and in oncogenes such as KRAS have major effects on tumorigenesis due to the central roles of these genes in cell proliferation and cell survival. Imprinted genes expressed from only one parental chromosome affect tumorigenesis if their monoallelic expression is lost or duplicated. Because imprinted loci are within regions deleted or amplified in cancer, the parental origin of genomic rearrangements could affect tumorigenesis. Gene polymorphisms can vary tumor incidence by affecting rate-limiting steps in tumorigenesis within tumor cells or surrounding stroma. In our mouse model of NF1, the incidence of tumors mutant for the tumor suppressor genes Nf1 and Trp53 is strongly modified by a linked imprinted locus acting epistatically on two unlinked polymorphic loci, Nstr1 and Nstr2. This interaction of an imprinted locus and polymorphic susceptibility loci has profound implications for human mapping studies where the parental contribution of alleles is often unknown.

Positive selection on MMP3 regulation has shaped heart disease risk.

BACKGROUND: The evolutionary forces of mutation, natural selection, and genetic drift shape the pattern of phenotypic variation in nature, but the roles of these forces in defining the distributions of particular traits have been hard to disentangle. To better understand the mechanisms contributing to common variation in humans, we investigated the evolutionary history of a functional polymorphism in the upstream regulatory region of the MMP3 gene. This single base pair insertion/deletion variant, which results in a run of either 5 or 6 thymidines 1608 bp from the transcription start site, alters transcription factor binding and influences levels of MMP3 mRNA and protein. The polymorphism contributes to variation in arterial traits and to the risk of coronary heart disease and its progression. RESULTS: Phylogenetic and population genetic analysis of primate sequences indicate that the binding site region is rapidly evolving and has been a hot spot for mutation for tens of millions of years. We also find evidence for the action of positive selection, beginning approximately 24,000 years ago, increasing the frequency of the high-expression allele in Europe but not elsewhere. Positive selection is evident in statistical tests of differentiation among populations and haplotype diversity within populations. Europeans have greater arterial elasticity and suffer dramatically fewer coronary heart disease events than they would have had this selection not occurred. CONCLUSIONS: Locally elevated mutation rates and strong positive selection on a cis-regulatory variant have shaped contemporary phenotypic variation and public health.

Susceptibility to astrocytoma in mice mutant for Nf1 and Trp53 is linked to chromosome 11 and subject to epigenetic effects.

Astrocytoma is the most common malignant brain tumor in humans. Loss of the p53 signaling pathway and up-regulation of the ras signaling pathway are common during tumor progression. We have shown previously that mice mutant for Trp53 and Nf1 develop astrocytoma, progressing to glioblastoma, on a C57BL/6J strain background. In contrast, here we present data that mice mutant for Trp53 and Nf1 on a 129S4/SvJae background are highly resistant to developing astrocytoma. Through analysis of F1 progeny, we demonstrate that susceptibility to astrocytoma is linked to chromosome 11, and that the modifier gene(s) responsible for differences in susceptibility is closely linked to Nf1 and Trp53. Furthermore, this modifier of astrocytoma susceptibility is itself epigenetically modified. These data demonstrate that epigenetic effects can have a strong effect on whether cancer develops in the context of mutant ras signaling and mutant p53, and that this mouse model of astrocytoma can be used to identify modifier phenotypes with complex inheritance patterns that would be unidentifiable in humans. Because analysis of gene function in the mouse is often performed on a mixed C57BL/6,129 strain background, these data also provide a powerful example of the potential of these strains to mask interesting gene functions.