Integration of a Cross-Ancestry Polygenic Model With Clinical Risk Factors Improves Breast Cancer Risk Stratification.

PURPOSE: To develop and validate a cross-ancestry integrated risk score (caIRS) that combines a cross-ancestry polygenic risk score (caPRS) with a clinical estimator for breast cancer (BC) risk. We hypothesized that the caIRS is a better predictor of BC risk than clinical risk factors across diverse ancestry groups. METHODS: We used diverse retrospective cohort data with longitudinal follow-up to develop a caPRS and integrate it with the Tyrer-Cuzick (T-C) clinical model. We tested the association between the caIRS and BC risk in two validation cohorts including > 130,000 women. We compared model discrimination for 5-year and remaining lifetime BC risk between the caIRS and T-C and assessed how the caIRS would affect screening in the clinic. RESULTS: The caIRS outperformed T-C alone for all populations tested in both validation cohorts and contributed significantly to risk prediction beyond T-C. The area under the receiver operating characteristic curve improved from 0.57 to 0.65, and the odds ratio per standard deviation increased from 1.35 (95% CI, 1.27 to 1.43) to 1.79 (95% CI, 1.70 to 1.88) in validation cohort 1 with similar improvements observed in validation cohort 2. We observed the largest gain in positive predictive value using the caIRS in Black/African American women across both validation cohorts, with an approximately two-fold increase and an equivalent negative predictive value as the T-C. In a multivariate, age-adjusted logistic regression model including both caIRS and T-C, caIRS remained significant, indicating that caIRS provides information over T-C alone. CONCLUSION: Adding a caPRS to the T-C model improves BC risk stratification for women of multiple ancestries, which could have implications for screening recommendations and prevention.

Whole-genome risk prediction of common diseases in human preimplantation embryos.

Preimplantation genetic testing (PGT) of in-vitro-fertilized embryos has been proposed as a method to reduce transmission of common disease; however, more comprehensive embryo genetic assessment, combining the effects of common variants and rare variants, remains unavailable. Here, we used a combination of molecular and statistical techniques to reliably infer inherited genome sequence in 110 embryos and model susceptibility across 12 common conditions. We observed a genotype accuracy of 99.0-99.4% at sites relevant to polygenic risk scoring in cases from day-5 embryo biopsies and 97.2-99.1% in cases from day-3 embryo biopsies. Combining rare variants with polygenic risk score (PRS) magnifies predicted differences across sibling embryos. For example, in a couple with a pathogenic BRCA1 variant, we predicted a 15-fold difference in odds ratio (OR) across siblings when combining versus a 4.5-fold or 3-fold difference with BRCA1 or PRS alone. Our findings may inform the discussion of utility and implementation of genome-based PGT in clinical practice.

Interspecific Gene Flow and the Evolution of Specialization in Black and White Rhinoceros.

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros.

External Reinfection of a Fungal Pathogen Does not Contribute to Pathogen Growth.

Chytridiomycosis is an emerging infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), which has led to devastating declines in amphibian populations worldwide. Current theory predicts that Bd infections are maintained through both reproduction on the host's skin and reinfection from sources outside of the host. To investigate the importance of external reinfection on pathogen burden, we infected captive-bred individuals of the highly susceptible Panamanian Golden Frog, Atelopus glyphus, and wild-caught glass frogs, Espadarana prosoblepon, with Bd. We housed the animals in one of three treatments: individually, in heterospecific pairs, and in conspecific pairs. For 8 weeks, we measured the Bd load and shedding rate of all frogs. We found that Atelopus had high rates of increase in both Bd load and shedding rate, but pathogen growth rates did not differ among treatments. The infection intensity of Espadarana co-housed with Atelopus was indistinguishable from those housed singly and those in conspecific pairs, despite being exposed to a large external source of Bd zoospores. Our results indicate that Bd load in both species is driven by pathogen replication within an individual, with reinfection from outside the host contributing little to the amplification of host fungal load.

Evaluating recovery potential of the northern white rhinoceros from cryopreserved somatic cells.

