cLD: Rare-variant linkage disequilibrium between genomic regions identifies novel genomic interactions.

Linkage disequilibrium (LD) is a fundamental concept in genetics; critical for studying genetic associations and molecular evolution. However, LD measurements are only reliable for common genetic variants, leaving low-frequency variants unanalyzed. In this work, we introduce cumulative LD (cLD), a stable statistic that captures the rare-variant LD between genetic regions, which reflects more biological interactions between variants, in addition to lack of recombination. We derived the theoretical variance of cLD using delta methods to demonstrate its higher stability than LD for rare variants. This property is also verified by bootstrapped simulations using real data. In application, we find cLD reveals an increased genetic association between genes in 3D chromatin interactions, a phenomenon recently reported negatively by calculating standard LD between common variants. Additionally, we show that cLD is higher between gene pairs reported in interaction databases, identifies unreported protein-protein interactions, and reveals interacting genes distinguishing case/control samples in association studies.

Analyses of allele age and fitness impact reveal human beneficial alleles to be older than neutral controls.

A classic population genetic prediction is that alleles experiencing directional selection should swiftly traverse allele frequency space, leaving detectable reductions in genetic variation in linked regions. However, despite this expectation, identifying clear footprints of beneficial allele passage has proven to be surprisingly challenging. We addressed the basic premise underlying this expectation by estimating the ages of large numbers of beneficial and deleterious alleles in a human population genomic data set. Deleterious alleles were found to be young, on average, given their allele frequency. However, beneficial alleles were older on average than non-coding, non-regulatory alleles of the same frequency. This finding is not consistent with directional selection and instead indicates some type of balancing selection. Among derived beneficial alleles, those fixed in the population show higher local recombination rates than those still segregating, consistent with a model in which new beneficial alleles experience an initial period of balancing selection due to linkage disequilibrium with deleterious recessive alleles. Alleles that ultimately fix following a period of balancing selection will leave a modest 'soft' sweep impact on the local variation, consistent with the overall paucity of species-wide 'hard' sweeps in human genomes.

Diverse African genomes reveal selection on ancient modern human introgressions in Neanderthals.

Comparisons of Neanderthal genomes to anatomically modern human (AMH) genomes show a history of Neanderthal-to-AMH introgression stemming from interbreeding after the migration of AMHs from Africa to Eurasia. All non-sub-Saharan African AMHs have genomic regions genetically similar to Neanderthals that descend from this introgression. Regions of the genome with Neanderthal similarities have also been identified in sub-Saharan African populations, but their origins have been unclear. To better understand how these regions are distributed across sub-Saharan Africa, the source of their origin, and what their distribution within the genome tells us about early AMH and Neanderthal evolution, we analyzed a dataset of high-coverage, whole-genome sequences from 180 individuals from 12 diverse sub-Saharan African populations. In sub-Saharan African populations with non-sub-Saharan African ancestry, as much as 1% of their genomes can be attributed to Neanderthal sequence introduced by recent migration, and subsequent admixture, of AMH populations originating from the Levant and North Africa. However, most Neanderthal homologous regions in sub-Saharan African populations originate from migration of AMH populations from Africa to Eurasia ∼250 kya, and subsequent admixture with Neanderthals, resulting in ∼6% AMH ancestry in Neanderthals. These results indicate that there have been multiple migration events of AMHs out of Africa and that Neanderthal and AMH gene flow has been bi-directional. Observing that genomic regions where AMHs show a depletion of Neanderthal introgression are also regions where Neanderthal genomes show a depletion of AMH introgression points to deleterious interactions between introgressed variants and background genomes in both groups-a hallmark of incipient speciation.

Genetic architecture and evolution of color variation in American black bears.

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Immediate Lymphatic Reconstruction during Axillary Node Dissection for Breast Cancer: A Systematic Review and Meta-analysis.

