Data from: Genotyping-by-Sequencing for Populus Population Genomics: An Assessment of Genome Sampling Patterns and Filtering Approaches

Continuing advances in nucleotide sequencing technology are inspiring a suite of genomic approaches in studies of natural populations. Researchers are faced with data management and analytical scales that are increasing by orders of magnitude. With such dramatic advances comes a need to understand biases and error rates, which can be propagated and magnified in large-scale data acquisition and processing. Here we assess genomic sampling biases and the effects of various population-level data filtering strategies in a genotyping-by-sequencing (GBS) protocol. We focus on data from two species of Populus, because this genus has a relatively small genome and is emerging as a target for population genomic studies. We estimate the proportions and patterns of genomic sampling by examining the Populus trichocarpa genome (Nisqually-1), and demonstrate a pronounced bias towards coding regions when using the methylation-sensitive ApeKI restriction enzyme in this species. Using population-level data from a closely related species (P. tremuloides), we also investigate various approaches for filtering GBS data to retain high-depth, informative SNPs that can be used for population genetic analyses. We find a data filter that includes the designation of ambiguous alleles resulted in metrics of population structure and Hardy-Weinberg equilibrium that were most consistent with previous studies of the same populations based on other genetic markers. Analyses of the filtered data (27,910 SNPs) also resulted in patterns of heterozygosity and population structure similar to a previous study using microsatellites. Our application demonstrates that technically and analytically simple approaches can readily be developed for population genomics of natural populations

Recovery Plan for Scots Pine Blister Rust Caused by Cronartium pini

Peer reviewedPostprin

Population structure and hybridization under contemporary and future climates in a heteroploid foundational shrub species (Artemisia tridentata)

Current and past climatic changes can shift plant climatic niches, which may cause spatial overlap or separation between related taxa. The former often leads to hybridization and introgression, which may generate novel variation and influence the adaptive capacity of plants. An additional mechanism facilitating adaptations to novel environments and an important evolutionary driver in plants is polyploidy as the result of whole genome duplication. Artemisia tridentata (big sagebrush) is a landscape-dominating foundational shrub in the western United States which occupies distinct ecological niches, exhibiting diploid and tetraploid cytotypes. Tetraploids have a large impact on the species’ landscape dominance as they occupy a preponderance of the arid spectrum of A. tridentata range. Three distinct subspecies are recognized, which co-occur in ecotones – the transition zone between two or more distinct ecological niches – allowing for hybridization and introgression. Here we assess the genomic distinctiveness and extent of hybridization among subspecies at different ploidies under both contemporary and predicted future climates. We sampled five transects throughout the western United States where a subspecies overlap was predicted using subspecies-specific climate niche models. Along each transect, we sampled multiple plots representing the parental and the potential hybrid habitats. We performed reduced representation sequencing and processed the data using a ploidy-informed genotyping approach. Population genomic analyses revealed distinct diploid subspecies and at least two distinct tetraploid gene pools, indicating independent origins of the tetraploid populations. We detected low levels of hybridization (2.5%) between the diploid subspecies, while we found evidence for increased admixture between ploidy levels (18%), indicating hybridization has an important role in the formation of tetraploids. Our analyses highlight the importance of subspecies co-occurrence within these ecotones to maintain gene exchange and potential formation of tetraploid populations. Genomic confirmations of subspecies in the ecotones support the subspecies overlap predicted by the contemporary climate niche models. However, future mid-century projections of subspecies niches predict a substantial loss in range and subspecies overlap. Thus, reductions in hybridization potential could affect new recruitment of genetically variable tetraploids that are vital to this species’ ecological role. Our results underscore the importance of ecotone conservation and restoration

A haploid pseudo-chromosome genome assembly for a keystone sagebrush species of western North American rangelands