The critically endangered northern white rhinoceros is believed to be extinct in the wild, with the recent death of the last male leaving only two remaining individuals in captivity. Its extinction would appear inevitable, but the development of advanced cell and reproductive technologies such as cloning by nuclear transfer and the artificial production of gametes via stem cells differentiation offer a second chance for its survival. In this work, we analyzed genome-wide levels of genetic diversity, inbreeding, population history, and demography of the white rhinoceros sequenced from cryopreserved somatic cells, with the goal of informing how genetically valuable individuals could be used in future efforts toward the genetic rescue of the northern white rhinoceros. We present the first sequenced genomes of the northern white rhinoceros, which show relatively high levels of heterozygosity and an average genetic divergence of 0.1% compared with the southern subspecies. The two white rhinoceros subspecies appear to be closely related, with low genetic admixture and a divergent time <80,000 yr ago. Inbreeding, as measured by runs of homozygosity, appears slightly higher in the southern than the northern white rhinoceros. This work demonstrates the value of the northern white rhinoceros cryopreserved genetic material as a potential gene pool for saving this subspecies from extinction.

Dead or alive? Viability of chytrid zoospores shed from live amphibian hosts.

Pathogens vary in virulence and rates of transmission because of many differences in the host, the pathogen, and their environment. The amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), affects amphibian hosts differently, causing extinction and population declines in some species but having limited effects on others. Phenotypic differences in zoospore production rates among Bd lineages likely contribute to some of the variation observed among host responses, although no studies have quantified the viability of zoospores shed from live animals. We compared host survivorship, infection intensity, shedding rates, and zoospore viability between 2 species of endangered tropical frogs, Hylomantis lemur and Atelopus zeteki, when exposed to a highly virulent lineage of Bd (JEL 423). We applied a dye to zoospores 30 to 60 min following animal soaks, to estimate shedding rate and proportion of live zoospores shed by different species. The average infection intensity for A. zeteki was nearly 17 times higher (31,455 ± 10,103 zoospore genomic equivalents [ZGEs]) than that of H. lemur (1832 ± 1086 ZGEs), and A. zeteki died earlier than H. lemur. The proportion of viable zoospores was ~80% in both species throughout the experiment, although A. zeteki produced many more zoospores, suggesting it may play a disproportionate role in spreading disease in communities where it occurs, because the large number of viable zoospores they produce might increase infection in other species where they are reintroduced.

Isotopic incorporation rates and discrimination factors in mantis shrimp crustaceans.

Stable isotope analysis has provided insights into the trophic ecology of a wide diversity of animals. Knowledge about isotopic incorporation rates and isotopic discrimination between the consumer and its diet for different tissue types is essential for interpreting stable isotope data, but these parameters remain understudied in many animal taxa and particularly in aquatic invertebrates. We performed a 292-day diet shift experiment on 92 individuals of the predatory mantis shrimp, Neogonodactylus bredini, to quantify carbon and nitrogen incorporation rates and isotope discrimination factors in muscle and hemolymph tissues. Average isotopic discrimination factors between mantis shrimp muscle and the new diet were 3.0 ± 0.6 ‰ and 0.9 ± 0.3 ‰ for carbon and nitrogen, respectively, which is contrary to what is seen in many other animals (e.g. C and N discrimination is generally 0-1 ‰ and 3-4 ‰, respectively). Surprisingly, the average residence time of nitrogen in hemolymph (28.9 ± 8.3 days) was over 8 times longer than that of carbon (3.4 ± 1.4 days). In muscle, the average residence times of carbon and nitrogen were of the same magnitude (89.3 ± 44.4 and 72.8 ± 18.8 days, respectively). We compared the mantis shrimps' incorporation rates, along with rates from four other invertebrate taxa from the literature, to those predicted by an allometric equation relating carbon incorporation rate to body mass that was developed for teleost fishes and sharks. The rate of carbon incorporation into muscle was consistent with rates predicted by this equation. Our findings provide new insight into isotopic discrimination factors and incorporation rates in invertebrates with the former showing a different trend than what is commonly observed in other animals.

More than skin deep: functional genomic basis for resistance to amphibian chytridiomycosis.