The objective of this study is to summarize the current body of evidence detailing the impact of immediate lymphatic reconstruction (ILR) on the incidence of breast cancer-related lymphedema (BCRL) following axillary node dissection (ALND). Methods: Medline and Embase databases were queried for publications, where ILR was performed at the time of ALND for breast cancer. Exclusion criteria included lymphaticovenous anastomosis for established BCRL, animal studies, non-breast cancer patient population studies, and descriptive studies detailing surgical technique. Meta-analysis was performed with a forest plot generated using a Mantel -Haenszel statistical method, with a random-effect analysis model. Effect measure was reported as risk ratios with associated 95% confidence intervals. The risk of bias within studies was assessed by the Cochrane Collaboration tool. Results: This systematic review yielded data from 11 studies and 417 breast cancer patients who underwent ILR surgery at the time of ALND. There were 24 of 417 (5.7%) patients who developed BCRL following ILR. Meta-analysis revealed that in the ILR group, 6 of 90 patients (6.7%) developed lymphedema, whereas in the control group, 17 of 50 patients (34%) developed lymphedema. Patients in the ILR group had a risk ratio of 0.22 (CI, 0.09 -0.52) of lymphedema with a number needed to treat of four. Conclusions: There is a clear signal indicating the benefit of ILR in preventing BCRL. Randomized control trials are underway to validate these findings. ILR may prove to be a beneficial intervention for improving the quality of life of breast cancer survivors.

Peripheral nerve surgical competency in plastic surgery and neurosurgery residents.

OBJECTIVE: To assess the case volume and self-perceived competence of current mandatory skills in peripheral nerve surgery. DESIGN: Cross sectional survey based study examining case volume and self-reported competence in peripheral nerve surgery. SETTING: Canadian Neurosurgery and Plastic Surgery accredited residency programs PARTICIPANTS: All Canadian Neurosurgery and Plastic Surgery senior trainees (PGY 3+) invited to participate RESULTS: Much variability exists in both exposure to cases and perceived senior resident competence for both plastic and neurosurgery residents. Confidence in surgical ability as perceived competency is lower in trainees for more advanced peripheral nerve procedures. Self- reported confidence increased with post-graduate experience. CONCLUSIONS: Overall, the findings in this study highlight the importance of increasing operative experience in complex peripheral nerve surgery among surgical residents.

The importance of including ethnically diverse populations in studies of quantitative trait evolution.

For many traits, human variation is less a matter of categorical differences than quantitative variation, such as height, where individuals fall along a continuum from short to tall. Most recent studies utilize large population-based samples with whole-genome sequences to study the evolution of these traits and have made significant progress implementing a broad spectrum of techniques. However, relatively few studies of quantitative trait evolution include ethnically diverse populations, which often harbor the highest levels of genetic and phenotypic diversity. Thus, our ability to draw inferences about quantitative trait adaptation has been limited. Here, we review recent studies examining human quantitative trait adaptation, and argue that including ethnically diverse populations, particularly from Africa, will be especially informative for our understanding of how humans adapt to the world around them.

Beare-Stevenson Syndrome With Blepharoptosis as a Complication of Front-Orbital Advancement and Remodeling.

Beare-Stevenson syndrome (BSS) is an extremely rare genetic disorder characterized by a broad range of congenital malformations including craniosynostosis, cutis gyrata, facial deformities, and abnormal genitalia. The authors report a case of a 7 month old female who developed a mechanical ptosis secondary to dermatochalasis as a complication of fronto-orbital advancement and remodeling (FOAR) surgery which subsequently required multiple lid surgeries to reverse ptosis. This is the first report of blepharoptosis correction in a child with BSS as a complication of FOAR.

An estimator of first coalescent time reveals selection on young variants and large heterogeneity in rare allele ages among human populations.