Increased ecological disturbances, species invasions, and climate change are creating severe conservation problems for several plant species that are widespread and foundational. Understanding the genetic diversity of these species and how it relates to adaptation to these stressors are necessary for guiding conservation and restoration efforts. This need is particularly acute for big sagebrush (Artemisia tridentata; Asteraceae), which was once the dominant shrub over 1,000,000 km2 in western North America but has since retracted by half and thus has become the target of one of the largest restoration seeding efforts globally. Here, we present the first reference-quality genome assembly for an ecologically important subspecies of big sagebrush (A. tridentata subsp. tridentata) based on short and long reads, as well as chromatin proximity ligation data analyzed using the HiRise pipeline. The final 4.2-Gb assembly consists of 5,492 scaffolds, with nine pseudo-chromosomal scaffolds (nine scaffolds comprising at least 90% of the assembled genome; n = 9). The assembly contains an estimated 43,377 genes based on ab initio gene discovery and transcriptional data analyzed using the MAKER pipeline, with 91.37% of BUSCOs being completely assembled. The final assembly was highly repetitive, with repeat elements comprising 77.99% of the genome, making the Artemisia tridentata subsp. tridentata genome one of the most highly repetitive plant genomes to be sequenced and assembled. This genome assembly advances studies on plant adaptation to drought and heat stress and provides a valuable tool for future genomic research.This research was made possible by 2 NSF Idaho EPSCoR grants (award numbers OIA-1757324 and OIA-1826801), as well as a Dovetail Genomics Tree of Life Award.Introduction Materials and methods Sample collection, in vitro tissue propagation, and biomass production Flow cytometry and genome complexity analysis PacBio and Omni-C sequence data generation PacBio long-read de novo assembly and validation Pseudomolecule construction with HiRise Genome annotation RNA sequencing Repeat identification Functional annotation Results and discussion Validation of genome assembly and annotation Genome complexity and evidence of past polyploidization Comparing the A. tridentata and A. annua genome assemblies Applications of the sagebrush reference genome Data availability Acknowledgments Literature cite

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although theMYConcogene has been implicated incancer, a systematic assessment of alterations ofMYC, related transcription factors, and co-regulatoryproteins, forming the proximal MYC network (PMN),across human cancers is lacking. Using computa-tional approaches, we define genomic and proteo-mic features associated with MYC and the PMNacross the 33 cancers of The Cancer Genome Atlas.Pan-cancer, 28% of all samples had at least one ofthe MYC paralogs amplified. In contrast, the MYCantagonists MGA and MNT were the most frequentlymutated or deleted members, proposing a roleas tumor suppressors.MYCalterations were mutu-ally exclusive withPIK3CA,PTEN,APC,orBRAFalterations, suggesting that MYC is a distinct onco-genic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such asimmune response and growth factor signaling; chro-matin, translation, and DNA replication/repair wereconserved pan-cancer. This analysis reveals insightsinto MYC biology and is a reference for biomarkersand therapeutics for cancers with alterations ofMYC or the PMN

Transcriptome characterization and polymorphism detection between subspecies of big sagebrush (Artemisia tridentata)

<p>Abstract</p> <p>Background</p> <p>Big sagebrush (<it>Artemisia tridentata</it>) is one of the most widely distributed and ecologically important shrub species in western North America. This species serves as a critical habitat and food resource for many animals and invertebrates. Habitat loss due to a combination of disturbances followed by establishment of invasive plant species is a serious threat to big sagebrush ecosystem sustainability. Lack of genomic data has limited our understanding of the evolutionary history and ecological adaptation in this species. Here, we report on the sequencing of expressed sequence tags (ESTs) and detection of single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers in subspecies of big sagebrush.</p> <p>Results</p> <p>cDNA of <it>A. tridentata </it>sspp. <it>tridentata </it>and <it>vaseyana </it>were normalized and sequenced using the 454 GS FLX Titanium pyrosequencing technology. Assembly of the reads resulted in 20,357 contig consensus sequences in ssp. <it>tridentata </it>and 20,250 contigs in ssp. <it>vaseyana</it>. A BLASTx search against the non-redundant (NR) protein database using 29,541 consensus sequences obtained from a combined assembly resulted in 21,436 sequences with significant blast alignments (≤ 1e^-15). A total of 20,952 SNPs and 119 polymorphic SSRs were detected between the two subspecies. SNPs were validated through various methods including sequence capture. Validation of SNPs in different individuals uncovered a high level of nucleotide variation in EST sequences. EST sequences of a third, tetraploid subspecies (ssp. <it>wyomingensis</it>) obtained by Illumina sequencing were mapped to the consensus sequences of the combined 454 EST assembly. Approximately one-third of the SNPs between sspp. <it>tridentata </it>and <it>vaseyana </it>identified in the combined assembly were also polymorphic within the two geographically distant ssp. <it>wyomingensis </it>samples.</p> <p>Conclusion</p> <p>We have produced a large EST dataset for <it>Artemisia tridentata</it>, which contains a large sample of the big sagebrush leaf transcriptome. SNP mapping among the three subspecies suggest the origin of ssp. <it>wyomingensis </it>via mixed ancestry. A large number of SNP and SSR markers provide the foundation for future research to address questions in big sagebrush evolution, ecological genetics, and conservation using genomic approaches.</p

Menstrual function among women exposed to polybrominated biphenyls: A follow-up prevalence study