The amphibian-killing chytrid fungus Batrachochytrium dendrobatidis (Bd) is one of the most generalist pathogens known, capable of infecting hundreds of species globally and causing widespread population declines and extinctions. However, some host species are seemingly unaffected by Bd, tolerating or clearing infections without clinical signs of disease. Variation in host immune responses is commonly evoked for these resistant or tolerant species, yet to date, we have no direct comparison of amphibian species responses to infection at the level of gene expression. In this study, we challenged four Central American frog species that vary in Bd susceptibility, with a sympatric virulent strain of the pathogen. We compared skin and spleen orthologous gene expression using differential expression tests and coexpression gene network analyses. We found that resistant species have reduced skin inflammatory responses and increased expression of genes involved in skin integrity. In contrast, only highly susceptible species exhibited suppression of splenic T-cell genes. We conclude that resistance to chytridiomycosis may be related to a species' ability to escape the immunosuppressive activity of the fungus. Moreover, our results indicate that within-species differences in splenic proteolytic enzyme gene expression may contribute to intraspecific variation in survival. This first comparison of amphibian functional immunogenomic architecture in response to Bd provides insights into key genetic mechanisms underlying variation in disease outcomes among amphibian species.

A fungal pathogen of amphibians, Batrachochytrium dendrobatidis, attenuates in pathogenicity with in vitro passages.

Laboratory investigations into the amphibian chytrid fungus, Batrachochytrium dendrobatidis (Bd), have accelerated recently, given the pathogen's role in causing the global decline and extinction of amphibians. Studies in which host animals were exposed to Bd have largely assumed that lab-maintained pathogen cultures retained the infective and pathogenic properties of wild isolates. Attenuated pathogenicity is common in artificially maintained cultures of other pathogenic fungi, but to date, it is unknown whether, and to what degree, Bd might change in culture. We compared zoospore production over time in two samples of a single Bd isolate having different passage histories: one maintained in artificial media for more than six years (JEL427-P39), and one recently thawed from cryopreserved stock (JEL427-P9). In a common garden experiment, we then exposed two different amphibian species, Eleutherodactylus coqui and Atelopus zeteki, to both cultures to test whether Bd attenuates in pathogenicity with in vitro passages. The culture with the shorter passage history, JEL427-P9, had significantly greater zoospore densities over time compared to JEL427-P39. This difference in zoospore production was associated with a difference in pathogenicity for a susceptible amphibian species, indicating that fecundity may be an important virulence factor for Bd. In the 130-day experiment, Atelopus zeteki frogs exposed to the JEL427-P9 culture experienced higher average infection intensity and 100% mortality, compared with 60% mortality for frogs exposed to JEL427-P39. This effect was not observed with Eleutherodactylus coqui, which was able to clear infection. We hypothesize that the differences in phenotypic performance observed with Atelopus zeteki are rooted in changes of the Bd genome. Future investigations enabled by this study will focus on the underlying mechanisms of Bd pathogenicity.

Pathophysiology in mountain yellow-legged frogs (Rana muscosa) during a chytridiomycosis outbreak.

The disease chytridiomycosis is responsible for declines and extirpations of amphibians worldwide. Chytridiomycosis is caused by a fungal pathogen (Batrachochytrium dendrobatidis) that infects amphibian skin. Although we have a basic understanding of the pathophysiology from laboratory experiments, many mechanistic details remain unresolved and it is unknown if disease development is similar in wild amphibian populations. To gain a better understanding of chytridiomycosis pathophysiology in wild amphibian populations, we collected blood biochemistry measurements during an outbreak in mountain yellow-legged frogs (Rana muscosa) in the Sierra Nevada Mountains of California. We found that pathogen load is associated with disruptions in fluid and electrolyte balance, yet is not associated with fluctuations acid-base balance. These findings enhance our knowledge of the pathophysiology of this disease and indicate that disease development is consistent across multiple species and in both laboratory and natural conditions. We recommend integrating an understanding of chytridiomycosis pathophysiology with mitigation practices to improve amphibian conservation.

Is chytridiomycosis an emerging infectious disease in Asia?

The disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has caused dramatic amphibian population declines and extinctions in Australia, Central and North America, and Europe. Bd is associated with >200 species extinctions of amphibians, but not all species that become infected are susceptible to the disease. Specifically, Bd has rapidly emerged in some areas of the world, such as in Australia, USA, and throughout Central and South America, causing population and species collapse. The mechanism behind the rapid global emergence of the disease is poorly understood, in part due to an incomplete picture of the global distribution of Bd. At present, there is a considerable amount of geographic bias in survey effort for Bd, with Asia being the most neglected continent. To date, Bd surveys have been published for few Asian countries, and infected amphibians have been reported only from Indonesia, South Korea, China and Japan. Thus far, there have been no substantiated reports of enigmatic or suspected disease-caused population declines of the kind that has been attributed to Bd in other areas. In order to gain a more detailed picture of the distribution of Bd in Asia, we undertook a widespread, opportunistic survey of over 3,000 amphibians for Bd throughout Asia and adjoining Papua New Guinea. Survey sites spanned 15 countries, approximately 36° latitude, 111° longitude, and over 2000 m in elevation. Bd prevalence was very low throughout our survey area (2.35% overall) and infected animals were not clumped as would be expected in epizootic events. This suggests that Bd is either newly emerging in Asia, endemic at low prevalence, or that some other ecological factor is preventing Bd from fully invading Asian amphibians. The current observed pattern in Asia differs from that in many other parts of the world.