Allele age has long been a focus of population genetic research, primarily because it can be an important clue to the fitness effects of an allele. By virtue of their effects on fitness, alleles under directional selection are expected to be younger than neutral alleles of the same frequency. We developed a new coalescent-based estimator of a close proxy for allele age, the time when a copy of an allele first shares common ancestry with other chromosomes in a sample not carrying that allele. The estimator performs well, including for the very rarest of alleles that occur just once in a sample, with a bias that is typically negative. The estimator is mostly insensitive to population demography and to factors that can arise in population genomic pipelines, including the statistical phasing of chromosomes. Applications to 1000 Genomes Data and UK10K genome data confirm predictions that singleton alleles that alter proteins are significantly younger than those that do not, with a greater difference in the larger UK10K dataset, as expected. The 1000 Genomes populations varied markedly in their distributions for singleton allele ages, suggesting that these distributions can be used to inform models of demographic history, including recent events that are only revealed by their impacts on the ages of very rare alleles.

OCMA: Fast, Memory-Efficient Factorization of Prohibitively Large Relationship Matrices.

Matrices representing genetic relatedness among individuals (i.e., Genomic Relationship Matrices, GRMs) play a central role in genetic analysis. The eigen-decomposition of GRMs (or its alternative that generates fewer top singular values using genotype matrices) is a necessary step for many analyses including estimation of SNP-heritability, Principal Component Analysis (PCA), and genomic prediction. However, the GRMs and genotype matrices provided by modern biobanks are too large to be stored in active memory. To accommodate the current and future "bigger-data", we develop a disk-based tool, Out-of-Core Matrices Analyzer (OCMA), using state-of-the-art computational techniques that can nimbly perform eigen and Singular Value Decomposition (SVD) analyses. By integrating memory mapping (mmap) and the latest matrix factorization libraries, our tool is fast and memory-efficient. To demonstrate the impressive performance of OCMA, we test it on a personal computer. For full eigen-decomposition, it solves an ordinary GRM (N = 10,000) in 55 sec. For SVD, a commonly used faster alternative of full eigen-decomposition in genomic analyses, OCMA solves the top 200 singular values (SVs) in half an hour, top 2,000 SVs in 0.95 hr, and all 5,000 SVs in 1.77 hr based on a very large genotype matrix (N = 1,000,000, M = 5,000) on the same personal computer. OCMA also supports multi-threading when running in a desktop or HPC cluster. Our OCMA tool can thus alleviate the computing bottleneck of classical analyses on large genomic matrices, and make it possible to scale up current and emerging analytical methods to big genomics data using lightweight computing resources.

Adaptive Landscape of Protein Variation in Human Exomes.

The human genome contains hundreds of thousands of missense mutations. However, only a handful of these variants are known to be adaptive, which implies that adaptation through protein sequence change is an extremely rare phenomenon in human evolution. Alternatively, existing methods may lack the power to pinpoint adaptive variation. We have developed and applied an Evolutionary Probability Approach (EPA) to discover candidate adaptive polymorphisms (CAPs) through the discordance between allelic evolutionary probabilities and their observed frequencies in human populations. EPA reveals thousands of missense CAPs, which suggest that a large number of previously optimal alleles experienced a reversal of fortune in the human lineage. We explored nonadaptive mechanisms to explain CAPs, including the effects of demography, mutation rate variability, and negative and positive selective pressures in modern humans. Many nonadaptive hypotheses were tested, but failed to explain the data, which suggests that a large proportion of CAP alleles have increased in frequency due to beneficial selection. This suggestion is supported by the fact that a vast majority of adaptive missense variants discovered previously in humans are CAPs, and hundreds of CAP alleles are protective in genotype-phenotype association data. Our integrated phylogenomic and population genetic EPA approach predicts the existence of thousands of nonneutral candidate variants in the human proteome. We expect this collection to be enriched in beneficial variation. The EPA approach can be applied to discover candidate adaptive variation in any protein, population, or species for which allele frequency data and reliable multispecies alignments are available.

Protein evolution depends on multiple distinct population size parameters.