BACKGROUND: Alteration in menstrual cycle function is suggested among rhesus monkeys and humans exposed to polybrominated biphenyls (PBBs) and structurally similar polychlorinated biphenyls (PCBs). The feedback system for menstrual cycle function potentially allows multiple pathways for disruption directly through the hypothalamic-pituitary-ovarian axis and indirectly through alternative neuroendocrine axes. METHODS: The Michigan Female Health Study was conducted during 1997–1998 among women in a cohort exposed to PBBs in 1973. This study included 337 women with self-reported menstrual cycles of 20–35 days (age range: 24–56 years). Current PBB levels were estimated by exponential decay modeling of serum PBB levels collected from 1976–1987 during enrollment in the Michigan PBB cohort. Linear regression models for menstrual cycle length and the logarithm of bleed length used estimated current PBB exposure or enrollment PBB exposure categorized in tertiles, and for the upper decile. All models were adjusted for serum PCB levels, age, body mass index, history of at least 10% weight loss in the past year, physical activity, smoking, education, and household income. RESULTS: Higher levels of physical activity were associated with shorter bleed length, and increasing age was associated with shorter cycle length. Although no overall association was found between PBB exposure and menstrual cycle characteristics, a significant interaction between PBB exposures with past year weight loss was found. Longer bleed length and shorter cycle length were associated with higher PBB exposure among women with past year weight loss. CONCLUSION: This study suggests that PBB exposure may impact ovarian function as indicated by menstrual cycle length and bleed length. However, these associations were found among the small number of women with recent weight loss suggesting either a chance finding or that mobilization of PBBs from lipid stores may be important. These results should be replicated with larger numbers of women exposed to similar lipophilic compounds

Postoperative outcomes in oesophagectomy with trainee involvement

BACKGROUND: The complexity of oesophageal surgery and the significant risk of morbidity necessitates that oesophagectomy is predominantly performed by a consultant surgeon, or a senior trainee under their supervision. The aim of this study was to determine the impact of trainee involvement in oesophagectomy on postoperative outcomes in an international multicentre setting. METHODS: Data from the multicentre Oesophago-Gastric Anastomosis Study Group (OGAA) cohort study were analysed, which comprised prospectively collected data from patients undergoing oesophagectomy for oesophageal cancer between April 2018 and December 2018. Procedures were grouped by the level of trainee involvement, and univariable and multivariable analyses were performed to compare patient outcomes across groups. RESULTS: Of 2232 oesophagectomies from 137 centres in 41 countries, trainees were involved in 29.1 per cent of them (n = 650), performing only the abdominal phase in 230, only the chest and/or neck phases in 130, and all phases in 315 procedures. For procedures with a chest anastomosis, those with trainee involvement had similar 90-day mortality, complication and reoperation rates to consultant-performed oesophagectomies (P = 0.451, P = 0.318, and P = 0.382, respectively), while anastomotic leak rates were significantly lower in the trainee groups (P = 0.030). Procedures with a neck anastomosis had equivalent complication, anastomotic leak, and reoperation rates (P = 0.150, P = 0.430, and P = 0.632, respectively) in trainee-involved versus consultant-performed oesophagectomies, with significantly lower 90-day mortality in the trainee groups (P = 0.005). CONCLUSION: Trainee involvement was not found to be associated with significantly inferior postoperative outcomes for selected patients undergoing oesophagectomy. The results support continued supervised trainee involvement in oesophageal cancer surgery

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although the MYC oncogene has been implicated in cancer, a systematic assessment of alterations of MYC, related transcription factors, and co-regulatory proteins, forming the proximal MYC network (PMN), across human cancers is lacking. Using computational approaches, we define genomic and proteomic features associated with MYC and the PMN across the 33 cancers of The Cancer Genome Atlas. Pan-cancer, 28% of all samples had at least one of the MYC paralogs amplified. In contrast, the MYC antagonists MGA and MNT were the most frequently mutated or deleted members, proposing a role as tumor suppressors. MYC alterations were mutually exclusive with PIK3CA, PTEN, APC, or BRAF alterations, suggesting that MYC is a distinct oncogenic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such as immune response and growth factor signaling; chromatin, translation, and DNA replication/repair were conserved pan-cancer. This analysis reveals insights into MYC biology and is a reference for biomarkers and therapeutics for cancers with alterations of MYC or the PMN. We present a computational study determining the frequency and extent of alterations of the MYC network across the 33 human cancers of TCGA. These data, together with MYC, positively correlated pathways as well as mutually exclusive cancer genes, will be a resource for understanding MYC-driven cancers and designing of therapeutics

Search for eccentric black hole coalescences during the third observing run of LIGO and Virgo

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass M&gt;70 M⊙) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities 0&lt;e≤0.3 at 0.33 Gpc−3 yr−1 at 90\% confidence level