Individual and combined effects of multiple pathogens on Pacific treefrogs.

In nature, individual hosts often encounter multiple pathogens simultaneously, which can lead to additive, antagonistic, or synergistic effects on hosts. Synergistic effects on infection prevalence or severity could greatly affect host populations. However, ecologists and managers often overlook the influence of pathogen combinations on hosts. This is especially true in amphibian conservation, even though multiple pathogens coexist within amphibian populations, and several pathogens have been implicated in amphibian population declines and extinctions. Using an amphibian host, Pseudacris regilla (Pacific treefrog), we experimentally investigated interactive effects among three pathogens: the trematode Ribeiroia sp. (hereafter, Ribeiroia), the fungus Batrachochytrium dendrobatidis (hereafter, BD), and the water mold Achlya flagellata. We detected no effects of A. flagellata, but did find effects of Ribeiroia and BD that varied depending on context. Low doses of Ribeiroia caused relatively few malformations, while higher Ribeiroia doses caused numerous deformities dominated by missing and reduced limbs and limb elements. Exposure to low doses of BD accelerated larval host development, despite there being no detectable BD infections, while exposure to higher BD doses caused infection but did not alter developmental rate. Hosts exposed to both Ribeiroia and BD exhibited the highest mortality, although overall evidence of interactive effects of multiple pathogens was limited. We suggest further research on the influence of multi-pathogen assemblages on amphibians, particularly under a variety of ecological conditions and with a wider diversity of hosts and pathogens.

Dynamics of an emerging disease drive large-scale amphibian population extinctions.

Epidemiological theory generally suggests that pathogens will not cause host extinctions because the pathogen should fade out when the host population is driven below some threshold density. An emerging infectious disease, chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd) is directly linked to the recent extinction or serious decline of hundreds of amphibian species. Despite continued spread of this pathogen into uninfected areas, the dynamics of the host-pathogen interaction remain unknown. We use fine-scale spatiotemporal data to describe (i) the invasion and spread of Bd through three lake basins, each containing multiple populations of the mountain yellow-legged frog, and (ii) the accompanying host-pathogen dynamics. Despite intensive sampling, Bd was not detected on frogs in study basins until just before epidemics began. Following Bd arrival in a basin, the disease spread to neighboring populations at approximately 700 m/yr in a wave-like pattern until all populations were infected. Within a population, infection prevalence rapidly reached 100% and infection intensity on individual frogs increased in parallel. Frog mass mortality began only when infection intensity reached a critical threshold and repeatedly led to extinction of populations. Our results indicate that the high growth rate and virulence of Bd allow the near-simultaneous infection and buildup of high infection intensities in all host individuals; subsequent host population crashes therefore occur before Bd is limited by density-dependent factors. Preventing infection intensities in host populations from reaching this threshold could provide an effective strategy to avoid the extinction of susceptible amphibian species in the wild.

Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis.

Global amphibian decline by chytridiomycosis is a major environmental disaster that has been attributed to either recent fungal spread or environmental change that promotes disease. Here, we present a population genetic comparison of Batrachochytrium dendrobatidis isolates from an intensively studied region of frog decline, the Sierra Nevada of California. In support of a novel pathogen, we find low diversity, no amphibian-host specificity, little correlation between fungal genotype and geography, local frog extirpation by a single fungal genotype, and evidence of human-assisted fungus migration. In support of endemism, at a local scale, we find some diverse, recombining populations. Therefore neither epidemic spread nor endemism alone explains this particular amphibian decline. Recombination raises the possibility of resistant sporangia and a mechanism for rapid spread as well as persistence that could greatly complicate global control of the pathogen.