That population size affects the fate of new mutations arising in genomes, modulating both how frequently they arise and how efficiently natural selection is able to filter them, is well established. It is therefore clear that these distinct roles for population size that characterize different processes should affect the evolution of proteins and need to be carefully defined. Empirical evidence is consistent with a role for demography in influencing protein evolution, supporting the idea that functional constraints alone do not determine the composition of coding sequences.Given that the relationship between population size, mutant fitness and fixation probability has been well characterized, estimating fitness from observed substitutions is well within reach with well-formulated models. Molecular evolution research has, therefore, increasingly begun to leverage concepts from population genetics to quantify the selective effects associated with different classes of mutation. However, in order for this type of analysis to provide meaningful information about the intra- and inter-specific evolution of coding sequences, a clear definition of concepts of population size, what they influence, and how they are best parameterized is essential.Here, we present an overview of the many distinct concepts that "population size" and "effective population size" may refer to, what they represent for studying proteins, and how this knowledge can be harnessed to produce better specified models of protein evolution.

Metallic intraocular foreign body as detected by magnetic resonance imaging without complications- A case report.

PURPOSE: To describe a case and present unique images of a metallic intraocular foreign body that was identified in a 12-year-old male patient who underwent routine magnetic resonance imaging (MRI) to assess neurodevelopmental delay. OBSERVATIONS: We present MRI and diagnostic imaging of a metallic intraocular foreign body in a young patient with no known history of trauma or reason for the existence of metal in the eye area. Computed tomography scan was performed to confirm the presence of the intraocular foreign body, followed by optical coherence tomography and electroretinogram to assess visual status. It was determined that no surgical intervention was currently required as no visual impairment or ocular toxicity was identified. The patient continues to be monitored. CONCLUSIONS AND IMPORTANCE: This case presentation highlights the novel imaging features of a metallic intraocular foreign body, unexpectedly detected with MRI.

Low-dose ranibizumab as primary treatment of posterior type I retinopathy of prematurity.

OBJECTIVE: The aim of this study was to evaluate the visual and structural outcomes of eyes that received ranibizumab as treatment for retinopathy of prematurity (ROP). METHODS: This was a retrospective case series of infants who received a 0.2 mg (0.02 mL) intravitreal injection of ranibizumab as the primary treatment for type 1 ROP. Outcome measures included regression or recurrence of ROP, complications of treatment, and assessment of visual acuity and refractive error. RESULTS: Forty-two eyes of 21 infants (13 male) were included. Mean gestational age and birth weight were 24.6 ± 1.3 weeks and 613 ± 91 g, respectively. Mean age at injection was 37.4 ± 2.2 weeks postmenstrual age, and mean follow-up period was 10.1 ± 7 months. Active neovascularization regressed rapidly, and anatomical outcomes were favourable in all eyes. Twelve eyes of 6 infants received supplemental laser photocoagulation at a mean post-menstrual age (PMA) of 72.0 ± 27.3 weeks when vascularization had not advanced beyond zone II. Visual acuity was measurable in 28 of 42 eyes. Mean visual acuity was 0.94 ± 0.36 logMAR. Mean spherical equivalent was +1.00. There were no ocular or systemic complications in these patients and no cicatricial complications were observed with no progression to stage 4 or 5 disease. CONCLUSIONS: A single intravitreal dose of 0.2 mg (0.02 mL) ranibizumab showed favourable anatomical and functional outcomes in eyes with type 1 ROP.

Genome-wide signature of local adaptation linked to variable CpG methylation in oak populations.

It has long been known that adaptive evolution can occur through genetic mutations in DNA sequence, but it is unclear whether adaptive evolution can occur through analogous epigenetic mechanisms, such as through DNA methylation. If epigenetic variation contributes directly to evolution, species under threat of disease, invasive competition, climate change or other stresses would have greater stores of variation from which to draw. We looked for evidence of natural selection acting on variably methylated DNA sites using population genomic analysis across three climatologically distinct populations of valley oaks. We found patterns of genetic and epigenetic differentiations that indicate local adaptation is operating on large portions of the oak genome. While CHG methyl polymorphisms are not playing a significant role and would make poor targets for natural selection, our findings suggest that CpG methyl polymorphisms as a whole are involved in local adaptation, either directly or through linkage to regions under selection.

Coselected genes determine adaptive variation in herbivore resistance throughout the native range of Arabidopsis thaliana.

The "mustard oil bomb" is a major defense mechanism in the Brassicaceae, which includes crops such as canola and the model plant Arabidopsis thaliana. These plants produce and store blends of amino acid-derived secondary metabolites called glucosinolates. Upon tissue rupture by natural enemies, the myrosinase enzyme hydrolyses glucosinolates, releasing defense molecules. Brassicaceae display extensive variation in the mixture of glucosinolates that they produce. To investigate the genetics underlying natural variation in glucosinolate profiles, we conducted a large genome-wide association study of 22 methionine-derived glucosinolates using A. thaliana accessions from across Europe. We found that 36% of among accession variation in overall glucosinolate profile was explained by genetic differentiation at only three known loci from the glucosinolate pathway. Glucosinolate-related SNPs were up to 490-fold enriched in the extreme tail of the genome-wide [Formula: see text] scan, indicating strong selection on loci controlling this pathway. Glucosinolate profiles displayed a striking longitudinal gradient with alkenyl and hydroxyalkenyl glucosinolates enriched in the West. We detected a significant contribution of glucosinolate loci toward general herbivore resistance and lifetime fitness in common garden experiments conducted in France, where accessions are enriched in hydroxyalkenyls. In addition to demonstrating the adaptive value of glucosinolate profile variation, we also detected long-distance linkage disequilibrium at two underlying loci, GS-OH and GS-ELONG. Locally cooccurring alleles at these loci display epistatic effects on herbivore resistance and fitness in ecologically realistic conditions. Together, our results suggest that natural selection has favored a locally adaptive configuration of physically unlinked loci in Western Europe.

Spatial localization of recent ancestors for admixed individuals.

Ancestry analysis from genetic data plays a critical role in studies of human disease and evolution. Recent work has introduced explicit models for the geographic distribution of genetic variation and has shown that such explicit models yield superior accuracy in ancestry inference over nonmodel-based methods. Here we extend such work to introduce a method that models admixture between ancestors from multiple sources across a geographic continuum. We devise efficient algorithms based on hidden Markov models to localize on a map the recent ancestors (e.g., grandparents) of admixed individuals, joint with assigning ancestry at each locus in the genome. We validate our methods by using empirical data from individuals with mixed European ancestry from the Population Reference Sample study and show that our approach is able to localize their recent ancestors within an average of 470 km of the reported locations of their grandparents. Furthermore, simulations from real Population Reference Sample genotype data show that our method attains high accuracy in localizing recent ancestors of admixed individuals in Europe (an average of 550 km from their true location for localization of two ancestries in Europe, four generations ago). We explore the limits of ancestry localization under our approach and find that performance decreases as the number of distinct ancestries and generations since admixture increases. Finally, we build a map of expected localization accuracy across admixed individuals according to the location of origin within Europe of their ancestors.

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation.

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, 'antagonistic' pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits 'adaptive' pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a 'drought escape' strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.

A new era of human population genetics.

The 1000 Genomes Project Consortium has recently published an important early contribution to a new generation of systematic surveys of rare human genetic variation.

A mixed-model approach for genome-wide association studies of correlated traits in structured populations.

Genome-wide association studies (GWAS) are a standard approach for studying the genetics of natural variation. A major concern in GWAS is the need to account for the complicated dependence structure of the data, both between loci as well as between individuals. Mixed models have emerged as a general and flexible approach for correcting for population structure in GWAS. Here, we extend this linear mixed-model approach to carry out GWAS of correlated phenotypes, deriving a fully parameterized multi-trait mixed model (MTMM) that considers both the within-trait and between-trait variance components simultaneously for multiple traits. We apply this to data from a human cohort for correlated blood lipid traits from the Northern Finland Birth Cohort 1966 and show greatly increased power to detect pleiotropic loci that affect more than one blood lipid trait. We also apply this approach to an Arabidopsis thaliana data set for flowering measurements in two different locations, identifying loci whose effect depends on the environment.